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PREFACE 


The following pages embody, substantially as delivered, a 
course of fifteen lectures given during the Spring of 1925 at the 
Moore School of Electrical Engineering of the University of 
Pennsylvania. 

After a brief introduction to the subject of electric circuit 
theory, the first five chapters are devoted to a systematic and 
fairly complete exposition and critique of the Heaviside Oper¬ 
ational Calculus, a remarkably direct and powerful method 
for the solution of the differential equations of electric circuit 
theory. 

The name of Oliver Heaviside is known to engineers the 
world over; his operational calculus, however, is known to, 
and employed by, only a relatively few specialists, and this 
notwithstanding its remarkable properties and wide applicability 
not only to electric circuit theory but also to the differential 
equations of mathematical physics. In the writer’s. opinion 
this neglect is due less to the intrinsic difficulties of the subject 
than to unfortunate obscurities in Heaviside’s own exposition. 
In the present work the operational calculus is made to depend 
on an integral equation from which the Heaviside Rules and 
Formulas are simply but rigorously deducible. It is the hope 
of the writer that this mode of approach and exposition will be 
of service in securing a wider use of the operational calculus by 
engineers and physicists, and a fuller and more just appreciation 
of its unique advantages. 

The second part of the present work deals with advanced 
problems of electric circuit theory, and in particular with the 
theory of the propagation of current and voltage in electrical 
transmission systems. It is hoped that this part will be of 
interest to electrical engineers generally because, while only a 
few of the results are original with the present work, most of the 
transmission theory dealt with is to be found only in scattered 
memoirs, and there accompanied by formidable mathematical 
difficulties. 



VI 


PREFACE 


While the method of solution employed in the second part is 
largely that of the operational calculus, I have not hesitated to 
employ developments and extensions not to be found in Heavi¬ 
side. For example, the formulation of the problem as a Poisson 
integral equation is an original development which has proved 
quite useful in the actual numerical solution of complicated 
problems. The same may be said of the Chapter on Variable 
Electric Circuit Theory. 

In view of its two-fold aspect this work may therefore be 
regarded either as an exposition and development of the oper¬ 
ational calculus with applications to electric circuit theory, or 
as a contribution to advanced electric circuit theory, depending 
on whether the reader^s viewpoint is that of the mathematician 
or the engineer. 

In’ the text I have not made reference to the literature of the 

c 

subject, now fairly extensive, except where the treatment in 
the text is inadequate and I wish to refer the reader to a fuller 
and more satisfactory discussion. In the bibliography, how¬ 
ever, there is furnished a list of original papers and memoirs. 
While this bibliography includes all pertinent papers known 
to the author, no claim to completeness can be made and it is 
possible that important papers have been omitted. The author 
would be grateful to any reader calling his attention to such 
omissions. 

It remains a pleasurable duty to acknowledge the assistance 
I have received in the preparation of this book. To Dr. Harold 
Pender, Dean of the Moore School of Electrical Engineering of 
the University of Pennsylvania I owe the invitation to deliver 
the lectures out of which this book grew, and stimulating interest 
and advice throughout their course. To my associate Mr. R. S. 
Hoyt I am indebted for a careful reading of the manuscript, 
the verification of most of the formulas, and valuable con¬ 
tributions and advice. For example, the extremely compact 
proof of formula (29) is a modification and improvement due 
to Mr. Hoyt of an earlier proof. To my assistant, Miss Dorothy 
Angell are due all numerical calculations and curves and the 
careful proofreading of the entire book. 

Finally to the Officials of the American Telephone and 
Telegraph Company and in particular to Mr. 0. B. Blackwell, 
Transmission Development Engineer, I would express appre- 
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ciation of their generous attitude toward the research under¬ 
lying this book and my indebtedness for the opportunity and 
encouragement to put this work into shape for publication. 

John E. Carson. 

New York, N. Y. 

October f 1926. 
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ELECTRIC CIRCUIT THEORY 
AND THE OPERATIONAL 
CALCULUS 

CHAPTER I 

THE FUNDAMENTALS OF ELECTRIC CIRCUIT THEORY 

While a knowledge, on the reader^s part, of the elements of 
electric circuit theory will be assumed, it seems well to start 
with a brief review of the fundamental physical principles of 
circuit theory, the mode of formulating the equations, and 
some general theorems which will prove useful subsequently. 

Electric Circuit Theory is that branch of general electro¬ 
magnetic theory which deals with electrical oscillations in linear 
electrical networks. An electrical network is a connected set of 
circuits or meshes each of which can be regarded as made up of 
resistance, inductance, and capacitance elements. If the number 
of independent closed meshes or circuits is finite, the network 
is said to have a finite number of degrees of freedom equal to 
the number of meshes, and the problem is mathematically 
formulated by a finite number of simultaneous linear differ¬ 
ential equations. If the number of closed circuits or meshes 
is made infinite and the size of the elements infinitesimal the 
problem is formulated by a partial differential equation, such 
as the telegraph equation. 

There are two distinct problems involved in electric circuit 
theory. First, the derivation of the differential equations from 
the physical laws and the constants and connections of the 
system, that is, the mathematical formulation of the physical 
problem: and secondly the solution of the differential equations. 

The equations of circuit theory may be established in a number 
of different ways. For example, they may be based on MaxwelPs 
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dynamical theory. In accordance with this method, the net¬ 
work forms a dynamic system in which the currents play the 
role of velocities. If we therefore set up the expressions for the 
kinetic energy, potential energy and dissipation, the network 
equations are deducible from general dynamic equations. 

The simplest, and for our purposes, a quite satisfactory basis 
for the equations of circuit theory are found in Kirchhoff^s 
Laws. These laws state that 

1, The total impressed force taken around any closed loop 
or circuit in the network is equal to the potential drop due to 
(a) resistance, (b) inductive reaction and (c) capacitive reactance. 

2. The sum of the currents entering any branch point in the 
network is always zero. 

Let us now apply these laws to an elementary circuit in order 
to deduce the physical significance of the circuit elements. 

Consider an elementary circuit consisting of a resistance 
element i?, an inductance element L and a capacity element 
C in series, and let an electromotive force E be applied to this 
circuit. If I denote the current in the circuit, the resistance 
drop is BI, the inductance drop is Ldl/dt and the drop across 
the condenser is Q/C where Q is the charge on the condenser. 
It is evident that Q and I are related by the equation I == dQ/dt 
or Q=^ J'Idt. Now apply Kirchhoff^s law relating to the drop 
around the circuit: it gives the equation 

RI+LdI/dt+Q/C = E. 

Multiply both sides by I: we get 

The right hand side is clearly the rate at which the impressed 
force is delivering energy to the circuit, while the left hand side 
is the rate at which energy is. being absorbed by the circuit. 
The first term RP is the rate at which electrical energy is being 
converted into heat. Hence the resistance element may be 
defined as a device for converting electrical energy into heat. 
d 1 

The second term ^ ^ increase of the magnetic 

energy. Hence the inductance element is a device for storing 

d 

energy in the magnetic field. The third term Q^/2C is the 
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rate of increase of the electric energy. Hence the condenser is a 
device for storing energy in the electric field. 

In the foregoing we have isolated and idealized the circuit 
elements. Actually, of course, every circuit element dissipates 
some energy in the form of heat and stores some energy in the 
magnetic field and some in the electric field. The analysis of the 
actual circuit element, however, into three ideal components is 
quite convenient and useful, and should lead to no misconception 
if properly interpreted. 

Now consider the general form of network possessing n inde¬ 
pendent meshes or circuits. Let us number these from 1 to n, 
and let the corresponding mesh currents be denoted by Ji, 
/2 . . . . /n. Let electromotive forces Eij E 2 ... . En 
be applied to the n meshes or circuits respectively. Let L/y, 
Rjj, Cjj denote the total inductance, resistance and capacity in 
series in mesh j and let Rjk, Cjk denote the corresponding 
mutual elements between circuits j and k. Now write down 
Kirchhoff^s equation for any circuit or mesh, say mesh 1; it is 

(Lu 1 +^ 11 +^ 1 + (L^, + 

■ ■ ■ ■ 

Corresponding equations hold for each and every one of the n 
meshes of the network. Writing them all down, we have the 
system of equations 

.( 1 ) 

(inl ^+Rnl + -^-J • • • + J 

The system of simultaneous differential equations (1) 
constitutes the canonical equations of electric circuit theory. 
The interpretation and solution of these equations constitute 
the subject of Electric Circuit Theory, and it is in connection 
>with their solution that we find the most direct and logical 
introduction to the Operational Calculus. 

As an example of the appropriate mode of setting up the 
circuit equations, consider the two mesh network shown in Fig. 
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1’. ‘ Writing down Kirchhoff’s Law for meshes 1 and 2, respec¬ 
tively, we have 


(^Li|+Ea+^-Jd{) 




C^ 

d 

dt 






1/2=402 


In this case the self and mutual coefficients are given by 


Lii^Li 
Cii = Ci 
Rii^Ri 


Ij22 =/^2 
C 22 = C^2 
jK22 = R 2 


Li2=L2\ == +il/ 

Cl2 = C2t = 0 

i ?12 = 1^21 = 0 . 


The conventions adopted for the positive directions of currents 
and voltages are indicated by the arrows. The. sign of the 
mutual inductance M will depend on the relative mode of 
winding of the two coils. 

Now write down Kirchhoff Law, or the circuital equations for 
the network of Fig 2. They are 


Comparison with equations (1) shows that 


Ln—Li-\-Ez 

i?ii= 

1 _1 1 

CxrCi Ca 


■L22 =Li-\-Ls 
R22 “ R^~]~Rz 
1 ^1 _1 
'C 22 


Li2 =Lzi = —La 

i2l2 = i?21 “ — 
1^1 

C 12 C 21 Ca 


It ,should he observed that the signs of the mutual coefficients 
Ria, Eli', <7i2 are a matter of convention. For example if the 
conventional directions of J 2 and E^ are reversed, the signs of 
the mutual coefficients are reversed. 

The system of equations (1) possesses two important properties 
which are largely responsible for the relative simplicity of classical 
electric circuit theory. First, the equations are linear in both 
currents and applied electromotive forces. Secondly, the 
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coefficients L/fc, iJ/t, C^k are constants. Important electro¬ 
technical problems exist, in which these properties no longer 
obtain. The solution, however, for the restricted system of 
linear equations with constant coefficients is fundamental and 
its solution can be extended to important problems involving 
non-linear relations and variable coefficients. These extensions 
will be taken up briefly in a later chapter. 



Another important property is the reciprocal relation among 
the coefficients; that is Ljk-Lki, Rjk = Rkj, and Cjk = Ckj. It is 
easily shown that these reciprocal relations mean that there are 
no concealed sources or sinks of energy. Again important cases 
exist where the reciprocal relations do not hold. Such excep¬ 
tions, however, while of physical interest do not affect the mathe- 


Ri Ct Li 




La C2 R2 
. —oj 



|E2 


Fig. 2. 


matical methods of solution, to which the reciprocal relation is 
not essential. 

Returning to equationTT) we shall now derive the equation of 
activity. Multiply the first equation by Ji, the second by Jo, 
etc. and add: we get 
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The right hand side is the rate at which the applied forces ai 
supplying energy to the network. The first term on the left i 
the rate of increase of the magnetic energy 




huIJk. 


while the second term is the rate of increase of the electric cnerg 



The last term, i»s the rate at which electromagnet; 

energy is being converted into heat in the network. CJonsequentl 
in the electrical network, the magnetic energy is a homogeneoi 
quadratic function of the currents, the electric energy is a hom( 
geneous quadratic function of the charges, and the rate < 
dissipation is a homogeneous quadratic function of tjie current 
In Maxwell’s dynamical theory of electrical networks, the? 
relations were written down at the start and the circuit equatioi 
then derived by an application of Lagrange’s dynamic equatioi 
to the homogeneous quadratic functions. 

Returning to equations (1), we observe that, due to tl 
presence of the integral sign, they arc integro-differential equ 
tions. They are, however, at once reducible to differenti 
equations by the substitution I = dQ/dtj whence they bccon 



Here, as a matter of convenience, we have written l/Cju — H 
It is often more convenient, at least at the outset, to deal wi 
equations (3) rather than (1). 

The Exponential Solution,—In taking u}) the inatlutmatic; 
solution of equations (1), we shall start with the expoiiantt 
solutum. This is of fundanumtal inqiortance, both th(3orotica] 
and practically. It serves as the most direct introduction to i 
Heaviside operational calculus, and in addition furnish(‘s t 
basis of the steady stale solution, or the tlu^ory of alternaii] 
currents. 

To derive this solution we set and put all the oth 

forces E 2 , , . En equal to zero. Tliis latter restriction is 
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mere matter of convenience, and, in virtue of the linear character 
of the equations, involves no loss of generality. 

Now, corresponding to Ei-=^Fie^\ let us assume a solution of 
the form 

(i = l,2, .. n) 

where Jj is a constant. So far this is a pure assumption, and 
its correctness must be verified by substitution in the differential 
equations. 

Now if = it follows at once that 

and 

Now substitute these relations in equations (1) and cancel the 
common factor We then get the system of simultaneous 
equations 

1“!“ ‘ * "b O^Ein--{-Rin-\-‘l/\Cin)J n~F 1, 
(XL21+i?21 + l AC2l)t/1+ • * +(XL2n+i22ri + lAC^27i)»^» — 
. (4) 

(\Lnl~\~Rnl~\~^/^Cnl)J l”l“ * * “b C^nn'^Rnn'^^/^Cjin)J n — 0. 

We note that this is a system of simultaneous algebraic equa¬ 
tions from which the time factor has disappeared. It is this that 
makes the exponential solution so simple, since we can immedi¬ 
ately pass from differential equations to algebraic equations. 
In these algebraic equations, 7i in number, there are n unknown 
quantities Ji, . . Jn. These can therefore all be uniquely 
determined. We thus see that the assumed form of solution is 
possible. 

The notation of equations (4) may be profitably simplified 
as follows: write 

XL jk + R jk +1 /X(7 3k ~ ^ jk (X) = Z jk 

and we have 

J l‘b^l2 J 

221<I^l‘b^22*/2*b ■ ' = /jrX 


ZnlJ 2+ * ‘ +^nn J'n -^0. 
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The solution of this system of cciiiations is 

I)^' 


y JLYX jXyy \i X' J 4 - 

//= ‘I'\cy=ry C^‘ 

jU /jji 

where D is the detonninant of the cf)cfficieiit.s, 


^13. . . * 

- Zu 

Z 21 Z 22 Zu • '' < 

, , . .z^n 

;^31 ^32..-- 

• • • ♦ ^3r* ^ 

Zn\ Zn*l .. • 

. , , ► Zmn 


and Mji is the cofactor, or minor with prnp(!r sign, of tin* jtli 
column and first row. 

I shall not attempt to discuss the theory of detenhin.'inls on 
which this solution is based. ^ We may note, however, one 
important property. Since Zjk = Zk„ Mjk^Mkj. From this the 
Reciprocal Theorem follows immediately. Thi.s may be .«tii(ed 
as follows: 

If a force Fe^‘ is applied in the jth mesh, or bruneij, of the 
network, the current in the /cth mesh, or branch, i.s by the 
foregoing 

j^lC . 


Now apply the same force in the /cth mesh, or brancii, tlam 
current in the jth mesh is 


tins 


Comparing these expressions and remembering that - M,,, 
it follows that the current in the /cth branch corresponding to an 
exponential impressed e.m.f. in the ith branch, is CKiual to the 
current in the jth branch corresponding to the .same e.m.f, in 
the /cth branch. This relation is of the greatest t<iclmicul 
importance. 


In many important technical problems we arc interfisted only 
in two accessible branches, such as the sending and nnadving. 
In such cases, where we are not concerned with the currents in 
the other meshes or branches, it is often convenient to eliminate 


^For a remarkably concise and complete discussion of the exponential 
Mlution by md of the theory of determinants, see Cisoidal Oscillations. 
Trans. A. I. E. E., 1911, by G. A. Campbell. 
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them from the equation. Thus, suppose that we have electro¬ 
motive forces El and E 2 in meshes 1 and 2 and are concerned 
only with the currents in these meshes. If we solve equations 
3,4, . . n, w—2 in number, for Is . . in terms of Ji and Is 
and then substitute in (1) and (2) we get 

Ziih-^-Ziili^E^. ( 8 ) 


The Steady State Solution.—The steady state solution, on 
which the whole theory of alternating currents depends, is 
immediately derivable from the exponential solution. Let us 
suppose that 

E 2 =-Ei= ■ ■ ■ =Sn=0 and that = cos (wi—e). 

Now by virtue of the well known formula in the theory of the 
complex variable, cos we can write 

=|'(cos 0—i sin 0)1^6“*-t-|(cos d-\-i sin 6)Fe~""\ (9) 

Now, by virtue of this formula, the applied electromotive force Ei 
consists of two exponential forces, one varying as e‘“‘ and the 
other as Hence it is easy to see that the currents are 

made up of two components, thus 

0 = 1,2 .. n) (10) 

and we have merely to use the exponential solution given above, 
substituting for X,fw and -iu respectively. That is, 

IF' r//_l F" 

2 Z,-x(-ico) 


or 


1 Fe-'^ i-i I 1 
— ct nr If 


Peio 


Q 


“2Zn{wT ''2Zni-M) 

The second term is the conjugate imaginary of the first, so that 


7,=B 


Fe—^ 




'Zjiiito) 

Z,-i(iw) 

\Zji(t()>)\ 

^ cos io3t — 6—<t>). 


'\Z(m)\ 


We thus arrive at the rule for the steady state solution. 
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If the applied e.m.f. is F cos (coi —0), substitute m for d/dt in 
the differential equations, determine the impedance function 

Z{i<^)^D{ic^)/M{io>) (11) 

by the solution of the algebraic equations, and write it in the 
form 

= ( 12 ) 

Then the required solution is 

This in compact form contains the whole theory of the symbolic 
solution of alternating current problems. 

The Complementary Solution.—So far in the solutions which 
we have discussed the currents are of the same type as the 
impressed forces: that is to say in physical language, the currents 
are “forced currents and vary with time in precisely the same 
manner as do the electromotive forces. Such currents are, 
however, in general only part of the total currents. In addition 
to the forced currents we have also the characteristic oscillations; 
or, in mathematical language, the complete solution must include 
both particular and complementary solutions. This may be 
shown as follows: Let Ji', .... 7/be solutions of the comple¬ 
mentary equations, 

. '( 14 ) 

Then if Ji . . . is a solution of (1), Ji-|-7/, . . . is 

also a solution. 

To derive the solution of the complementary system of equa¬ 
tions (14), assume that a solution exists of the form 
I/ = J/ e ^ (j-=l,2 . . ?^) 

so that d/dt^\ and S di — l/\. Substitute in equations (14) 
and cancel out the common factor e}\ Then we have 


i2fll(X)t7|'+ • * ’ • +2;in(X)/n'=0, 

^^nl(X)Ji'+ .... +^„„(\)7/== 0. 


(15) 


This is a system of n homogeneous equations in the unknown 
quantities J/, . . Jn. The condition that a finite solution 
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shall exist is that, in accordance with a well known principle of 
the theory of equations, the determinant of the coefficients 
shall vanish. That is, 


^11 (X) • • ■ • Zir, (X) 


J5(X) = 


= 0. (16) 


^nl (X) • • • • Znn (X) 


Consequently the possible values of X must be such that this 
equation is satisfied. In other words, X must be a root of the 
equation jD(X) = 0. Let these roots be denoted by Xi, X 2 , . . X,n. 
Then, assigning to X any one of these values, we can determine 
the ratio J//JiJ from any (n —1) of the equations. That is to 
say, if we take 

= (17) 

substitution in any (n — 1) of the equations determines J 2 ', . . U- 
The m constants Ci, . . Cm are so far, however, entirely arbitrary, 
and are at our disposal to satisfy imposed boundary conditions. 

This introduces us to the idea of boundary conditions which 
is of the greatest importance in circuit theory. In physical 
language the boundary conditions denote the state of the system 
when the electromotive force is applied or when any change in 
the circuit constants occurs. The number of independent 
boundary conditions which can, in general, be satisfied is equal 
to the number of roots of the equation Z>(X) = 0. Evidently, 
therefore, it is physically impossible to impose more boundary 
conditions than this. On the other hand, if this number of 
boundary conditions is not specified, the complete solution is 
indeterminate: that is to say, the problem is not correctly set. 
As an example of boundary conditions, we may specify that the 
electromotive force is applied at time ;f==0, and that at this time 
all the currents in the inductances and all the charges on the 
condensers are 55ei*o. 

So far we have been following the classical theory of linear 
differential equations. We have seen that the forced exponential 
solution and the derived steady state solution are extremely 
simple and are mere matters of elementary algebra. The 
practical difficulties in the classical method of solution begin 
with the determination of the constants Ci, . . Cm oi the 
complementary solution as well as the roots Xi, . . X^ of the 
equation D(X)==0. It is at this point that Heaviside broke 
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with classical methods, and by considering special boundary 
conditions of great physical importance, and particular types of 
impressed forces, laid the foundations of original and powerful 
methods of solution. We shall therefore at this point follow 
Heaviside’s example and attack the problem from a different 
standpoint. In doing this we shall not at once take up an 
exposition of Heaviside’s own method of attack. We shall 
first establish some fundamental theorems which are extremely 
powerful and will serve us as a guide in interpreting and ration¬ 
alizing the Heaviside operational calculus. 



CHAPTER II 


THE SOLUTION WHEN AN ARBITRARY FORCE IS 
APPLIED TO THE NETWORK IN A STATE OF 
V EQUILIBRIUM 

In engineering applications of electric circuit theory there are 
four outstanding problems: 

(1) The steady state distribution of currents and potentials 
when the network is energized by a sinusoidal electromotive 
force. This problem is the subject of the theory of alternating 
currents v^ich forms^he basis of our calculations of power lines 
and the more elaborate networks of communication systems. 

(2) The distribution of currents and potentials in the network 
in response to an arbitrary electromotive force applied to the 
network in a state of equilibrium, i.e., applied when the currents 
and charges in the network are identically zero. 

(3) The effect on the distribution of currents and potentials of 
suddenly changing a circuit constant or connection, such as 
opening or closing a switch, while the system is energized. 

(4) The free oscillations of the network from an arbitary initial 
configuration. This problem will not be considered in detail 
for several reasons. In the first place it is of minor technical 
importance. Furthermore, if the initial state is producible by 
electromotive forces applied at accessible terminals of the net¬ 
work, the problem is reducible to and can be dealt with as a 
special case of problem (3). Finally, by reason of the arbitrary 
initial boundary conditions, the formulation and solution of the 
problem is inherently incapable of reduction to the elegant sim¬ 
plicity of those problems directly amenable to the operational 
calculus. The formal solution of problem (4) is, however, 
straightforward and is dealt with in treatises on dynamics. 
After evaluating the roots of, the equation D (X) = 0 (see last 
paragraph of Chap. I), the constants of integration Ci, . . Cm 
of the complementary solution must be determined to satisfy 
the initial conditions or configuration. Heaviside, by aid of 
his conjugate Energy Theorem, has given an elegant and rela- 
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lively simple expression for these constants in terms of the initial 
currents in the inductances and the initial charges on the con¬ 
densers of the network.^ 

We shall base our further analysis of circuit theory on the solu¬ 
tion of problem (2), for the following reasons: 

(A) It is essentially a generalization of the Heaviside problem 
and its solution will furnish us a key to the correct understanding 
and interpretation of operational methods and lead to an auxiliary 
formula from which the rules of the operational calculus are 
directly deducible. 

(B) The solution of problem (2) carries with it the solution of 
problem (3) and also serves as a basis for the theory of alternating 
currents. 

(C) The solution of problem (2) leads directly to an extension of 
circuit theory to the case where the network contains variable ele¬ 
ments: i.e., circuit elements which vary with time and in which 
non-linear relations obtain. 

Problem (2) is therefore the fundamental problem of circuit 
theory and the formula which we shall now derive may be termed 
the fundamental formula of circuit theory. 

Consider a network in any branch of which, say branch 1, a 
unit e.m.f. is inserted at time t = 0, the network having been 
previously in equilibrium. By unit e.m.f. is meant an electro¬ 
motive force which has the value unity for all positive values of 
time (i>0). Let the resultant current in any branch, say branch 
rif be denoted by Ani{t). A^iit) will be termed the indicial admit¬ 
tance of branch n with respect to branch 1—or, more fully, the 
transfer indicial admittance. 

The indicial admittance, aside from its direct physical signifi¬ 
cance, plays a fundamental rble in the mathematical theory of 
electric circuits. In words, it may be defined as follows: The 
indicial admittance, AniCO? is equal to the ratio of the current 
in branch n, expressed as a time function, to the magnitude of 
the steady e.m.f. suddenly inserted at time i = 0 in branch 1. 
It is evidently a function which is zero for negative values of 
time and approaches either zero or a steady value (the d.c. admit- 

^ The reader who wishes to study this problem further is referred to the 
following papers: 

K. W. Wagner (E. N. T., Heft 11, Band 2, 1925). 

M. S. Vallarta (A. I. E. E., Apr,, 1926, “Heaviside^s Proof of his Expansion 
Theorem”). 
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tance) for all actual dissipative systems, as t a^pproaches infinity. 
It may be noted that, aside from its mathematical determination, 
which will.engage our attention later, it is an experimentally 
determinable function. 

We note, in passing, an important property of the indicial 
admittance Ajk{i)j which is deducible from the reciprocal theo¬ 
rem this is that A jk (0 = Akjit) . That is to say, the value of the 
transfer indicial admittance is unchanged by an interchange of 
the driving point and receiving point. It is therefore immaterial 
in the expression A^kit) whether the e.m.f, is inserted in branch 
j and the current measured in branch k, or vice-versa. In 
general, unless we are concerned with particular branches, the 
subscripts will be omitted and we shall simply write A{i), it 
being understood that any two branches or a single branch (for 
the case of ^equal subscripts) may be under consideration. 

From the linear character of the network, it is evident that if a 
steady e.m.f. E — E^is inserted at time i=r, the network being in 
equilibrium, the resultant current is 

Er A 

Generalizing still further, suppose that steady e.m.fs. Ei, 
jE 2 , . . . En are impressed in the same branch at the respective 
times 0, ri, t 2 . . . rn; the resultant current is evidently 

n 

E, A it) +Ei A it-Tx) + --+EUit -T„) = '^EiA it -r,-). (18) 

j 

To apply the foregoing to our problem we suppose that there is 
applied to the network, initially in a state of equilibrium, an 
e.m.f. E(;t) which has the following properties. 

1. It is identically zero for i{<0, 

2. It has the value E($) for 0<t<At 

3. It has the value E($)+AiE for At<t<2At 

4. It has the value eI%)+AiE+A 2 E for 2Ai<t<SAtr 


In other words it has the increment AjE at time i—jAtJ 
Evidently then the resultant current I (t) is 

EQA(t)'{-A.iEA(t —AQ-f-* * A'A^nE,A(t — nAt). 

^ Exceptions to this relation exist where the network contains sources of 
energy such as amplifiers. These need not engage our attention here. 
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Now evidently if the interval At is made shorter and shorter, 
then in the limit At-^dt and jAt = r and 

AiE=~E(r)dT. 

Passing to the limit in the usual manner this summation becomes 
a definite integral and we get , 

I(t) =E(,0)Ait)+£A{t-r)f^ E{r)dT. (19)'^ 

Finally by obvious transformations of the expression 
at the fundamental formula of circuit theory^ 

7(i)=| £Ait~r)E{r)dr, 

-sfm-^Wr)dr. 

For completeness we write down the following equivalents^ of 
(20) and (20a) 


7(0 =A{0)Eit)+ £A'{t-r)Eir)dr, 

(20b) 

= A{Q)Eit) + A'(t)N(« -r)dT, 

(20c) 

= Ei0)A (0 +£^'(f (^)^’'> 

(20d) 

= E{0)A(t)+j 'E'ir)Ait-r)dr. 

(20e) 


where the primes denote differentiation with respect to the 
argument. Thus A'(t)==d/dtA(t ). 

These equations are the fundamental formulas which mathe¬ 
matically relate the current to the type of applied electromotive 
force and the constants and connections of the system, and 

1 This important theorem, independently derived and published by the 
author, is actually the equivalent of a much older theorem in dynamics due 
to Duhamel. Substantially the same proof of the theorem as given in the 
text was iadependently communicated to the author by H. W. Nichols and 
by Stuart Ballantine. 

2 In case either A(t) or E{t) has discontinuities for positive values of 
we integrate through these discontinuities in the usual manner. When tiik 
is borne in mind, equations (20b) to (20e) are not restricted to continuou 
functions. 


we arrive 

( 20 ) 

(20a) 
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constitute the first part of the solution of our problem. The 
most important immediate deductions from these formulas are 
expressed in the following theorems. 

1. The indicial admittance of an electrical network completely 
determines, within a single integration, the behavior of the net¬ 
work to all types of applied electromotive forces. As a corollary, 
a knowledge of the indicial admittance is the sole information 
necessary to completely predict the performance and character¬ 
istics of the system, including the steady state. 

2. The applied e.m.f. and the indicial admittance are similarly 
and coequally related to the resultant current in the network. 
As a corollary the form of the current may be modified either by 
changing the constants and connections of the network or by 
modifying the form of the applied e.m.f. 

3. Since the applied e.m.f. may be discontinuous these formulas 
determine not only the building up of the current in response to 
an applied e.m.f. but also its subsidence to equilibrium when the 
e.m.f. is removed and the network left to itself. In brief, formu¬ 
las (20) reduce the whole problem to a determination of the 
indicial admittance of the network. In addition, as we shall 
see, they lead directly to an integral equation which determines 
this function. 

It is of interest to show the relation between formulas (20) 
and the usual steady state equations. To do this let the e.m.f., 
applied at time / = 0, be E sin (coi+0). Substitution in for¬ 
mula (20c) and rearrangement gives 


I{t)—A(0)E sin (cof-|-0) 

-\-E sin {(at-j-O) J cos corA'(T)dT 

— j? cos (a;Z+0) sin corA'(r)dr 


d 




where 

Now this can be resolved into two parts 

E sin (co^+0) I A(0)+ J' cos COT A'(r) dr 


^00 

I sir 
Jo 


E cos (cof+d) ^ I sin corA^(r)dr 


( 21 ) 


(22) 
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which is the final steady state, and 


-E sin (cot+6) 


00 

COS coTA^(T)dr 


noo 

+E cos (cot+&) J sin corA^(T)dr 


( 23 ) 


which is the transient distortion, which ultimately dies away for 
sufficiently large values of time. 

To correlate the foregoing expressions for the steady state with 
the usual formulas we observe that if the symbolic impedance of 
the network at frequency a)/27r be denoted by Z(ioi), and if we 
write 

^ = «(co)+i^(a>) 

then the steady state current is 

E [q:(co) . sin (aj^+0)+i3(co) . cos (co^+0)]. 

Comparison^with_(22)- gives_at once 

^00 

a(a))=A(0)+| COS cor A'(r) dr, (24) 

^00 

jS (co) = — J sin oirA \r)dr. (26) 

The Integral Equation for the Indicial Admittance.—So far 
we have tacitly assumed that the indicial admittance is known. 
As a matter of fact its determination constitutes the essential 
part of our problem. It is, in fact, the Heaviside problem, and its 
investigation, to which we now proceed, will lead us directly to 
the operational calculus. 

Heaviside^s method in investigating this problem was intuitive 
and “experimental.’^ We, however, shall establish a very gen¬ 
eral integral equation from which we shall directly deduce his 
methods and extensions thereof. 

Let us suppose that an e.m.f. where p is either a positive real 
quantity or complex with real part positive, is suddenly impressed 
on the network at time t==0. It follows from the foregoing theory 
that the resultant current 7(0 will be made up of two parts, (1) a 
forced exponential part which varies with time as and (2) a 
complementary part which we shall denote by y(t). The expo¬ 
nential or “forced’' component is simply e^^!Z(p), where Z(p) 
is functionally of the same form as the usual symbolic or complex 
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impedance Z(ico), It is gotten from the differential equations 
of the problem, as explained in a preceding chapter, by replacing 
by p^j cancelling out the common factor solving 

the resulting algebraic equations. The complementary or char¬ 
acteristic component, denoted by y(J;)j depends on the constants 
and connections of the network, and on the value of p. It does 
not, however, contain the factor and it dies away for suffi¬ 
ciently large values of tj in all actual dissipative systems. Thus 




(26) 


Now return to formula (20a) and replace E(t) by We get 


( A(T)e-^-dT 




which can be written as 
d 
dt 




A(T)e-P’-dT— A(T)e-’^dT 
Carrying out the indicated differentiation this becomes 


m 


^00 ^00 

= peP^ J A(r)e-^^dr-pe^^ J A(T)e~^^dT+A{t). (27) 


Equating the two expressions (26) and (27) for I (t) and dividing 
through by we get 

‘;^\+y(t)e~^^ = p r A(r)e'~^^dr—p f u4(r)e~’^’’dr+j4(^5)c~^^ (28) 

^{P) Jo Jt 


This equation is valid for all values of t. Consequently if we set 
03 and if the real part of p is positive, only the first term on 
the right and the left hand side of the equation remain, the rest 
vanishes, and we get 


pZ(p) 



(29) 


This is an mtegral equation^ valid for all positive real values of p, 
which completely determines the indicial admittance A{t). It is 


^ An integral equation is one in which the unknown function appears under 
the sign of integration. (29) is an integral equation of the Laplace type. If 
Z{p) is specified, A{t) is uniquely determined. Methods for solving the 
integral equations are considered in detail later, in connection with the 
exposition of the operational calculus. The phrase all positive values of p ’ ’ 
will be understood as meaning all values of p in the right half of the complex 
plane. 
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on this equation that we shall base our discussion of operational 
methods and from which we shall derive the rules of the oper¬ 
ational calculus. Equations (20) and (29) constitute a complete 
mathematical formulation of our problem, and from them the 
complete solution is obtainable without further recourse to the 
differential equations, or further consideration of boundary 
conditions. 



CHAPTER III 


THE HEAVISIDE PROBLEM AND THE OPERATIONAL 

EQUATION 

By virtue of the analysis and formulas of the preceding 
chapters, we are now in a position to take up a systematic 
exposition and development of the Heaviside operational calculus. 
This analysis may be reviewed and summarized as follows: 

The determination of the current in a network in response to 
an electromotive force E{t)j impressed on the network at reference 
time i = 0 has been reduced to the mathematical solution of two 
equations: first the integral equation 



and second, the definite integral 

( 20 ) 

In equation (29), Z{p) is the impedance function of the network 
and the characteristics, i.e., constants and connections, of the 
network enter into the problem only through Zip), and, therefore, 
Zip) completely specifies the network. If p is replaced by icj, 
where o^/2ir is the frequency, then Ziiof) is the complex steady 
‘state or symbolic impedance of alternating current theory. 
Thus, if Zip) is specified, equations (20) and (29) constitute a 
complete mathematical formulation of the problem, and their 
solution automatically takes care of the equilibrium boundary 
conditions. That is to say, contrary to the classical method of 
solution, the explicit determination of integration constants is 
completely dispensed with. 

It will be observed that in deducing these equations we have 
merely postulated (1) the linear and invariable character of the 
network and (2) the existence of an exponential solution of the 
type e^^/Zip) for positive values of p. Consequently, while we 
have so far discussed these formulas in terms of the determination 
of the current in a finite network, they are not limited in their 

21 
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application to this specfic problem. In this connection it may 
be well to call attention explicitly to the following points. 

The formulas and methods deduced above apply not only to 
finite networks, involving a finite system of linear equations, but 
to infinite networks and to transmission lines, involving infinite 
systems of equations, and partial differential equations: in fact 
to all electrical and dynamical systems in which the connections 
and constants are linear and invariable. 

Secondly, the variable determined by formulas (20) and (29) 
need not, of course, be the current. It may equally well be the 
charge, potential drop, or any of the variables with which we 
may happen to be concerned. This fact may be explicitly recog¬ 
nized by writing the formulas as: 


pHivri 

<i)=‘Jt£Kt-r)E{r)dr. 


(30) 

(31) 


Here E(it) is the applied e.m.f., x{t) is the variable which we desire 
to determine (charge, current, potential drop, etc.), and 

x^E/H{v) (32) 

is the operational equation. E (p) therefore corresponds to and is 
determined in precisely the same way as the impedance Z{'p), but 
it may not have the physical significance or the dimensions of an 
impedance. Similarly, in character and function, h{t) corre¬ 
sponds to the indicial admittance, though it may not have the 
same physical significance. It is a generalization of the indicial 
admittance and may be appropriately termed the Heaviside 
Function, Similarly Hi'p) may be termed the generalized 
impedance function. 

The Operational Equation.—The physical problem which 
Heaviside attacked and which led to his operational calculus was 
the determination of the response of a network or electrical sys¬ 
tem to a '^unit e.m.f.(zero before, unity after time f=0) with, 
of course, the understanding that the system is in equilibrium 
when the electromotive force is applied. His problem is therefore, 
essentially that of the determination of the indicial admittance. 
In our exposition and critique of Heaviside's method of dealing 
with this problem we shall accompany an account of his own 
method of solution with a parallel solution from the corresponding 
integral equation of the problem. 



HBAVIBIDE PROBLEM AND OPERATIONAL EQUATION 23 

Heaviside^s first step in attacking this problem was to start 
with the differential equations, and replace the differential opera¬ 
tor d/dt by the symbol p, >nd the operation Sdi hy h/p, thus 
reducing the equations to^an algebraic form. He then wrote the 
impressed e.m.f. as 1 (unity), thus limiting the validity of the 
equations to values of ^>0. The formal solution of the algebraic 
equations is straightforward and will be written as 

= (33) 

whe^e h is the ^^generalized indicial admittance,*^ or Heaviside 
function (denoting current, charge, potential or any variable 
with which we are concerned) and i?(p) is the corresponding 
generalized impedance. Thus, if we are concerned with the 
current in any part of the network, we write 

A = l/Z(p). (34) 

The more general notation is desirable, however, as indicating 
the wider applicability of the equation. 

The equations 

A = l/i?(p) 

A^UZiv) 

are the Heaviside operational equations. They are, as yet, purely 
symbolic and we have still the problem of determining their 
explicit meaning and in particular the significance of the operator 
V* 

Comparison of the Heaviside operational equations with the 
integral equations (29) and (30) of the preceding chapter leads to 
the following, fundamental theorem. 

The Heaviside operational equations 

A^l/Z{p) 

h^l/H{p) 

are merely the symbolic or short-hand equivalents of the correspond¬ 
ing integral equations 

The integral equations^ therefore, supply us with the meaning and 
significance of the operational equations, and from them the rules 
of the operational calculus are deducible. 
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By virtue of this theorem, we have the advantage, at the 
outset, of a key to the meaning of Heaviside's operational 
equations, and a means of checking and deducing his rules of 
solution. This will serve us as a guide throughout our further 
study. 

Eeturning now to Heaviside's own point of view and method of 
attack, his reasoning may be described somewhat as follows:— 
The operational equation 

is the full equivalent of the differential equations of the problem 
and must therefore contain the information necessary to the 
solution provided we can determine the significance of the sym¬ 
bolic operator p. The only way of doing this, when starting 
with the operational equation, is one of induction: that is, we 
must compare the operational equation with known solutions of 
specific problems and thus attempt to infer by induction general 
rules for interpreting the operational equation and converting 
it into the required explicit solution. 

The Poiuer Series Solution 

Let us start with the simplest possible problem; the current in 
response to a “unit e.m.f." in a circuit consisting of an inductance 
L in series with a resistance R. 

The differential equation of the problem is 

Lj^A+RA = l, t>0, 


where A is the indicial admittance. Consequently replacing 
d/dt by p, the operational equation is 

pL+R 

The explicit solution is easily derived: it is 
A = 


where a=R/L. Note that this makes the current initially zero, 
so that the equilibrium boundary condition at « = 0 is satisfied. 

Now suppose that we expand the operational equation in 
inverse powers of p: we get, formally. 


i ^ 1 o: 1 _1 O' 

pLl+a/p Rpl+cx/p~R"p 




by the Binomial Theorem. 
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Now expand the explicit solution as a power series in it is 


[ at 

{aty 



2 ! 

^ 3 ! 


Comparing the two expansions we see at once that the operational 
expansion is converted into the explicit solution by assigning 
to the symbol 1/p” the value It was from this kind of 

inductive inference that Heaviside arrived at his power series 
solution. 

Now there are several important features in the foregoing which 
require comment. In the first place the operational equation is 
converted into the explicit solution only by a particular kind of 
expansion, namely an expansion in inverse powers of the operator 
p. For example, if in the operational equation 


A 


1 afp 
R l+ce/p 


we replace 1/p by iJ/l! we get 


1 at 

Rl+at 


which is incorrect. Furthermore, if we expand in ascending 
instead of descending powers of p, namely 


A 


-^l^l-(p/a) + (p/ay 


no correlation with the explicit solution is possible and no signifi¬ 
cance can be attached to the expansion. We thus infer the gen¬ 
eral principle, and we shall find this inference to be correct, 
that the operational equation is convertible into the explicit 
solution only by the proper choice of expansion of the impedance 
function, or rather its reciprocal. 

In the second place we notice that in writing down the opera¬ 
tional equation and then converting it into the explicit solution 
no consideration has been given to the question of boundary 
conditions. This is one of the great advantages of the opera¬ 
tional method: the boundary conditions, provided they are those 
of equilibrium^ are autoriiatically taken care of. This will be 
illustrated in-the next example: 

Let a ^‘unit e.m.f.’’ be impressed on a circuit consisting of 
resistance 12, inductance L, and capacity C : required the resultant 
charge on the condenser. 
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The differential equation for the charge Q is 

<> 0 . 

Consequently the operational formula is 

Q = _?:_ 

^ Lp^+Rp+l/C 

“““ '"ro- 

This can be expanded by the Binomial Theorem as 



Performing the indicated operations and collecting in inverse 
powers of p, the first few terms of the expansion are:— 

1 ^ j_1 

Lp^\ V p® j 

where . ci = a 

C2 — h—a^ 

C3 = 2a6~-a® 

C4 = 6"-3a2&+a'' 

C5 = 3a6^“*4a^6+a® 

C 6 - 62+5^4 6-a® 


We infer therefore that in accordance with the rule of replacing 
1/p” by f^/n\ the solution is:— 



5^2 j^C 

^j-Cl3-,-C2j-,+C3g-j + C4g-,- ■ ■ 


Owing to the complicated character of the coefficients in the 
expansion, the series cannot be recognized and summed by inspec¬ 
tion. If, however, we put J? = 0 then a=0, and the series becomes 


C 


kUmi-iUm) 




whence 


Q = C'{l-cos (f/v'hC')}. 


We have still to verify this solution by comparison with the 
explicit solution of the differential equation. This is of the form 
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where kx and k 2 are constants which must be chosen to satisfy the 
boundary conditions and Xi, X 2 are the roots of the equation 

LX2+SX+1/C = 0. 

Now since we have two arbitrary constants we satisfy the 
equilibrium conditions by making Q and dQ/dt zero at ^==0. 
whence 

= 0 , 

XiAji“hX2 /i‘2 “ 0, 

and 

kx = X2 C/ (Xi ” X2), 

= Xi ( 7 / (X2—Xi). 

We have also 



Writing down the power series expansion of 

Q = C +/cie^^^+/c2e^®^, 

then 

+ (feiXl^ + fc2X2^)^j+ • • * ■ 

Introducing the values of /ci, /c 2 , Xi; X 2 given above and comparing 
with the power series derived from the operational solution we see 
that they are identical term by term. 

This example illustrates two facts. First the power series 
expansions may be complicated, laborious to derive and of such 
form that they cannot be recognized and summed by inspection. 
In fact in arbitrary networks of a large number of meshes or 
degrees of freedom the evaluation of the coefficients of the power 
series expansion is extremely laborious. 

On the other hand, in such cases, the solution by the classical 
method presents difficulties far more formidable—in fact insuper¬ 
able difficulties from a practical standpoint. First there is the 
location of the roots of the function /f(X)j which in arbitrary 
networks is a practical impossibility without a prohibitive amount 
of labor. Secondly there is the determination of the integration 
constants to satisfy the imposed boundary conditions: a process, 
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which, while theoretically straightforward, is actually in practice 
extremely laborious and complicated. We note these points in 
passing; a more complete estimate of the value of the power 
series solution will be made later. 

To summarize the preceding: Heaviside, generalizing from 
specific examples otherwise solvable, arrived at the following 
rule:— 

Expand the right hand side of the operational equation 
h^l/H{p) 

in inverse powers of p: thus 

h'^ao+ai/p+a2/p^-\ -b On/p^ H- 


and then replace — by t'^/n !. The operational equation is thereby 
converted into the explicit power series solution :— 


h'=aQ-\-ait/l l-\-a2t^/2\’-\- * * -\-ant^/nl’^ • * • ( 35 ) 

As stated above, this rule was arrived at by pure induction and 
generalization from the known solution of specific problems. It 
cannot, therefore, theoretically be regarded as satisfactorily 
estabhshed. The rule can, however, be directly deduced from 
the integral equation 


To its derivation from this equation we shall now proceed. 

First suppose we assume that h(t) admits of the power series 
expansion 

hij-\-hi tjX ^^2 !“f-. 

Substitute this assumed expansion in the integral, and integrate 
term by term. The right hand side of the integral equation 
becomes formally 

ho/p+Ai/p2+/i2/2)3+. 

by virtue of the formula 


Now expand the left hand side of the integral equation asymp¬ 
totically in inverse process of p: it becomes 

ao/p+ai/p^+a2/p^+ • •; • • 


where 


(io+cii/p+a2/p^H - 
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is the asymptotic expansion of Comparing the two 

expansions and making a term by term identification, we see that 
hn = and 

^/l !“f'• • * * 

which agrees with the Heaviside formula. 

This procedure, however, while giving the correct result has 
serious defects from a mathematical point of view. For example, 
the asymptotic expansion of 1/H{p) has usually only a limited 
region of convergence, and it is only in this region that term by 
term integration is legitimate. Furthermore we have asmmed 
the possibility of expanding h(t) in a power series: an assumption 
to which there are serious theoretical objections, and which, 
furthermore, is not always justified. A more satisfactory deriva¬ 
tion, and one which establishes the condition for the existence of a 
power series expansion, proceeds as follows:— 

Let 1/H{p) be a function which admits of the formal asymp¬ 
totic expansion 

CO 

0 i 

and let it include no component which is asymptotically repre¬ 
sentable by a series all of whose terms are zero, that is a function 
(l>(p) such that the limit, as oo, of <j5)(p) is zero for every value 
of n. Such a function is 6“^. With this restriction understood, 
start with the integral equation, and integrate by parts: we get 

where denotes Now let p approach infinity: 

in the limit the integral vanishes and by virtue of the asymptotic 
expansion 

• 00 

10 | 

l/Hip) approaches the limit Oo- Consequently 

A(0) = ao. 

Now integrate again by parts: we get 

p(l/i?(p) -oo) = A(«(0)+ 
use Lib B’lore ° 

621.31921 N26 


30 


MLECTRIC CIRCUIT THE6RY 


Again let p approach infinity: in the limit the left hand side oJ 
the equation becomes oi and we have 

Proceeding by successive partial integrations we thus establisl: 
the general relation 

¥''K0)=a.. 

But by Taylor’s theorem, the power series expansion of hit) if 
simply 

Kt) =A(0) +Af«(0)i/l!+/i^«(0)«V2!+. 

whence, assuming the convergence of this expansion, we get 

00 

A(i) = ao+Oif/lH" ctn t^/n ! (35) 

0 

which establishes the power series solution. It should be 
carefully noted, however, that it does not establish the con¬ 
vergence of the power series solution. As a matter of fact, 
however, I know of no physical problem in which H{p) satisfies 
the conditions for an asymptotic expansion, where the power 
series solution is not convergent. On the other hand many 
physical problems exist, including those relating to transmission 
lines, where a power series solution is not derivable and does not 
exist. 

The process of expanding the operational equation in such a 
form as to permit of its being converted into the explicit solution 
is what Heaviside calls “algebrizing” the equation. In the 
case of the power series solution the process of algebrizing 
consists in expanding the reciprocal of the impedance function 
in an asymptotic series,^ thus 

1 /H (p) ~oo -f ai/p -faa/p^ 4-. 

Regarded as an expansion in the variable p, instead of as a 
purely symboHc expansion, this series has usually only a limited 
region of convergence. This fact need not bother us, however, as 
the series we are really concerned with is 

aa-\-axt/l\-\-aiP/'2i\-\- • • • • 

* The symbol frequently employed in the following pages, is to be 
understood as signifying ‘‘is asymptotically equal to.” In other words the 
two sides of the equation are not, in general, identities but are asymptotically 
equivalent or equal. 
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It is interesting to note in passing that the latter series is what 
Borel, the French mathematician, calls the associated function 
of the former, and is extensively employed by him in his 
researches on the summabilifcy of divergent series. 

The process of ^^algebrizing,^^ as in the examples discussed 
above, may often be effected by a straightforward binomial 
expansion. In other cases the form of the generalized impedance 
function H(p) will indicate by inspection the appropriate pro¬ 
cedure. A general process, applicable in all cases where a power 
series exists, is as follows. Write 

1/H{p) = 1/h(^=G(,x). ( 36 ) 

Now expand G(x) as a Taylor^s series: thus formally 


where 


(?W(0) = 


dx' 




( 37 ) 


:c=0. 


Qin) (Q) 

Denote —by an, replace by 1/p”, and we have 


G{x) = 1/H{p) = ao+ai/p+a 2 /p^A - 

This process of “algebrizing^^ is formally straightforward and 
always possible. As implied above, however, in many problems 
much shorter modes of expansion suggest themselves from the 
form of the function H(p).. 

We note here, in passing,*^ that the necessary and sufficient 
conditions for the existence of a power series solution is the 
possibility of the formal expansion of G(x) as a power series in x. 

At this point a brief critical estimate of the scope and value of 
the power series solution may be in order. As stated above, in a 
certain important class of problems relating to transmission 
lines, a power series does not exist, though a closely related 
series in fractional powers of t may often be derived. Con¬ 
sequently the power series solution is of restricted applicability. 
Where, however, a power series does exist, in directness and 
simplicity of derivation it is superior to any other form of solu¬ 
tion. Its chief defect, and a very serious defect indeed, is that 
except where the power series can be recognized and summed, 
it is usually practically useless for computation and interpreta¬ 
tion except for relatively small values of the time This 
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disadvantage is inherent and attaches to all power series solu¬ 
tions. For this reason I think Heaviside overestimated the 
value of power series as practical or working solutions, and that 
some of his strictures against orthodox mathematicians and 
their solutions may be justly urged against the power series 
solution. He was quite right in insisting that a solution must 
be capable of either interpretation or computation and quite 
right in ridiculing those formal solutions which actually conceal 
rather than reveal the significance of the original differential 
equations of the problem. On the other hand, the following 
remark of his indicates to me that Heaviside had a quite exagger¬ 
ated idea of the value and fundamental character of power 
series in general: regret that the result should be so compli¬ 

cated. But the only alternatives are other equivalent infinite 
series, or else a definite integral which is of no use until it is 
evaluated, when the result must be the series (135), or an equiv¬ 
alent one.” As a matter of fact the properties of most of the 
important functions of mathematical physics have been investi¬ 
gated and their values computed by methods other than series 
expansions. I may add that in technical work the power series 
solution has proved to be of restricted utility, while definite 
integrals, which Heaviside^ particularly despised, have proved 
quite useful. 

The Expansion Theorem Solution^ 

We pass now to the consideration of another extremely 
important form of solution. Heaviside gives this solution with¬ 
out proof: we shall therefore merely state the solution and then 
derive it from the integral equation. 

Given the operational equation 

h^l/H(p) 

which has the significance discussed above: i.e., the response of 
the network to a '^unit e.m.f.” The explicit solution may be 
written as 

‘-gro)+2w 

1 

^ Vide a remark of his to the effect that some mathematicians took refuge 
in a definite integral and called that a solution. 

* This terminology is due to Heaviside. A more appropriate and physi¬ 
cally significant expression would be “The Solution in terms of normal or 
characteristic vibrations/' 
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where pi, 
and 


Pn are the n roots of the equation 

/f(p) = 0 


H'ip.) 


±H{p) 

_dp 


(39) 


As remarked above, this solution, referred to by him as The 
Expansion Theorem^ was stated by Heaviside without proof; 
how he arrived at it will probably always remain a matter of 
conjecture.^ Its derivation from the integral equation is, 
however, a relatively simple matter, though in special cases 
troublesome questions arise. 

The derivation of the expansion solution from the integral 
equation 


follows immediately from the partial fraction expansion 

n 

_1_^_1_ I _1_ 

pH ip) pHlo) ^{p ~Pi)pi H'ipi) 

3 = 1 


(40) 


where pi, pi . . . are the roots of the equation 2?(p)=0, 
and 


H'iPj) 


[dp 


H(j>) 


JP«Pj 


(41) 


Partial fraction expansions of this type are fully discussed in 
treatises on algebra and the calculus and the conditions for their 
existence established. Before discussing the restrictions imposed 
on H(p) by this expansion, we shall first, assuming its existence, 
derive the expansion theorem solution. 

By virtue of (40) the integral equation is 

n 

WW) 2 'W) ~ X 

1 


The expansion on the left hand side suggests a corresponding 
expansion on the right hand side: that is, we suppose that 

7i(^) ==/io(i)+^i(0+^2(0 +* * *+^n(0 (43) 

^ Since this was written it has been independently shown by K. W. 
Wagner and M. S. Vallarta (see footnote to p. 14) that a proof based on his 
^‘conjugate theorem” can be reconstructed from scattered papers of 
Heaviside. 
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and specify that these component functions shall satisfy the 
equations 

i,- (44) 




It follows at once from (43) and direct addition of equations (44) . 
and (45) that (42) is satisfied and hence is solved provided 
Ao, * . hn can be evaluated from (44) and (45). 

Now since 

(46) 


provided the real part of X is not positive (a condition satisfied 
in all network problems), we see at once that equations (42) 
and (43) are satisfied by taking 

7 . 7 1 /A}-T\ 






y = l,2,. .n. 


Consequently from (43) and (47) it follows that 




PiH'iPi) 


which establishes the expansion theorem solution. 

As implied above, the partial fraction expansion (40), on 
which the expansion theorem solution depends, imposes certain 
restrictions on the impedance function H(p). Among these are 
that H(p) must have no zero root, no repeated roots, and 1/H(p) 
must be a proper fraction. In all finite networks these conditions 
are satisfied, or by a slight modification, the operational equation 
can be reduced to the required form. The case of repeated 
roots, which may occur where the network involves a unilateral 
source of energy such as an amplifier, can be dealt with by 
assuming unequal roots and then letting the roots approach 
equality as a limit. Without entering upon these questions in 
detail, however, we can very simply and directly establish the 
proposition that the expansion theorem gives the solution 
whenever a solution in terms of normal or characteristic vibrations 
exists. The proof of this proposition proceeds as follows. 
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It is known from the elementary theory of linear differential 
quations that the general solution of the set of differential 
quations, of which the operational equation is h = l/H{p)^ is of 
ae form 

n 

h{t) = Co+'^Cie^^‘ 

I 

rhere 'Pi is the jth root of if(p)=0, and Co, . . . Gn are 
onstants of integration which must be so chosen as to satisfy 
he system of differential equations and the imposed boundary 
onditions. The summation is extended over all the roots of 
I (p) which is supposed not to have a zero root or repeated roots. 

Now substitute this known form of solution in the integral 
quation of the problem and carry out the integration term by 
erm. We get 

letting p = 0, we have at once 

Co = l/f/(0). (50) 

To determine Cj let p = where qi^ 2 . small quantity ulti- 
nately to be set equal to zero, and write the equation as 

= (51) 

f now v-'pi'^-q and q approaches zero, this becomes in the 


imit 



)r 


(52) 

c — ^ ) 

(53) 

vhence 


(54) 

vhich is 

the Expansion Theorem Solution. 



We shall not attempt to discuss here cases where the expansion 
solution breaks down though such cases exist. In every such 
jase, however, the breakdown is due to the failure of the imped- 
mce function H{p) to satisfy the conditions necessary for the 
partial fraction expansion- (40), and correlatively the non- 
3 xistence of a solution in normal vibrations. Furthermore, it is 
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usually possible by simple modifications to deduce a modified 
expansion solution. It may be added here, that while the 
proof given above is also limited implicitly to finite networks, 
the expansion solution is valid in most transmission line problems. 

Let us now illustrate how the expansion solution works by 
applying it to a few simple examples. Take first the case 
considered in the preceding chapter in connection with the 
power series solution. Required the charge Q on a condenser 
C in series with an inductance L and resistance R in rseponse 
to a ^^unit e.m.f.^' The operational equation is 

^ Lp^+Rp+l/C 

“ ^ Lp^+2apW 

where a = i2/2L and = 1/LC. 

The roots of the equation H{p)—0 are the roots of the equation 

p^+2ap'i-(a^=0 

whence 

Pi= "~a+\/a2-“C0^= — 

Also H'(p) = 2L{p+a), so that 

ir'(px)=2/3L 

H'(p2)==-2IIL 

and 

l/H(0)==l/Lco2 = C. 

Inserting these expressions in the Expansion Theorem Solution 
(38), we get 

It is now easy to verify the fact that this solution satisfies the 
differential equations and the boundary condition 

Q = 0 and dQ/dt = 0 at time ^ = 0. 

If w > O', /3 is a pure imaginary 

/? = ico\/l — (oj/co) ^ — ito' 
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In connection with this problem we note two advantages of the 
expansion solution, as compared with the power series solution: (1) 
it is much simpler to derive from the operational equation, and (2) 
its numerical computation is enormously easier, A table of expo¬ 
nential and trigometric functions enables us to evaluate Q for any 
value of t almost at once whereas in the case of the power series 
solution the labor of computation for large values of i is very great. 
A third and very important advantage of the expansion solution 
in this particular problem is that without detailed computation 
we can deduce by mere inspection the general character of the 
function and the effect of the circuit parameters on its form: 
an advantage which never attaches to the power series solution. 

This last property of the particular solution above is extremely 
important* The ideal form of solution, particularly in technical 
problems, is one which permits us to infer the general character 
and properties of the function and the effect of the circuit con¬ 
stants on its form, without detailed computation. A solution 
which possesses these properties, even if its exact computation is 
not possible without prohibitive labor, is far superior to a solution 
which, while completely computable, tells us nothing without 
detailed computation. It is for this reason that some of the 
derived forms of solution, discussed later, are of such importance. 
In fact a solution which requires detailed computation before it 
yields the information implied in it is merely equivalent to an 
experimentally determined solution. 

Unfortunately the advantages attaching to the expansion 
solution of the specific problem just discussed, do not, in general, 
characterize the expansion solution. The following disadvan¬ 
tages should be noted. First, the location of the roots of the 
impedance function J?(p) is practically impossible in the case of 
arbitrary networks of more than a few degrees of freedom. In 
the second place, when the number of degrees of freedom is large 
it is not only impossible to deduce the significance of the solution 
by inspection, but the computation becomes extremely laborious. 
In such cases, the practical value of the expansion solution 
depends, just as in the power series solution, on the possibility of 
recognizing and summing the expansion. This will be clear 
in the case of transmission lines, where the roots of Hip) are 
infinite in number and the direct computation of the expansion 
solution (except in the case of the non-inductive cable) is quite 
impossible. 



CHAPTER IV 


SOME GENERAL FORMULAS AND THEOREMS FOR THE 
SOLUTION OF OPERATIONAL EQUATIONS 

We have seen that the operational equation 
is the symbolic or short-hand equivalent of the integral equation 



and from the latter we have deduced two very important forms of 
the Heaviside solution. In recognizing the equivalence of these 
two equations we have a very great advantage and are able, in 
fact, to base the operational calculus on deductive instead of 
inductive reasoning. In this chapter we shall employ this equiv¬ 
alence to establish certain general formulas and theorems for 
the solution of operational equations. That is to say, we shall 
make use of the principles that (1) any method applicable to 
the solution of the integral equation supplies us with a corre¬ 
sponding method for the solution of the operational equation, 
and (2) a solution of any specific integral equation gives at once 
the solution of the corresponding operational equation. We 
turn therefore to a brief discussion of the appropriate methods 
for solving the integral equation. 

It has been noted above that any solution of the integral equa¬ 
tion supplies a solution of the corresponding operational equation. 
This principle enables us to take advantage of the fact that a very 
large number of infinite integrals of the type 



have been evaluated. The evaluation of every infinite integral of 
this type supplies us, therefore, with the solution of an operational 
equation. 

Of course, not all the operational equations so solvable have 
physical significance. Many however, do. Below is a list of 
infinite integrals with their known solutions, accompanied by 
the corresponding operational equation and its explicit solution. 

38 
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All of these solutions are directly applicable to important tech¬ 
nical problems. It may be remarked in passing that the infinite 
integrals have for the most part been evaluated by advanced 
mathematical methods which need not concern us here. 


Table of Infinite Integrals, the Corresponding Operational Equa¬ 
tions, and Their Explicit Solutions 


(a) 

(b) 

(c) 

(d) 

(e) 

(f) 

(g) 

(h) 


/: 




P+\ 


h = -^ = e-K 

p+\ 


00 

0^ 


-pt 






pn ^ 


/: 


1 , 1 
— dt=- 


6-^-—7=-.dJ = — 
■VTri VP 

h = '\/p^lf\/TvL 


£ 


A= 


-pi . 


(2iy 


dt 


1.3.5* •* (2?2 — 1) v^Tj.^ p^^-\/p 
Vp_ (2t)" 1 


p” 1.3.5-• •(2n-l)V7rt 
I 6“*" e~^‘ dt — 7— ; -rT.- 

Jo n! (p-fX)"+i 

/j=_?._=?!!.. g-w 

(p+x)”+' nr • 


h==pe-~^’^^P= yj 

poo 

e- 
Jo 


-4t = e-2V\p^ 


p-\!t 

■p‘i — dt= 


■^tVi 

g——2 "v/ Xp 


V^t 


Vv 

^ g-XA 

h = Vp 6”^ 'Vxp = — . 

V xi 


I 




sin Xi dt 


p2+X^ 


p2+X2 

sin X^. 
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(i) 


^00 

I cos \t — 

Jo 




= cos Xi. 


J '»00 


■Pt Q-pi gQg _ 




(P+m)HX2 


cos Xf. 


<« /; 


(p+m)“+v 
e-p‘ e~>‘* sin M dt = 


(?>+m)^+X^ 


( 1 ) 


/; 


pX 

(P+m)"+'x" 
e~^‘ Jo0^t)dt- 


=6“"* sin X<. 


Vp^+X^ 


(m) re-p{JB{Vi^^)dt = 


Vp*+i* 


/i = _..^__g-xVp“+i = 0 for i<X 

Vp^+1 

=/o(V«“-X2) for «>X. 


(n) 


(P) 


/: 


e~'^^Jn0^t)dt 


=F,(?)'- 


p^+x^, 


r\ X 


=/n(X0. 


r 


7i = 


e“®*e“^‘/o(X<)(if= 
1 


Vl+2X/p 


a/ p2+2Xp 

==e-^^Jo(XO. 


In formulas (1), (m), (n), Jn{x) denotes the Bessel function of 
order n and argument x. In formula (p), Io(x) denotes the Bessel 
function Jo{ix) where z = V— 1- 

This list might be greatly extended. As it is, we are in posses¬ 
sion of a set of solutions of operational equations which occur in 
important technical problems and which will be employed later. 

The foregoing emphasizes the practical and theoretical impor¬ 
tance of recognizing the equivalence of the integral and operational 
equations. With this equivalence in mind, the solution of an 
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operational equation is often reduced to a mere reference to a table 
of infinite integrals. Heaviside did not recognize this equivalence. 
As a consequence many of his solutions of transmission line 
problems are extremely laborious and involved and in the end 
unsatisfactory because expressed in involved power series. 

Not all the infinite integrals corresponding to the operational 
equations of physical problems have been evaluated or can be 
recognized without transformation. This statement corresponds 
exactly with the fact that a table of integrals is not always suffi¬ 
cient but must be supplemented by general methods of integra¬ 
tion. We turn, therefore, to stating and discussing some general 
theorems applicable to the solution of operational equations. 

In the derivation of the operational theorems, which constitute 
the general rules of the operational calculus, the following proposi¬ 
tion, due to Borel and known as BorePs theorem, will be fre¬ 
quently employed.^ 

If the functions f(t)^ /i(0; €indf^{t) are defined by the integral 
equations 

^ oo 

^’(p) = Jo 

/'\(p) = J^ fi{t)erp‘dt 

and if the functions F, Fi and F% satisfy the relation 
F(.p)=Fiip).Fi{p) 

then 

fit) = jjiir)f2it-T)dr 

\ir)fiit-r)dT. 

The operational theorems will now be stated and briefly proved 
from the integral equation identity. 

I For a proof of this important theorem the reader is referred to Borel, 
“Lecons sur les Series Divergentes” (1901), p. 104; to Bromwich, “Theory 
of Infinite Series,” pp. 280-281; or to Ford, “Studies on Divergent Series and 
Summability,” pp. 93-94 (being Vol. II of the Michigan University Science 
Series, published by Macmillan). The proof depends on Jacobi’s trans¬ 
formation of a double integral: see Edward’s “Integral Calculus, 1922, 
Vol. II, pp. 14r-15. 
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Theorem I 

If in the operational equation 

h = l/H{v) 

the generalized impedance function H{p) can be expanded in a 
sum of terms, thus 

Hip) Hi(p)^H,ip)^ ^H„ip) 

and if the auxiliary operational equations 

^^°mip) 

can be solved^ then 

h = hi-\-h2'{‘ • • * ~hhn» 

This theorem is too obvious to require detailed proof: in fact 
it is self-evident. The power series and expansion theorem 
solutions are examples of its application. In general, however, 
the appropriate form of expansion of l/H(p) will depend on the 
particular problem in hand. The theorem, as it stands is a 
formal statement of the fact that solutions can often be obtained 
by an appropriate expansion whereas the equation cannot be 
solved as it stands. 

Theorem II 

If h=li(t) and g=g(t) are defined by the operational equations 
h^l/H{p) 
g=^l/pH{p) 

then 

git)^j\ir)dr. 

To prove this theorem we start with the integral equations 



The second of these is in form for an immediate application of 
BoreTs theorem since 

1 1 
p^H(p) p pH{p) 
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The functions fi and /a of Borel’s theorem then satisfy the 
equations 


1 

V 


1 

VH{f) 






It follows at once that 


/xW = l 

whence by Borel’s theorem 


Thisorem III 


'=X 


h(r)dr. 


If h=h(t) and g=g{i) are defined by the operational equations 

h = l/H(p) 
g = p/H(p) 

then 

provided /i(0) =0. 

The integral equations of the problem are 

Integrating the first of these by parts we have, 

-^-r = h(0)+X 

pH{p) P ^ ' pjo 
where h'{t)=^d/dt h(t). 

If A(0)==0, we have at once 

I h'{t)e ^^dt 

I-iiP) Jo 

Comparison with the integral equation for g{t) shows at once 
that g(i)~h'it), since the integral equation determines the 
function uniquely. 

Theorems II and III establish the characteristic Heaviside 
operations of replacing 1/p byj' dt and p by d/di. 
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Theorem IV 

If in the operational equation 

h = l/H{v) 

the generalized impedance function can be factored in the form 
H{p)=H,{p)-H,{p) 

and if the auxiliary operational equations 

h^=^l/H^{p) 

hi=l/E2{p) 

define the auxiliary variables hi and h^, then 


hit) 


-it. 
£ 


A 

dt, 


hi{r)h^{t—r)dr 


This theorem is immediately deducible from BoreFs theorem and 
theorems II and III, as follows. * 

The integral equations are 

pH ip) ~ ‘PpH lip)' pH 2 ip)~ hit)e ^ dt 


lip) Jo 

1_ _ r 

^2(P) Jo 


it)e-p‘dt 




Now define an auxiliary function g(t) by the operational 
equation 

1 


Then 


vim 

pHi(^'pH^~ I ^ 
and by BoreFs theorem 

^(0 = hi{r)h2{t-r)dr 

From this equation it follows that ^(0) ~0, and hence comparing 
the operational equations for h and g, we have by aid of Theorem 
III 

A ( 0=|(/(0 
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and hence 

h{t)=~ {\^{j)hi{t-r)dr 
dtJo 

= 5 rh2{r)h,{t-T)dr. 

dtJo 

This theorem is extremely important, although not stated or 
employed by Heaviside himself. We shall make use of it in estab¬ 
lishing two important general theorems and shall have frequent 
occasion to employ it in specific problems occurring in connection 
with the subsequent discussion of transmission theory, 

Theokem V 

If h = h{t) and g — g{t) are defined hy the operational equations 

, 1 

Mv) 

1 

^ /i(p+X) 

where X is a positive real.parameter, then 

To prove this theorem we start with the integral equations 

pirc^xrl 

In the first of these equations replace the symbol p by (?+X: 
we get 

and then to preserve our original notation replace the symbol q 
by p, whence 

(a) 

The integral equation in g{t) can be written as 

(i+^)(p+x)»(f+xrX ® 
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Comparing equations (a) and (b) it follows at once from 
theorems I and II that 

gf(0 = ^l+xj' d^h{t)e-^K 

From the foregoing, the following auxiliary theorem is immedi¬ 
ately deducible. 

Theorem Va 

If and g = g(i) are defined by the operational equations 

. JL 

H{v) 

_ V 
^ (p+X)H(p+X) 

then 

g{t) =hit)e-^‘. 

The proof of this theorem will be left as an exercise to the reader. 
Theorem VI 

If h — hit) and g=g{t) are defined by the operational equations 

h = l/H(p) 
g=l/H(Kp) 

where X is a positive real parameter, then 

g{t) = hit/\). 

We start with the integral equations 

and in the first of these equations we replace p by 'Kq and t by 
r/X, whence it becomes 

Now replacing the symbols q and r by p and t respectively, we 
have 

29F(Xp)“/p 

whence by comparison with the integral equation in ^(0 it 
follows at once that 


?W = A(^/X). 
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This theorem is often useful in making a convenient change in 
the time scale and eliminating superfluous constants. 


Theorem VII 


If h — h(t) and g = g(f) are defined by the operational equations 


, 1 
H{p) 


where X is a positive real quantity, then 
g{t)=Q for t<\ 


— h{t—\) for t>\. 

This is a very important theorem in connection with trans¬ 
mission line problems where retardation, due to finite velocity 
of propagation, occurs. Its proof proceeds as follows: 

If the auxiliary function h=^h{t) is defined by the operational 
equation 

A; = 6“'^^ 


then by Theorem IV, 

Kr)Ht-r)dT. 


(a) 


Now, corresponding to the operational equation /c = e 
have the integral equation 





we 


The solution of this integral equation, which is easily verified 
by direct substitution in the infinite integral, is 


k(t)=^0 iort<\ 
= 1 for t>\. 


Hence equation (a) becomes 

git)-0 ioTt<\ 
d 

^di \ for^>X 

=A(i—X) fori>X. 

Theorem IV, employed in the preceding proof, as stated above, 
is extremely important and we shall have frequent occasion to 
employ it in specific problems. We shall now apply it to deduce 
an important theorem which extends the operational calculus 
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to arbitrary impressed forces, whereas heretofore the operational 
equation h = l/H(p) applied only to the case of a '^unit e.mi/’ 
impressed on the system. 

It will be recalled from a previous chapter that if x{t) denotes 
the response of a network to an arbitrary force /(O? impressed 
at time i = 0, and if h{t) denotes the corresponding response to a 
^^unit e.m.f.,’^ then 


h(r)f(t-T)dT 

(31) 

1 

vHipri 

(30) 


Now /(O may be of such form that the infinite integral 





can be evaluated and has the value F{p)/p: thus 





(55) 


This is possible, of course, for many important types of applied 
forces, including the sinusoidal. 

It follows at once from Theorem IV that x(t) satisfies and is 
determined by the integral equation 

We have thus succeeded, by virtue of Theorem IV in expressing 
the response of a network to an arbitrary e.m.f. impressed at 
time i==0, by an integral equation of the same form as that 
expressing the response to a “unit e.m.f.That is to say 
we have, at least formally, extended the operational calculus 
explicitly to the case of arbitrary impressed forces. 

We now translate the foregoing into the corresponding 
operational theorem. 


Theobem VIII 
If the operational equation 

h^l/Hip) 
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expresses the response of a network to a ^^unit and if an 

arbitrary e.m,f. E impressed at time t = 0, is expressible by the 
operational equation 

E=V{p) 


or the infinite integral 



then the response x of the network to the arbitrary force is given by 
he operational equation 

Hip) 


and xH) is determined by the integral equation 


1 YSp1= 

pH ip) 



Theorem IX 
If the operational equation 

h = l/Hip) 


is reducible to the form 


h = 


m 

l+\Kip) 


where X is a real parameter, and if the auxiliary functions /=/(0 
and k = kit) are defined by the auxiliary operational equations 


f^Fip) 

k=Kip) 


then hit) is determined by the Poisson integral equation 
hit)=fit)—\^-j.J hiT)kit—T)dr. 

This theorem is of considerable practical importance in 
connection with the approximate and numerical solutions of 
operational equations when the operational equation and the 
equivalent Laplace integral equation prove refractory. In such 
cases, as will be shown later, the numerical solution of the Poisson 
integral equations can often be rapidly and accurately effected, 
and in many cases the qualitative properties of hit) can be 
deduced from it without detailed numerical solution. 
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The proof of this theorem proceeds as follows: 

By virtue of the relation /i = l/i?(p) the operational equation 




m 


can be written as 


h-{-\ 


K{v) 

H(v) 


=F{p) 




Kjy) “ 

H{p) 


A direct application of BorePs theorem or Theorem IV gives at 
once the explicit equivalent 

h{r)k(t-T)dT. 

The preceding theorems, together with the power series and 
expansion theorem solutions formulate the most important rules 
of the operational calculus, and are constantly employed in 
the solution of electrotechnical problems. On the other 
hand, the table of infinite integrals furnishes the solution of a set 
of operational equations, which are of the greatest usefulness in 
the systematic study of propagation phenomena in transmission 
systems which will engage our attention. Before taking up 
this study, however, we shall first solve a few specific prob¬ 
lems which will serve as an introduction to asymptotic and 
divergent solutions involving Heaviside^s so-called ^Tractional 
differentiation.^^ 

Problem A: Current Entering the Non-Inductive Cable.— 
The non-inductive cable is a smooth line with distributed resist¬ 
ance R and capacity C per unit length; for the present we neglect 
inductance and leakage. A consideration of cable problems 
leads to some of the most interesting questions relating to 
operational methods, particularly to questions regarding diver¬ 
gent expansions. It would seem best to allow specific problems 
to' serve as an introduction to these general questions. 

The differential equations of the cable are 





A 

dx 


I 


( 57 ) 
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where x is the distance, measured along the cable from any fixed 
point, I is the current at point x, and V the corresponding 
potential. 

Replacing d/dt by the operator p, we have 





(58) 


Eliminating, successively, V and I from these equations, we 


get 

and 

pRO-V-^V 

These equations have the general solution 

1-—— 

(59) 


(60) 

where 

7 = V pRC. 

(61) 


The term in represents the direct wave and the term in 
the reflected wave. Fi and V 2 are constants which must be 
so chosen as to satisfy the imposed boundary conditions at the 
terminals of the cable. 

For the present we shall assume that the line is infinitely long 
so that the reflected wave is absent. We shall also assume that 
a voltage E is impressed directly on the cable at x == 0: we have 
then, 

V === (62) 

(63) 

where a denotes x^RC. 

To convert these to operational equations let us suppose that 
E is a ^^unit e.m.f.^^ (zero before, unity after time t = 0). We 
have then, in operational notation 

y ^Q-Vap 


( 64 ) 
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Now suppose that x = 0 so that a = 0, in other words consider 
a point at the cable terminals. Then 
7=1 _ 

( 66 ) 

The first of these equations means that V is simply the impressed 
voltage, zero before, unity after time i5 = 0, as of course, it should 
be from physical considerations. 

Corresponding to the operational equation 


we have the integral equation 

/Cl r“ ^ 


(67) 


The solution of this is known (see formula (c) of the preceding 
table of integrals): it is 




Triit 


( 68 ) 


Heaviside arrived at this solution from considering the known 
solution of the same problem in the theory of heat flow. He 
therefore inferred that the operational equation 

has the explicit solution 

I = I/a/ 7rt. 


This is correct; we, however, have derived it directly from the 
integral equation of the problem and the known integral 


Vp Jo ' 


>-pt. 


dt 

\/ irt 


(69) 


We thus see from the foregoing that, if a “unit e.m.f.” is 
impressed on the cable terminals, the current entering the cable 
is initially infinite and dies away in accordance with the formula 
's/Cj'K'Rt The case is, of course, idealized and the infinite 
initial value of the current results from our ignoring the dis¬ 
tributed inductance of the cable, which, no matter how small, 
keeps the initial current finite, as we shall see later. 

Now let us go a step farther; suppose that in addition to 
distributed resistance JS and capacity C, the cable also has 
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distributed leakage G per unit length. The differential equations 
are now 


R.I=-±Y 

*dx 

(Cp+a)F.-|i. 


(70) 


Consequently it follows that in the operational equation for the 
current entering the cable we need only replace C-p by Cp+G. 
Therefore, when leakage is included, equation (66) is to be 
replaced by 

1 = 

where 'k=G/C. 

The corresponding integral equation is, of course, 

We shall give two solutions of this problem; first the solution of 
the integral equation, and second the typical Heaviside solution 
directly from the operational equation. 

Equation (72) may be written as 


yC+G 
R 


=J|v^x 


(71) 


\G (1+ X/p) 

■Vs v^x 


-f 




(73) 


Now suppose that J{t) is the solution of the equation 

it follows at once from Theorems (I) and (II) of the preceding 
chapter that 

/(«) = -J|(n-x£d«)/(0. (75) 


Also from formula (c) of the table of integrals and Theorem 
(Va) the solution of (74) is 


/(0 = 


'v /irt 


( 76 ) 


[C 


1 e 


r-\t rt 


\ ^ 


r 






whence 
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The integral appearing in (77) can not be evaluated in finite 
terms; it is easily expressible as a series, however, by repeatec 
integration by parts. Thus 

r r e~''^^d'\/t==2'\/i e~^^+2\ f 

Jo Vi Jo Jo 

Proceeding in this way by repeated partial integration we ge 
for the integral term of (77) 



2\t (2Xi)2 

1.3*^ 1.3.5 


(78 


The straightforward Heaviside solution is obtained by expanc 
ing the operational equation as follows: 

/=aJ|v'p+X 



1+^ 

p 


1/2 


1 + 


1 x_ 

2p 


Vp 

i_AV 

2.4 \??/ 


JL3 /xV 
'^2A:Q\p) 



Identifying -v/p with ll\^t (from known solutions of alli( 
problems) and substituting for 1/p" multiple integrations of tl 
nth order we get 


1 = 




(2X0^ 1.3(2Xi)’ 
2.3.4 ■''2.3.4.5.6 


(7 


It can be verified that this solution is convergent and equivale 
to (77). 

This problem, while simple and of minor technical intere 
will serve to introduce us to the very important and interesti 
question of asymptotic series solutions. 

An asymptotic series, for our purposes, may be defined as 
series expansion of a function, which, while divergent, may 
used for numerical computation, and which exhibits the behav 
of the function for sufficiently large values of the argume 

Let us return to equation (77). We observe that the sei 
solution (78) of the definite integral becomes increasingly lab 
ious to compute as the value of t increases. This remark app 
with even greater force to the Heaviside solution (79) on accoi 
of the alternating character of the series. Right here we h; 
an excellent example of what I regard as Heaviside's exaggera 
sense of the importance of series solutions as compared v 
definite integrals. Consider the solution in the form of (77] 
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compared with Heaviside^s series solution (79). The former is 
incomparably easier to interpret and to compute, either by 
numerical integration or by means of an integraph or planimeter. 
In fact the series (79) is practically unmanageable except for 
small values of t. 

Returning to the question of an asymptotic expansion of the 
solution (77), we observe that the definite integral appearing in 
that equation can be written as, 



provided X is positive, as it is in this case. Now the value of the 
infinite integral is known; it is 's/tt/X. Consequently 


furthermore, 




-x« 


Vt 


dt; 


(81) 


1 pc 

Ji Vt XJ( 


Vt 




1 

X V t 


-ki 


tVi 


dt. 


Integrating again by parts we get 

i'Vt 2XHVt'^2V^Jt tWt 


Continuing this process, we get 


I 



e-w 

Wt 


1 1.3_ 1.3.5 

2Xr (2Xt)~2 (2Xt)* 


_j_.. 


„ 1 .3.5.. .( 2ra-l) ~ 

(2\ty 


(-1^ 1.3.5. . .(2«+l) f” e--^ ' 

-X 2(2XF“ J, f+Vil 


(82) 


Now this series is divergent, that is, if we continue out far 
enough in the series the terms begin to increase in value without 
limit. On the other hand, if we stop with the nth term the error 
is represented by the integral term in (82) and this is less than 

.Cli)” l-3.5...( 2n-l) 

X-\7 i (2Xi)" 

Consequently the error committed in stopping with any term in the 
series is less than the value of that term. Therefore if we stop with 
the smallest term in the series, the error is less than the smallest 
term and decreases with increasing values of t. 
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We can therefore write the solution (77) as 


jxC I C I 1 1.3 1.3.5 

\ R 1 2Xi (2X0'"^(2Xi)^ 


( 84 ) 


The first term, since \=0/C, is simply -s/gJR, the d.c. admit- 
tance of the leaky cable. The divergent series shows how the 
current approaches this final steady value. 

In this particular problem no asymptotic solution is derivable 
directly from the operational equation, at least by the straight¬ 
forward Heaviside processes. Asymptotic solutions, however, 
constitute a large and important part of Heaviside's transmission 
line solutions. We shall therefore discuss next a problem for 
which Heaviside obtained both convergent and divergent series 
expansions. 

Problem B: Terminal Voltage on Cable with ‘^Unit 
Impressed on Cable through Condenser.—We now take up a 
problem for which Heaviside obtained a divergent solution, and 
which will introduce us to the theory of his divergent solutions 
and so-called ^'fractional differentiation." We suppose a “unit 
e.m.f." impressed on an infinitely long cable of distributed resist¬ 
ance R and capacity C per unit length through a condenser of 
capacity Co: required the voltage V at the cable terminals. 
The operational equation of the problem is derived as follows:— 

We know from the problem just discussed that the current 
entering the cable whose terminal voltage is 7, is, in operational 
notation 



But the current flowing into the condenser is 

Cop(l-F) 

since the voltage across the condenser is 1 — 7. Equating these 
two expressions we get 


7==. 


pCo 


which is the operational equation of the problem. 
This may be written as 


( 85 ) 


V=^ 


1 + 


W' 

->0 \ . 


c 1 
R 


l-f-Va/p 
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where 



Ip p p\p 
V \Py> 

■(^+ 1+1 

-1 + -^,+-^+* 


Now expanding this by the binomial theorem 

\p V p\p 

{aty 

" r 

{2aty 


2at (2aty 

r+X3T 


1.3.5. 


/■s/irai 


( 86 ) 


by the usual Heaviside rules of "algebrizing.” 

It is worth while verifying this from the integral equation of the 
problem. We have 


1_1 

V l+\/a lp 



(87) 


The left hand side can be written as 

1_1^ la 

p—a p—a\p 

and by the formulas and theorems given in a preceding section the 
solution can be recognized at once as;— 

V(t)=e‘‘‘-.ke‘“ r^dr. (88) 

Jo Vr 

This can also be written as 

If the definite integral of (88) is evaluated by successive partial 
integrations it will be found in agreement with the Heaviside 
solution (86). 

Now the solution (86) is in powers of i and while absolutely 
convergent becomes progressively more difficult to interpret 
and compute as the value of i increases. From (89), however, 
we can derive a divergent or asymptotic solution applicable both 
for interpretation and computation, when the value of t is suffi- 
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ciently large. As in the example discussed before, the asymptotic 
expansion results from repeated partial integrations; thus 


Jt Vr 


and finally 


Vt 
1 


2a 


r-x/1 


dr 


a \/1 

^ 1 1 . 

a\/t rVr 

, 1.3 

+ 


e~ 


a\/t 2a^t\^t 


3 r 


V. 


dr 


1~ 


1.3 


1.3.5 

-y/t ] ^ ' (2a0^ (2a0^ 


(90) 


The series (90) is divergent just as is (82) of a preceding problem 
and the error committed by stopping with the smallest term, is of 
the same character and subject to the same discussion. With this 
understanding we write the solution (89) as 


v{ty 


'\/'Kat 1 


1 1.3 

2at (2aiy 


1.3.5 

"(274'’^' 


(91) 


For large values of t {at >5) this series is accurately and rapidly 
computable. Furthermore it shows by mere inspection the 
behavior of 7(^) for large values of t, and that it ultimately 
approaches zero as 

Let us now see how Heaviside attacked this problem and how 
he arrived at a divergent solution from the operational formula. 
Returning to the operational equation (85), it can be written as 


V 'T/a 

l + 's/'p/a 


(92) 


Now expand the denominator by the binomial theorem: we geh 
formally 
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aviside’s procedure at this point was as remarkable as it was 
jcessful. He first discarded the second series in integral powers 
V as meaningless. He then indentified a/p with 
d replaced p” by d^ldt^ in the first series, getting 

\ 1 
/ •%/jTfli 

carrying out the indicated differentiation, 

1 1 I 1-3 1-3.5 

2ar(2at)^ (2at) 


H 


1+^—— + 
a dt 


(94) 


F- 


;+• 


lich agrees with (91). 

This is a typical example of a Heaviside divergent solution for 
lich he offered no explanation and no proof other than its 
actical success. His procedure in this respect is quite unsatis- 
ctory and in particular his discarding an entire series without 
planation is intellectually repugnant. We shall leave these 
lestions for the present, however; later we shall make a syste- 
atic study of his divergent solutions and rationalize them in a 
tisfactory manner. First, however, we shall take up a specific 
‘oblem for which Heaviside obtains a divergent solution without 
scar ding any terms. 

Problem C: Current Entering a Line of Distributed L, R and C. 
onsider a transmission line of distributed inductance L, 
sistance jB, and capacity C per unit length. The differential 
luations of current and voltage are 


4''- 


dx 

dx 


replacing d/dt by p, we get 

{pL+R)I-- 


Cp-V—I-L 


(95) 


(96) 


Equations (96) correspond exactly with (58) for the non- 
iductive cable, except that we must replace i? by pL+B. For 
he infinitely long line, therefore, the operational formula for 
he current entering the line is 


4 


pC 

pL+B 


Fo 


(97) 
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where 7o is the voltage at the line terminals, 
e.m.f.” we have, as our operational equation, 

If this is a “unit 

I- 1 

\pL+B 

(98) 

which can be written as 


I- P ^ 

yJ^Vl+2X/p 

(99) 


where X—2J/2L. 

The corresponding integral equation is 


f ( 100 ) 

\LVp^+2Xp Jo 

From either equation (99) or (100) and formula (p) of the table of 
integrals, we see at once that the solution is 

I = yJ~e-^Io(\t) ( 101 ) 


where Io(Xi) is the Bessel function JoCiXi), where ^ = 1- 

(The function is, however, a pure real.) 

Heaviside^s procedure, in the absence of any correlation 
between the operational equation and the infinite integral, was 
quite different. Remarking, with reference to equation (99), 
that “the suggestion to employ the binomial theorem is obvious, 
he expands it in the form 



X , 1,3 


1.3.5/X\ 
3! \pj 



( 102 ) 


aad replaces 1/p’* by i^/nl in accordance with the rule discussed 
in preceding sections. The explicit solution is then 


4 


1.3 

1—Xt+^2!)2^ 


(X<)^ 


1.3.5 

(3!)= 


(X<)»+- 


(103) 


a convergent solution in rising powers of t. As yet, however, he 
does not recognize this series as the power series expansion of 
(101), which it is. He does, however, recognize the practical 
impossibility of using it for computing for large values of i, 
and remarks “But the binomial theorem furnishes another way 
of expanding the operator (operational equation), viz. in rising 
powers of p.” Thus, returning to (99), it can be written as. 



Vp/2\ 

\/l+p/2X 


( 104 ) 
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Now expand the denominator by the binomial theorem: we get 





4X 




He now identifies \/p/2\ with l/\/ 27 r\t and replaces p” in the 
series by d^/dt^^ thus getting finally 


1 = 


IC 1 

L ’\/2TKt 


1 (1.3)^ (1.3.5)^ 

8Xr 2!(8X0" 


^ (106) 


This series solution is divergent: Heaviside recognizes it, how¬ 
ever, as the asymptotic expansion of the function and 

thus arrives at the solution 

( 101 ) 


which we have obtained from our table of integrals. 

Now the divergent expansion (106) is the well known asymp¬ 
totic expansion of the function 0”^*7o(XO; which is usually derived 
by difficult and intricate processes. The directness and simplic¬ 
ity with which Heaviside derives it is extraordinary. 

We note in this example that no integral powers of p appear 
in the divergent expansion: consequently no terms are discarded. 
Otherwise Heaviside's process is as startling and remarkable as 
in the example discussed in the preceding section. 

We shall later encounter many problems in which asymptotic 
solutions are derivable as in the preceding example. We have 
sufficient data, however, in these two typical examples to take 
up a systematic discussion of the theory of Heaviside's divergent 
solution of the operational equation.^ 

^ Quite recently an entirely different treatment of this whole problem, 
based on the Fourier Integral identity, has been developed by Dr. Norbert 
Wiener. A brief report on this work was made before the American Mathe¬ 
matical Society in May, 1925; it has since been published in Math. Annalen, 
Bd. 95, Heft 4, 1926. 



CHAPTER V 


THE THEORY OF THE ASYMPTOTIC SOLUTION OF 
OPERATIONAL EQUATIONS 

A study of Heaviside’s methods, as exemplified in the preceding 
examples and in many problems dealt with in his Electromagnetic 
Theory, Vol II, shows that they may be divided into two classes 
(I) those of which the operational equation is of the form 

l^FivWp (I^ 

and (II) those of which the operational equation is of the form 

a=<^(pVp) (ii: 

where h is an integer. 

Heaviside himself does not distinguish between the two classes, 
but employs the following rule* for obtaining asymptotic expan¬ 
sion solutions: 

If the operational equation 

h^l/Hiv) 

can be expanded in the form 

h=^aQ+aiP+a2 p^-\ -p” d- 

(ho+bip+iip^i -b&nP”H- )^/ pf (107; 

a solutionj usually divergentj is obtained by discarding the firs\ 

expansion entirely ^ except for the leading constant terms ao, replacing 
\t'p by l/VS and p'^ by d^/dt^ in the second expansion, whena 
an explicit series solution results. 

( 108 ) 

The present chapter is devoted to the mathematical theory o: 
this type of expansion solution of the operational equation, anc 
to an investigation of its practical and theoretical significance 
^ It does not appear that Heaviside makes a formal statement of this rule 
It is a generalization, however, of his procedure in a large number of specifi( 
problems. 
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• Heaviside himself gives no information which would serve us as 
a guide in informing us when the rule is applicable and when it is 
not. Consequently it becomes a matter of practical importance, 
not only to investigate the underlying mathematicar philosophy 
of the rule and to establish it on the basis of orthodox math¬ 
ematics, but also to develop if possible a criterion of its appli¬ 
cability. In this investigation we shall have recourse to the 
integral equation of the problem. 

We shall take up first the type of problem (Class I) in which the 
operational equation is 

( 110 ) 

and assume that P(p) admits of the formal power series expansion 

F('p)=bo+iiV+^2P^+h:iP^-{ - ( 111 ) 


The corresponding integral equation is 

( 112 ) 

■ Vp Jo 


We now assume the existence of an auxiliary function k(t)j 
defined and determined by the auxiliary integral equation 


Now since 


/ri' 


k(t)e~p‘dt. 


'-nt 

v Ttt 


(113) 


(114) 


it follows from (112), (113), and (114) and Borel’s Theorem, or 
Theorem IV, that 

fc(r) 


m 


= ^ f 

VttJo 


dr. 


(115) 


/O V t—T 

Now if we differentiate (113) repeatedly with respect to p and put 
p=0, it follows from the expansion (111) that 

» yi f n 

(116) 


(-!)’• 


This equation presupposes, it should be noted, the convergence of 
the infinite integrals for all values of n, and therefore imposes 
severe restrictions on k(i) and hence on F{p). We shall suppose 
that these restrictions are satisfied and discuss them later. 

Now (115) can be written as:— 

h{t) = AJ^drMr)il-r/t)-yK 


(117) 
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It can be shown that, if k{t) satisfies the restrictions underlying 
(116), the integral (117) has an asymptotic solution obtained 
as follows:—Expand the factor (1 —t/<)~^''' by the binomial 
theorem, replace the upper limit of integration by oo , and integrate 
term by term: thus 




i: 




Finally from (116) we get 

1 I 


k(ty 


*>/ 7rt 


r’" t 


10 f2 

+(2oq„ 

..(118) 

1 

, 1 ,, 1.3 , 1.3.5 , 1 

(119) 


which agrees exactly with the Heaviside rule for this case. 

The foregoing says nothing regarding the asymptotic character 
of the solution. It is easy to see qualitatively, however, that 
(IIS) and therefore (119) does represent the behavior of the 
definite integral (117) for large values of I, provided k(t) con¬ 
verges with sufficient rapidity. 

The foregoing analysis may now bo summarized in the following 
proposition; which will be called Hulk I. 

If the operational equation h = l/H(p) u reducible to the form 

h—F{p)^p 


and if F{p) admits of power seines expansion in p : thus 
F(p) =bQ+bip+hop^+ • ■ • +hnp^-i - 


so that formally, 

h=(bo+bij)+bip^+ • ■ • • )'\/p 


an explicit series solution, usually asymptotic, is obtained by replac- 
ing s/p by V-s/Tri and p’'(n integral) by d^/dt^, whence 


h{ty 
1 




(bo~bi^l 


j 1-3 I 


•Y-L 


1.3.5 


‘( 20 ^ “ ( 20 * 


•) 


provided the function k^k(t), defined by the opcralional equation 
k—F(p), and the infinite integrals 


i' 


P*k{t)di (n — 1,2,... .) 


exist 
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We shall now apply the foregoing theory to a physical problem 
discussed in the last section: namely, the current entering an 
infinitely long line of inductance L, resistance R and capacity C 
per unit length. It will be recalled (see equation (100)) that the 
integral equation of this problem is 

IC 1 


Y 




= e 


-pt 


I(t)dt 


where \ = R/2L, and that the solution is 

‘7o(XO- 


'-4 


We can derive the solution in another form appropriate for our 
purposes by writing 

_ Y— = r e-p‘Iit)dt 


Now since 


and 


1 

Vp 


Y 

"Y p-j-2X 

it follows from Borel’s theorem that 


^00 

‘J. *■’ 

^00 

ri 


dt 

"VS* 




“2Xr 


-dr. 


Vt- _ 

Now subject this definite integral (omitting the factor -y/CfL) 
to the same process applied to (117): we get 


irYilJo Yi !’• 


‘e-^^di 

f 

(20do 


2! 


■2U 


dt+' 


The infinite integrals are known and have been evaluated. Sub¬ 
stituting their values this series becomes:- 


J 


l+sl:+ 


12 . 3 ^ 5 * 


8Xi'^2!(8Xi)* 3!(8X«)® 


+ ■ 


which is in fact the well known asymptotic expansion of the 

^T^slcond Example may be worth while. Consider tlie 
of an e.m.f. impressed at time i=0 on a cable of distributed 





66 


ELECTRIC CIRCUIT THEORY 


resistance R and capacity C: required the current entering the 
cable. The required formula is^ 

Cf d 
I ttR dt^ 


/ = 


! rv 

Uo 


—T, 

Vr 


dr 


C 

rR iVi 


vr"i’';7P 


dr 


( 120 ) 


by obvious transformations. 

Asymptotic expansion of the definite integral as in the preced¬ 
ing example gives the asymptotic formula 

1 


'=~T^ 


1.3 ,1.3.5 

irRt 1 ' 


The operational formula of the problem is 
IC p 


/ = 


-4 


Vp 
V V 


Id pH-X 

c p/\ 

R 1+p/X 

lcjp_/py 

^ ' r\^ 


+ 


• • • IVp- 

i 


Applying the Heaviside Rule, wo get the asymptotic expansion 
C f 1 , 1.3 


I=- 


,1.3.5 

■jrRt\2\r {2uy-y2\ty 


which agrees with the preceding formula, derived from the 
definite integral. 

We shall now discuss a specific problem in which the Heaviside 
Rule breaks down. For example let us take the preceding prob¬ 
lem, and replace the applied e.m.f. by sin cct The formula 
corresponding to (120) is now 

'c r 




cos COT 


ttEJo \/ f 


dr. 


( 121 ) 


If we now attempt to expand the definite integral of (121) in the 
same way as that of (120), we find that the process breaks dowr 
because each component of the infinite integral is now itselJ 
infinite. In fact no asymptotic solution of this problem exists 
Let us, however, start with the operational formula; since 

£0 


/: 


sin o)idt = 




1 The derivation of the formulas in this problem is left as an exercise for thi 
reader. 
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it is 


J = 


C up 


R p^+w 


Wp- 


Now expand this in accordance with the Heaviside Rule: we get, 
operationally^ 


^-4 

and explicitly 


7 = 



Vp 



1.3.6 

{2o,ty 


+ • 


which is quite incorrect.^ The incorrectness of the result will 
be evident when we remember that the final value of the current 
is the steady state current in response to sin o)tj or 



(cos wi+sin oit). 


( 122 ) 


This result can be derived directly from (121) by writing it as 
T I ^ ^ f • ^ r^sin cot 

[ Jo Jo Vi 

If the time is made indefinitely great the upper limits of the 
integrals may be replaced by infinity. The infinite integrals 
are known: substitution of their known values gives (122). 

This example illustrates the care which must be used in apply¬ 
ing Heaviside’s rules for obtaining divergent solutions and the 
importance of having a method of checking the correctness of 
his processes and results. 

We now take up the discussion of the asymptotic expansion 
solution of operational equations of the type 

h = <^>(p^'\/p) (k integral). (123) 

The Heaviside Kule for this type of operational equation may 
be formulated as follows (Rule II): 

If the operational equation h—l/II{p) is reducible to the form 

h=<t>(vWv) 

and if ^ admits of power series expansion in the argument^ thus 

/?. = ao+o^iP^ a/29+o^2P^^‘^^+c& 3P®^'^’’ '\/p + * * * 

^ While this series is incorrect as an asymptotic expansion of the current it 
has important significance, as we shall see, in connection with the building up 
of alternating currents. 
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a s&ries solution, usually divergent and asymptotic, is obtained by dis¬ 
carding integral powers of p, and writing 

/i = ao+(ai2)*’+a3P*''+^+a6?'“+*H-)Vp- 


The explicit series solution then results from replacing \/p by 
\Is/ret, and p” by d^fdt^, whence 


/ d'‘ 

Jee^aQ -^( ai j,,^~\~az I ^5 ' 




'■dt'‘ ' ^dt^'‘+^ ' “"“dP+z 


Js/rrt 


~ao- 


(-lyr 1.3. ■ .(2^- 1)1.3...( 6fe+l) 
Vrrt \ ^ 


/2t/ (2<)“+i^ 

The theory of this series solution will be based on the follow- 
ing proposition, deducible from the identity I —di — 1 /\/p. 

Jo V TTif 

If the function F{p) of the integral equation 




^’(p) = Jo me-^^di 

approaches l/\^p as p approaches zero^ thenf{t) ultimately behaves 
as l/V^ri: that is^ if P(p)->l/\/p as p—>0, ikenf{t)^l/\/Trtas 
t—> CO f provided that f(i) converges to zero, and contains no term or 
factor which is ultimately oscillatory. 

To illustrate what this condition means suppose that 

a , 6 cos 


then 


X 


vs 

w 

f(t)e-p‘dt—^a/-\/p as p~^0. 


and the oscillatory term, in f(t) converges to a higher order. The 
presence of such oscillatory terms vitiate, therefore, the Heaviside 
Rule: in the following discussion we shall assume that they are 
absent. 

We are now prepared to discuss the operational equation 
h = 4>(p’‘ -y/ p) 


and for convenience shall assume that Ic = 0 so that the operational 
equation becomes 

h = ct>(\/p) 


of which the 


corresponding or equivalent integral equation is 


^<^(Vp) = 

Jr 


03 

h(t)e~^^dt. 


(123b) 
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We assume that <^»('\/p) admits of formal power series expansion 
in the argument: thus 

^( Vp ) =QJo+ai Vp +<^2 2?+<^3P'\/p+a4p^d- 

without, however, implying anything regarding the convergence 
of this expansion. 

We now introduce the series of auxiliary functions, 
gr 3 , . . . defined by the following scheme 
g{t)=^h{t)-aQ 



Successive substitutions in the integral equation (123b) "“and 
repeated differentiations with respect to lead to the set of 
formulas, 


j g{t) h_ as p 0 

Jo V p 

r dt^ -^^as p-^0 

Jo 2v p 

J tg 2 it)e-p^ dt^~ as p 0 

J tgz{t)e-^^dt‘^-^^^- P'->0 


(123d) 


Now assuming that hit) satisfies the restrictions stated in the 
preceding proposition, it follows from that proposition, that 
g{t)'^ai/'\/'Kt as 00 

C123e) 

/ V 1.3.5 a^ , 
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From the set of equations (123d) and (123e) it follows by 
successive substitutions that 

^ /a 


h{ty 




1 , 1.3 


1.3.5 , 

-U7-(2^y3+- 


which agrees with the series gotten by applying the Heaviside 
Rule. 

The defect of this derivation, which, however, appears to be 
inherent, is that it requires us to know or assume at the outset 
that 'h{t) satisfies the required restrictions. Consequently an 
automatic application of the Heaviside Rule may or may not 
give correct results. On the other hand if we know that an 
expansion solution in inverse fractional powers of t exists, the 
Heaviside Rule gives the series with extraordinary directness 
and simplicity. 

The type of expansion solution just discussed will now be 
illustrated by some specific problems. The first problem is that 
of the propagated voltage in the non-inductive cable in response 
to a ^^unit e.m.f.^' It will be recalled that in a preceding chapter 
we derived the operational formula 

(124) 

where a-x^RCj for the voltage at distance x from the terminal 
of a non-inductive cable of distributed resistance R and capacity 
C, in response to a “unit e.m.f.^’ impressed at point a; = 0. Heavi- 
side^s solution of this operational equation proceeds as follows: 

Expansion of the exponential function in the usual power 
series gives 


7-1- 


ap .ap 
1! ^2!' 


aps/ap .{apY 
8! 4! “ 


which may be rearranged as 


7 = 1- 


^+3!^^ 5! + 


•\/ ap 




ap {afY 

2!"’' 4! 6! 


+ • 


(125) 


Heaviside now discards the series in integral powers of p entirely, 
replaces Vp by l/-\/wi and p" by in the first series, and 

thus gets 


F = l- 


a d 




1 

'3! 


H- 


1.3 

5! 


a 

wt 

a 

2t 


7! 


+ • 


(126) 
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or 


7 = 1- 




+ 


5.2!V4i 


fa 

'7.3l[4t 


+ • 


(127) 


This solution is correct, as will be shown subsequently. 

A rather remarkable feature of this solution—a point on 
which Heaviside makes no comment—^is that it is absolutely 
convergent. In other words, a process of expansion which in 
other problems leads to a divergent or asymptotic solution, here 
results in a convergent series expansion. 

To verify this solution we start with the corresponding integral 
equation of the problem 

► 00 




P 


■f 


V(t)e'-P^di. 


(128) 


It follows from this formula and theorem (II) that 

V(l)= fVodi 
Jo 

where </>(«) is determined by the integral equation 
^ 00 

g-Vip= 


Now from formula (f) of the table of integrals 



2\irJo 


whence 

♦(f)-2^ 

ja 





and finally 


. (Itj where t* = U/a, 

^Jo Vr 

(129) 

To convert this to the form of (127) we write 



II 

(‘“e-Vr 1 

do rs/r Virjr rjr 

(130) 


The value of the infinite integral is known to be \/^ so that 


7 = 1- 


1 r<^e-Hr 

r-\/T 


dr. 


(131) 


Now in the integral term of (131) expand in the usual 
exponential power series and then integrate term by term; the 
series solution (127) results. This series, while absolutely con- 
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vergent, is difficult to compute for small values of ^; an asymptotic 
expansion, which can be employed for computation for small 
values of t is gotten as follows: 

Write (129) as 




2\/'jrJo Vr 


Repeated partial integrations of this type lead to the series 




(132) 


It is interesting to note, in passing, that an asymptotic solution 
of this type does not appear to be directly deducible from the 
operational equation. We observe also that, in this problem, 
the series in inverse powers of t is convergent while the series 
in ascending powers of i is divergent: the converse is the case in 
the problems discussed previously. 

A second specific problem may be stated as follows: 

Let a “unit e.m.f.” be impressed on an infinitely long non- 
inductive cable of distributed resistance R and capacity C per 
unit length through a terminal resistance i2o: required the voltage 
V on the cable terminals. The formulation of the operational 
equation of this problem is very simple. It will be recalled 
that the operational formula for the current entering the cable 
with terminal voltage V is But the current is clearly 

also equal to (1 —y')/i2o*. equating these expressions we get 


whence 

7=-7-J- (133) 

v^pA+i 

where This is the required operational 

formula. 

To derive the Heaviside divergent expansion, expand (133) by 
the binomial theorem: thus 


7=1- VpA+Ci’A) - (23 A)“^^+i: • 

= 1 —(l+p/X4-(p/X)^H-• •)-\/pA (134) 

+(pA+(pA)*+(pA)'‘+---). 
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Discard the second series in integral powers of p ; replace V "p 
by I/VtS and p^ by d'^^/dP in the first series, thus getting 


"‘-vU‘-2S+(2W-) 


which is the asymptotic solution of the problem. 

To verify this solution we shall consider the 
operational equation 

h— y (n integral) 

p”Vp + l 


more general 


(137) 


a form of equation to which a number of fairly important 
problems is reducible. (The parameter X of equation (133) can 
be eliminated from explicit consideration by means of theorem 
VI.) 

Multiplying numerator and denominator of equation (137) 
by p^\/p — l, it becomes 


/i = 


_ a/p — 1 




pn 

2n (-1_ 


iVp- 


p 




(138) 


and by direct partial fraction (ixpansion, this is equivalent to 


n+l 


2 n 

2n+l-^p—p,n 2n+i^^p—p„ 

0 nie=i 0 


Vv 


where 

Write, for convenience. 


,_2TO7I^ 

= (m = 0,b2. . ,2a)- 


2n 


and consider the operational equation 


2n+lVp-p 


71-f 1 

P-P, 


(139) 


(140) 


By the rules of the operational calculus, fully discussed in 
preceding chapters, the solution of this is 




( 141 ) 
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We have now to distinguish two cases: (1) when the real pan 
of prn is positive, and (2) when the real part is negative. 

Taking up case (1) first, the preceding can be written 




1 






«4-l 


P 

V' 


r 

Jo \/r 

fpn+l 

iV^Jo Vr 




(142) 


(1431 


(144) 


Repeated integration by parts of the definite integral leads to an 
asymptotic series, identical with that obtained by applying the 
Heaviside Rule to the operational equation (140). 

If, on the other hand, the real part of pm is negative, we write 
(141) as 




1 

2ra-f-l 



'v/tt 



(145) 


The term ultimately dies away, and the definite integral 
can be expanded asymptotically in accordance with the theory 
discussed under Rule 1, again leading to an asymptotic series 
identical with that given by direct application of the Heaviside 
Rule to the operational equation. 

Consequently since the operational equation in can be 
asymptotically expanded by means of the Heaviside Rule, the 
operational equation in h — ^hm is similarly asymptotically 
expandible, and the Heaviside Rule is verified for equation (137). 

We have now covered, more or less completely, the theoretical 
rules and principles of the operational calculus in so far as they 
can be formulated in general terms. We shall now apply these 
principles and rules to the solution of important technical 
problems relating to the propagation of current and voltage 
along lines. In doing so, while we shall take advantage of our 
table of integrals with the corresponding solutions of the oper¬ 
ational equation, we shall also sketch Heaviside^s own methods 
of solution. 

We shall close this discussion of divergent and asymptotic 
expansions with a general expansion solution of considerable 
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theoretical and practical importance in the problem of the 
building-up of alternating currents. It will be recalled from 
Theorem VIII that the response of a network of generalized 
operational impedance H(p) to an e.m.f. E(t) impressed at time 
i = 0 is given by the operational formula 


H{p) 


where E=V(p) is the operational equation of the applied e.m.f.: 
that is, analytically 


1 

V 


F(P): 


i' 




Now suppose that the impressed e.m.f. is sin ojt: then by formula 
(h) of the table of integrals 




op 


and denoting x by Xg 


Xa 


'^p‘^+oHi(py 


(146) 


(147) 


If, on the other hand, the impressed e.m.f. is cos cot, then by 
formula (i) 

r)2 

(148) 


p^+co‘^ 




and 


x=x, 


p® 1 


p^+co^H(p) 


(149) 


Now let us consider the operational expansion suggested by the 
Heaviside processes: 


and 



(150) 



( 151 ) 
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Now’let us identify l/H{v) with h{t) and replace by d^/dt^: 


we get 

x,= 

ld_ 

o)dt 

1 d^ 1 _' 

co^dt^ co^di^ j 

• Ht) 

(152) 

and 

Xc— ‘ 

1 d^ 
co^dt^' 

i 

+ 

1 3 

1 

■hit). 

(153) 


We have now to inquire into the significance of equations (152) 
and (153), derived from the operational equations for the response 
of the system to an e.mi. sin co^ and cos respectively, impressed 
at time i=0. From the mode of derivation of these expansions 
from the operational equations it might be inferred that they are 
the divergent or asymptotic expansions of the operational 
equations (147) and (149). This would certainly not be an 
unreasonable inference in the light of the Heaviside expansions 
we have just been considering. This inference is, however, not 
correct: on the other hand, the series (152) and (153) have a 
definite physical significance, as we shall now show from the 
expheit equations of the problem. 

By equation (31), the explicit equation for given operation¬ 
ally by (147), is 



sin — 


T)dT — ^ Bhxo){t—T)h\T)dr+h{0) sin o)t 


(154) 


where h\t) =d/dt h(t). By a well known trigonometric formula, 
this is 


Xs = Bin oot r cos o)t.h\t)dI — cos o)t fsin a)t.h'(t)di+h(())Bm cot. 
Jo Jo 


Writing 


this becomes 
= sin 


I dt= I dt— I dt 

Jo Jo Jt 


cos coLh\t)dt — Q.o^ o)t 


^00 

j sin ut.h' 


(i)di 


X CO 

sin cu(i-T)A'(r)dr. (155) 


The first three terms are simply the steady state response to 
the impressed e.m.f. sin cot: that is, they represent the ultimate 
steady state value of when the transient oscillations have 
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died away. The last term, which we shall denote by Tt 
sents the transient oscillations which are set up when th 
is applied. Thus 

X oo 

sin o)(t--r)h\r)dr. 

Now from (156) 

I /i'(r).d.cos a)(T —i) 

and integrating by parts 

Repeating the process of partial integration, we get: 

Repeating the process again 

This process can be repeated indefinitely, and we get 



1 1 

i sma,(r-i)^A(r)d- 


The series expansion (159), except for the remainder 
identical with the series expansion (152) derived direct 
the operational equation. This series may be either cor 
or divergent, depending on the frequency o}/2v and the c] 
of the indicial admittance function h{t). In the iir 
problems of the building-up of alternating currents ii 
and lines we shall see that, even when divergent, the 
of an asymptotic character and can be employed for comp 

We thus arrive at the following theorem: 

If an e.m.f. sin ut is impressed at time i{=0 on a net 
system of generalized indicial admittance h{i), and if the i 
distortion, T„ is defined as the instantaneous difference 
the actual response of the system and the steady state r 
then T, can be expressed as the series, 
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with a remainder term 


If the impressed e.m.f. is cos oit^ the corresponding series for the 
transient distortion, Tcy is 


\a)2di52 *’ 

with a remainder term 


(161) 


(^l)n ^27H-1 

■“2n““" COS CO(r if)^-^2n+l ('^) 

The second part of this theorem, relating to the transient 
distortion, Tc, in response to an e.m.f. cos coif, is derived from 
formula (31) by processes precisely analagous to those employed 
above in deriving the series expansion for Ts. The derivation 
will be left to the reader. 

To summarize the preceding discussion of the divergent 
solution of operational equations, it may be said that the theory 
is as yet rather unsatisfactory. To the physicist it is unsatis¬ 
factory because he requires an automatic rule giving a correct 
asymptotic expansion by purely algebraic operations without 
investigations of remainder terms or auxiliary'functions. Fur¬ 
thermore, the precise sense in which the expansion asymptotically 
represents the solution cannot be stated in general, but requires 
an independent investigation in the case of each individual 
problem. 

On the other hand, when an asymptotic expansion is known 
to exist, the Heaviside Rule finds this expansion with incom¬ 
parable directness and simplicity, the problem of justifying the 
expansion being a purely mathematical one, which usually need 
not trouble the physicist. Furthermore, on the purely mathe¬ 
matical side, the Heaviside Rule is of large interest and should 
lead to interesting developments in the theory of asymptotic 
expansions. 


APPENDIX TO CHAPTER V 


It will have been observed that the asymptotic expansions 
derivable from the operational equation by the direct applica¬ 
tion of the Heaviside Rules, as discussed in the preceding chapter, 
are all of the form: 
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In other words, they are all simple series in inverse fractional 
powers of t. In addition to this limitation, the direct applica¬ 
tion of the Heaviside Rules, as we have seen, does not always 
lead to correct results. In this appendix we shall briefly discuss 
some extensions and modifications of the direct Heaviside method 
whereby the limitations, inherent in the direct application of the 
Heaviside Rules, may often be successfully avoided. 

The asymptotic solution of operational equations is of great 
importance from the standpoint of both pure and applied mathe¬ 
matics. Its importance in practical applications to technical 
and physical problems does not require emphasis, while from the 

mathematical standpoint it not only is of value in the solution 

^ 00 

of infinite integrals of the type I but promises to be a 

Jo 

useful tool in the general theory of asymptotic expansions. 
I. We start with the operational equation 


h = 


Fiv) 

G{p) 


( 1 ) 


and assume that F(p) admits of the formal scries expansion 

F(p) =ho+bip+h2P^+h:iP^-\ - (2) 

It is evident that equation (1) is a generalization of (110) of the 
preceding chapter, into which it degenerates when 1/G(p) is 
put equal to '\/p. A consideration of (1) will therefore lead to a 
formal generalization of Rule I. 

Let the function g~g(t) be defined by the operational equation 

g=^l/a(p) (3) 

and an auxiliary function h = h{p) by the operational equation 


h — p*F{p). 


(4) 


Assuming that these functions exist, wo have by BorePs theorem 



fc(r)( 7 (c—r)dr 


(5) 


Now returning to (1) let us write, by virtue of the expansion 

( 2 ), 

h = (ao+aipH-a2 


( 6 ) 
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Now, guided by the usual Heaviside procedure, let us replace 
l/(?(p) by g{t), in accordance with (3), and replace p" by 
in the series factor. The right hand side of (6) then becomes 

(ao+ai^+a2^+'• OffCO . (7) 

We have now to inquire into the significance of this expansion. 

Expand the function g{i—r) as a power series in t, and substi¬ 
tute in (6): we get 

h{t) =S'(0j^ 

' fj/c(r)dr (8) 

Jo 

+ ^\{T)Rn{t-T)dr 

where 2?^(i5—r) is the remainder in the expansion of g{t--r). 

We now suppose that h{t) converges in such a way that the 
integrals 

J^“^,/c(r)dT (n=1,2,3.. n) 

exist. On this assumption it follows from the integral equation 

noo 

J k(t)e^^^dt=F(p) 

and the series expansion (2), that 


Now writing 


Jo Jo Jl 


equation (8) becomes by virtue of (9), 

hit) = (bo+h~+h^,+ - ■ +bn^g{t) 
-fj* h{T)Rn{t—r)dT 
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The first series of (10) is identical with (7), while the rest of 
the right hand side of (10) represent remainder terms. The 
question as to whether the series (7) represents the function 
h{t), asymptotically or otherwise, can only be answered by an 
independent investigation of the remainder terms in each specific 
problem. The expansion (10) therefore formulates a general¬ 
ization of the problem underlying Rule I, into which it degener¬ 
ates when g(t) = l/\/Tri. The fact that no general statement 
can be made as to the way in which the series (7) represents 
h{l) indicates the care which must be taken in the interpretation 
of the operators. Often, however, the expansion (7) is valid 
and quite useful: in addition it possesses considerable theoretical 
interest. 

II. If the variable x(t) is defined by the definite integral 

a;(<) = 4 (2-1) 

dtjQ 


it is known that, provided the real part of X is not negative, it 
may be written as 

( 2 . 2 ) 

Furthermore, in a number of important cases the component 
y{t) admits of asymptotic expansion. On the other hand, the 
operational equation 

P _L 

corresponding to (2.1) does not admit of direct asymptotic 
expansion and the direct application of the Heaviside Rule 
is invalid. We can, however, set up an operational equation 
for y{t) which leads to the required asymptotic expansion, where 
such an expansion exists. The procedure is as follows: 

From (2.2) and (2.3) we have the usual integral equation: 


x = - 


(2.3) 


/: 


dte-p'' 


yit)+i 




Integrating the second term in the integral and rearranging 
gives 


I 




1 _ 1 

,0 ■p"”-xL//(p') F(X) 

and the corresponding operational equation is 


V’- 


■p-X 


' 1 
H(p)‘ 


1 

H(X) 


(2.4) 
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As an example, let \=iu and = The problem is 

then that discussed on page 66 and formulated by equation 
(121), concerning which it was remarked that the straight¬ 
forward procedure gave incorrect results. In this case 

i/H(v)==Vp; im-K) = V^ 

and the operational equation is 

= _ p _ I 

'\/p-h’\/io) '\/icx)\ 'v/ io)/ 

or, expanding in powers of ^/ p, 

\ L pVp I f _fy y_\ 

■\/m\ ■\/iai (^'“) (iu)-\/ic>i 

Dropping the integral powers of p, according to Rule II (p. 67) 
and replacing Vp by l/Viri! and p” by d’^/df^ gives 
■ fid 1 11 


which is correct. 

Note that, if in the last equation we generalize l/\Ar^ 
h(t), the series becomes the equivalent of (160) and (161) for the 
transient distortion in the case of sinusoidal impressed forces. 

III. In a number of important problems, the function h, 
defined by the operational equation ^ = l/J/(p), has an asymp¬ 
totic expansion of the form 

h'-^-^r^^{ao+ai/t+a2/t^-\ -) (3.1) 

In such cases the direct application of the Heaviside rules 
breaks down, due to the presence of the exponential factor. If, 
however, the form of the asymptotic expansion is known, the 
following procedure enables us to derive the asymptotic expansion. 
by operational methods. 

Let us write 

and substitute in the integral equation: we get 


+ 


1 ^ 
?df-' 


1 d-* 
co^di^^’ 



(2.5) 



(3.2) 
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and, by Theorem Va, 

r® 1 

J. 

whence, operationally, 




V 


(p+X)f/(p+X) 


(3.3) 

(3.4) 


If the assumed form of expansion exists, then g must admit of 
expansion in inverse fractional powers of i, and this expansion 
must be derivable from the operational equation by the direct 
application of the Heaviside rules. 

Example 

As an interesting example of the preceding we shall derive 
the usual asymptotic expansion of the Bessel function J’o(XO 
from the operational equation (n) and the corresponding integral 
equation 

(3.5) 

Jo 


It is known that the function /o(Xi) can be represented as 

u{t) cos Xi!+i;(iS) sin \t (3.6) 

and that u and v admit of asymptotic expansions in inverse 
fractional powers of t. 

Write (3.6) as 

[u(t) —w(t)] + [u{t) +w{t)] (3.7) 


and substitute in (3.5): we get 


1 r 
2J0 

+i/. 


dt 6“*’* e*’'* [u{t) —iv{t)] 

dt e~^* [w(i) +iv(t)] = —~ > 

0 Vp^+X® 


(3.8) 


whence, replacing p—i\ by q and then replacing q by p, we get: 


dt e~^‘ [uit) 


1 r“ 

2J0 

1 r“ 1 

+w(0]= 


(3.9) 


-\-i2\p 


Now the function 




is oscillatory and does not admit of asymptotic expansion in 
inverse fractional powers of f. Consequently if the right hand 
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side of (3.9) is expanded in fractional powers of p, we shall get 
therefrom by the procedure discussed in Chapter V the asymp¬ 
totic expansion of the function 

^[uit)-w{t)] 


That is to say, in an asymptotic sense only, we can set up the 
operational equation 




2p 


\/ p^-\-i2\p 
2 

\/i2\ 


i^- 

2X 


Vp 


-‘i |2 

Wx 


/- 

1 -#' 

h 

L 2Xj 


-H 


■s/p 


(3.10) 


Thefactore suggests that m—w be replaced bye ^Q), 

whence 


P-iQ = 



Vp 


(3.11) 


Expanding the right hand side by the Binomial Theorem, replac¬ 
ing p" by (PfdP and \/p by l/s/irt, we get: 

/2Ti_- 1 _ 

'VtW "sXi 2!(8M)^^ 


P-iQ- 


whence 


+i 


,1132.52 , 12.32.52.72 




2 

irXi 


3!(8X0' 4!(8X«)' 


1*.3" 


2!(8X<)*“ 


Q'- 



1 , 

“ 

/ TrXi 

C>5 

00 

CO 

00 



and finally 


(3.12) 


(3.13) 


J’o(XO~^’ cos ^Xf—sin^Xi—(3.14) 


which is the usual asymptotic expansion. 

The foregoing discussion is tentative and far from complete. 
It does, however, suggest that important developments of the 
theory of asymptotic expansion may result from a study of the 
operational equation or rather its integral equation equivalent. 




CHAPTER VI 


PROPAGATION OF CURRENT AND VOLTAGE IN THE 
NON-INDUCTIVE CABLE 

The principal practical applications of the operational calculus 
in electrotechnics are to the theory of the propagation of current 
and voltage along transmission systems. Of such transmission 
systems the simplest is the non-inductive cable. The theory of 
the non-inductive cable is not only of great historic interest, 
relating as it does to Kelvin’s early work on the possibility of 
transatlantic telegraphy, but is also of very considerable practical 
importance today, and serves as a basis for the theory of submarine 
telegraphy over long distances. We shall therefore consider the 
propagation phenomena in the non-inductive cable in some detail. 

The propagation phenomena in any type of transmission 
system are isolated and exhibited in the clearest possible manner 
when we confine attention to the infinitely long line, with voltage 
applied directly to the line terminals. Furthermore, as we 
shall see later, the solution for the infinitely long line is funda¬ 
mental and can be extended to the more practical case of the 
finite line with terminal impedances. We therefore, in this 
chapter, shall confine our attention to the case of the infinitely 
long cable with voltage applied directly to the cable terminals. 

Consider a cable of distributed resistance R and capacity C 
per unit length, extending from a; = 0 along the positive x axis. 
From a previous chapter (see equations (64) and (65), we are in 
possession of the operational equations of voltage and current; 
they are, for the infinitely long line, 

7 = e-v'orp7„^ ^ (102) 

I= 7o = Vo, (163) 

where a^x^BC, and Fo is the terminal cable voltage at x = 0. 
Let us now assume that the terminal voltage Fo is a ^^unit 
e.m.f.;’^ then ^ 

7 = e-Vap^ (104) 

I=-~V^pe~^o‘P. (165) 

85 
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The solution of (164) for V was considered in some detail in the 
preceding chapter; it is, by (120) 


F = 




(166) 


where r = 4:tja = ^ilx^RO. Series expansions of this solution were 
also given. Another equivalent form is, by (131) 

2 /*!/ a/t 

F = (167) 

VttJo 

This last form, recognizable also from inspection of the series 
expansion (132), is useful because the integral term is what is 
called the error function and has been completely computed and 
tabulated. 

Before discussing these formulas and the light they throw on 
propagation phenomena in the non-inductive cable, we shall 
derive the solution for the current. A very simple way of doing 
this is to make use of the differential equation (57) 


Now from (166) and the relation 

d _ dr d 
dx dxdr 

we get 

dx ^ T'\/ T X^RCj 

2 

X's/t 's/ r 

whence 

y 2 6 “^/^ I "o' 


It is worthwhile verifying the formula by direct solution from 
the operational equation (165). From formula (g) of the table 
of integrals, we have 

jc 

Vvt^lR 
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Comparison with the operational equation shows that they are 
identical, within a constant factor provided we put X = 
Consequently the solution of (165) is 


/ = . ~ g-l/r 


which agrees with (168). This, it may be remarked, is an excel¬ 
lent example of the utility of the table of integrals in solving 
operational equations. 

This formula is easily calculated for large values of i by 
expanding the exponential function; it is 


2 J 

RX TTT 



The propagation phenomena of the non-inductive cable are 
therefore determined by the pair of equations 


f 

\/irJo 


tq-I/t 2 

— ar == 1—7 

rVT V' 


-r 


1/ Vt 


'^dr 


and 


where T=^4d/a^ 


1 = 
At 


-X/r 


_ G 

's/txR "s/r \7rRt 




-1/r 


(169) 


(170) 


Now an important feature of these formulas is that the voltage 

4 

at point a: is a function only of t; that is, of 4< divided hy the 


total resistance and capacity of the cable from a;=0 to a:=a;. The 
same statement holds for the/om of the current wave; its mag¬ 
nitude, however, is inversely proportional to xR, or the total 
resistance of the cable up to point x. Consequently a single 
curve, with proper time scale serves to give the voltage wave at 
any point on the cable. Similarly a single curve, with proper 
time and amplitude scales, serves to depict the current wave at 
any distance from the cable terminals. These curves are given 
in Figs. 3 and 4. 

Referring to the curve depicting the current wave, we observe 
that it is finite for all values of 4>0; consequently, in the ideal 
cable, the velocity of propagation is infinite. This is a con¬ 
sequence, of course, of the fact that the distributed inductance 
of the cable is neglected. Actually, of course, the velocity of 
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propagation cannot exceed the velocity of light. The error, 
however, in neglecting the inductance in the case of long cables 
is appreciable only near the head of the wave provided we confine 
attention to d.c. or low frequency voltages. This point will be 
discussed and explained more fully in connection with the 
transmission line. 



Fig. 3 . —Current in non-inductive cable ((? — 0 ) unit e.m.f. applied. 



0 t.O 2.0 3.0 AQ 5.0 6.0 7.0 tQ 60 10,0 

Fig. 4. —Voltage in non-inductive cable (G^~0) unit e.m.f. applied. 


The current, while finite, is negligibly small until r reaches th( 
value 0.2. In the neighborhood of this point it begins to buik 
up rapidly; reaches at t = 2 its maximum value 

2 2 

and then begins to decrease, ultimately dying away in accordance 
with the formula 
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Its subsidence to its final zero value is very slow; for example^ 
when r = 100 its value is still 

Turning to the voltage curve, Fig. 4, we see that it is negligibly 
small until t reaches the value 0.25, at which point it begins to 
build up. Its maximum rate of building up occurs when r = 2/3, 



after which it builds up more and more slowly. Its approach to 
its final steady value is in accordance with the formula 

y = l- ^ j-i A_\ 

V'^ttA 3j+2!5r== I 

Even, therefore, when r is as great as 100, V differs sensibly 
from its ultimate value, unity, its value being 0.8876. 

x^RC 

Since the actual time is —it follows that the speed of 

building up is inversely proportional to the square of the length 
of the cable. 

The power curve VT is given in Fig. 5. V-I is the rate at 
which energy is being transmitted past the point x of the cable. 

The fact that the form of the current and voltage waves 
depends only on it/x^RC is at the basis of Kelvin’s famous “KR” 
law, long applied to cable telegraphy and sometimes incorrectly 
applied to telephony. When the first transatlantic telegraph 
cable was under consideration, Kelvin attacked the problem of 
propagation along the non-inductive cable and arrived at for- 
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mulas equivalent to (169) and (170). From these formulas he 
announced the law that the speedof the cable, i.e., the number 
of signals transmissible per unit time, is inversely proportional 
to the product of the total capacity and total resistance of the 
cable (KR in the English notation). To see just what this 
means requires a little digression into the elementary theory of 
telegraph transmission. 

Telegraph signals are transmitted in code by means of ^^dots^' 
and “dashes,” The “dot” is the signal which results when a 
battery is impressed on the cable for a definite interval of time, 
after which the cable is short circuited. A “dash” is the same 
except that the time interval during which the battery is con¬ 
nected to the cable is increased. The “dots” and “dashes” are 
separated by intervals, called “spaces,” during which the cable 
is short circuited. Now when the cable is short-circuited we may 
imagine a negative battery impressed on the cable in series with 
the original battery. Consequently the current in the cable, 
corresponding to a signal composed of a series of dots, dashes and 
spaces, will be represented by a series of the form 

J(0-J(^-^l)+/(^-^2)~I(^-^3)+^(^-M“• * • (171) 

where, in the cable under consideration, I{t) is given by (168). 
ti is the duration of the first impulse, of the first space, 

Iz'-h of the second impulse, etc. 

Now by (168) 

2 2 .. 

T is, of course, 4:t/x^CR-4:t/KR (in the English notation). Now 
suppose that 

4^2 4, 

Then, the signal can be written as 

■~^-r--{<t>(r)~tl>(T—ri) + cl>(T—Ti) -}. (172) 

xR-Vir 

Now if the relative time intervals n, n . . . are kept constant 
(as the length of the cable is varied), the actual time intervals 
ti, h . . . are proportional to x^CE or to KR, and the wave 
form of the total signal is independent of KR, when referred to 
the relative time scale r. Hence, if T is the total time of the 
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signal, T is proportional to x^CR (or to KR), That is to say, if 
the duration of the component dots, dashes, spaces of the signal 
are proportional to the of the cable, the wave form of the 

received signal, referred to the r time scale, is invariable, and the 
total time required to transmit the signal is proportional to 
the of the cable. Now the maximum theoretical speed of 

transmission on the cable is limited by the requirement that the 
received signal shall bear a recognizable likeness to the original 
system of dots and dashes: in other words there is a maximum 
allowable departure in wave form between received and trans¬ 
mitted signals. If, therefore, the actual speeds of two cables 



Fig. 6.—Elementary telegraph signals in non-inductive cable. 


are inversely proportional to their ^^KRs/^ the wave form will be 
the same. This establishes Kelvin^s KR ” law. As a corollary, 
if the length of the cable is doubled the speed of signaling is 
reduced to one-quarter, assuming the same definition of signals. 

The foregoing will be somewhat clearer, perhaps, if we refer to 
curves, 4, 5, 6 , 7 of Fig. 6 which illustrate the distortion suffered 
by elementary dot signals in cable transmission. Curve 4 shows 
the dot signal produced by a unit battery applied to the cable 

x^BC 

terminals for a time interval f = 2 ^ while curves 5, 6 and 7 


are the corresponding dot signals when the battery is applied 


for the time intervals 


x^RC Ix^C ,lx^RC 

-^, - 2 - 4 -and 4 --^— 


Any further 


decrease in the duration of the impressed dot, beyond that 
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shown in curve 7, does not affect the shape of the transmitted dot, 
which means that the cable speed has reached its theoretical 
maximum. These curves, it should be observed, can be inter¬ 
preted in two ways. First, we can regard the length x of the 
cable as fixed and the duration of the impressed dot as varied. 
On the other hand, we can regard the actual duration of the 
impressed dot as constant and the length of the cable as varied. 
From the latter standpoint the curves illustrate the progressive 
distortion of the signal as it is transmitted along the cable. 

The dot signal of relative duration T can be written as 
D=/(t), r<T 
= r>T 

and the second expression can be expanded in a Taylor^s series, 
giving 

If T is sufiSciently short this becomes 

(173) 

Hence when the dot signal is of sufficiently short relative duration 
jT, the wave shape of the received signal is constant, /^(r), and 
its amplitude is proportional to the relative duration of the dot. 

This can be generalized for any type of transmission system: 
Let the dot signal be produced by an e.m.f./(0 of actual duration 
T. Then the received dot signal, by formula (31), is 

t<T 

= t>T. 

For t>T this becomes 

D^r{t) f /(r)dr-~r'(0 fV(T)dr+... 

Jo Jo 

and for sufficiently short duration T, we* have approximately, 

D = I'{t)^J{r)dr. (174) 

lienee for a sufficiently short duration of the impressed e.m.f. 
the received dot signal is of constant wave form, independent of 
the shape of the impressed e.m.f., and its amplitude is propor¬ 
tional to the time integral of the impressed e.m.f. These princi¬ 
ples are of considerable practical importance in telegraphy. 
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The leaky cable, that is, a cable with distributed leakage 
conductance G in addition to resistance R and capacity C, is of 
some interest. The differential equations of the problem are 
given in equations (70); the operational formulas for the case of 
voltage directly impressed on the terminals of the infinitely long 
line are 

/ = ^^ + 1 

Writing CRx^ — a and RGx^ — ^, G/C=^\ and assuming a ‘'unit 
e.m.f.” impressed on the cable, this becomes 

7=6““^^^, (175) 

/ = ^^Vp+Xe'“'^“2>+i3. (176) 

These equations are readily solved by means of the table and 
formulas given in a preceding chapter. 

But first let us attempt to solve the operational equation (175) 
for the voltage by Heaviside methods, guided by the solution of 
the operational equation 

(124) 

of the preceding chapter. Expand the exponential function in 
(175) in the usual power series; it is 

F=1 - .. (177) 


Now discard the integral terms and write 
F=l- 


Vap+jS. 


(178) 


We have now to interpret the expression -v/ap+iS- We have by 
ordinary algebra 


•\/ap+i3 = (l+^j Voep = \l + ~J V cep 


' /■^2p 2!(2pj ■^'3fV2pJ 


+ ■ 


v" ap. 


(179) 


Now identify Vp with l/Viri in accordance with the Heaviside 
rule, and 1/p with J'di. We get 


V«P+^ = ^^ 


/“I 

(\ty 

14 ] 


3! ^ 

h g 1 [M) 1 


(180) 




94 


ELECTRIC CIRCUIT THEORY 


Now in the terms of the expansion (178) identify p" with d^/df^ 
and substitute (180); we get 


y=i- 1+ 






^■^i! 3! ' 


.14(>^0:_, 

5! 


(181) 


This series is hopelessly complicated to either interpret or com¬ 
pute. It is, in fact, an excellent illustration of the grave dis¬ 
advantages under which many of Heaviside’s series solutions 
labor. We shall therefore attack the solution by aid of the 
theorems and formulas of Chapter IV. The simplicity of the 
solution which results is remarkable. 

The_operational formula for the voltage is 

7=e-V^P. (175) 

Now the operational formula for the voltage in the non-leaky 
cable is (see equation (164) ) 

7=g-\/a!). 


In order to distinguish between the two cases, let us denote the 
voltage in the latter case by F®; thus 

Fo = e-v'«P. (182) 

Now by theorem (Va) and equation (182) we have 


■yog->.t—^]P - g-Vaip+x), 

p+X 

= — p— V ap + j3 

p+x 


Now write (175) as 


Y =. _JL.. g- V<xi>+^ 
V p+x 




g-V«P+/5 

p+x® 


(183) 


(184) 


It follows at once by comparison with (183) and the rule that 
1/p is to be replaced by ydi, that 


F = 



(185) 


By a precisely similar procedure with the operational formula 
(176) for the current, we get 

l-t-xJ^diy»e-^‘ 



(186) 
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where Z® is the current in the non-leaky cable. Now by formulas 
(169) and (170) 


y®=-4- ® ,=dt, 

VTT Jo Wt 

(169) 

1 

o 

(170) 


which completes the formal solution of the problem. 

Formulas (185) and (186) are extremely interesting, first as 
showing the superiority of the definite integral to the series 
expansion—compare (185) with the series expansions (181)— 
and secondly as exhibiting clearly the effect of leakage on the 
propagated waves of current and voltage. We see that in both 
the current and voltage the effect of leakage is two-fold; first it 
attenuates the wave by the factor {\=G/C), and secondly it 
adds a component consisting of the progressive integral of the 
attenuated wave. This, it may be remarked, is the general 
effect of leakage in all types of transmission systems. Its effect 
is, therefore, easily computed and interpreted. 

Formulas (185) and (186) are very easy to compute with the 
aid of a planimeter or integraph; or, failing these devices, by 
numerical integration. However, for large values of t, the 
character of the waves is more clearly exhibited if we make use 
of the identity 

I di= I dt— I dt 


<•/ U v/ C 

whence 



F = (l-|-x| 

f Voe-'^dt 

%) i 

(187) 

\ 0 / 

and 



Z = (^l+xJ^ d<y«e-^‘-X 

^ 00 

j Pe-^dt. 

(188) 

The first two terms of these formulas 

are clearly the ultimate 


steady state values of the voltage and current waves, and can 
be determined by evaluating the infinite integrals. A far simpler 
and more direct way, however, is to make use of the fact that the 
ultimate steady values of V and I are gotten from the operational 
formulas by setting p = 0. That this statement is true is easily 
seen if we reflect that the steady d.c. voltage and current are 
gotten from the original differential equations of the problem by 
assuming a steady state and setting d/dt = 0, 
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From the operational formulas we get, therefore 

= = (189) 

d?jP e e - =^1 e- (190) 

Introducing these expressions into (187) and (188) respectively, 
we get 

7=e-=:V5ff_x I VU-^dt, (191) 

(192) 

The definite integrals can be expanded by partial integration; 
thus 

^05 f*CO 

-X 7‘>e-^«di= 7“de-w 

Jt Jt 

e-'^^-VOdt. 
dt 

Continuing this process we get 

F-e , .)f' (193) 

7. .(194) 

Using the values of V° and P, as given by (169) and (170), it is 
extremely easy to compute V and J, for large values of J, from 
(193) and (194). 

So far we have considered the current and voltage waves in 
response to a ^'unit e.m.f.,” impressed on the cable at x — 0. 
It is of interest and importance to examine the waves due to 
sinusoidal e.m.fs., suddenly impressed on the cable, particularly 
in view of proposals to employ alternating currents in cable 
telegraphy. 

We start with the fundamental formula 
^{t) =^J^/(^~"r)A(r)d:r 

= JW-r)ii'(T)dr, 

provided A(0) = 0, which is the case in the cable. 


. ■p'Og-Xf, 


■X 
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If /(i{) =sin o)t, we write 


Xsit) —sin a)tj cos o)t.h^(t)dt 

— cos sin o)t.h\t)dt. 
i. ii 

aje(0 = cos o^t f cos cot.h'(i)dt 

Jo 

+sin o)fJ sin oiUi'{t)dt. 


Similarly, if the impressed e.mi. is cos 


(194-a) 


(194-b) 


The investigation of the building-up of alternating currents and 
voltages, therefore, depends on the progressive integrals 


iS- 


^ cos o)t.h'(t)dtj 

n 

I sin o)th^(f)dt, 

Jo 


(194-c) 


For the case of the voltage wave on the non-inductive, non-leaky 
cable these integrals, by aid of equations (169), become, if we 
write oj' — aco/4, 


C = 



Q-l/r (jQg 

r-x/r 


rx/r 


dr, 

dr, 


(194-d) 


where, as before, 

For the current wave we have, by (170), 



-l/r 


cos co r 
r\/r 
sin to'r 


rx/r 


dr. 

dr. 


(194-e) 


For small values of r and w' these integrals can be numerically 
evaluated without great labor. Mechanical devices, such as the 
Coradi Harmonic Analyzer, are here of great assistance. In fact 
the Coradi Analyzer gives these progressive integrals auto¬ 
matically. It may be said, therefore, that a complete mathe¬ 
matical investigation of the building-up of alternating current 
and voltage waves in the non-inductive cable presents no serious 
difficulties, although the labor of computation is necessarily 



98 


ELECTRIC CIRCUIT THEORY 


considerable. One fact makes the complete investigation 
much less laborious than might be supposed. This is, if the 
foregoing integrals are calculated for a given value of oj', the 
results apply to all lengths of cable and all actual frequencies 
w/ 27 r, such that aco is a constant. Thus if we double the length 
of the cable and quarter the frequency, the integrals are 
unaffected. 

The solid curve of Fig. 7 shows the building-up of the cable 
voltage in response to an e.m.f. cos impressed at time ^ = 0. 
The frequency ci;/27r is so chosen that w' = aw/4 = 27r, and the curve 
is calculated from equations (194-b) and (194-e). The dotted 



curve shows the corresponding steady state voltage on the cable; 
that is, the voltage which would exist if the e.m.f. cos o)t had 
been applied at a long time preceding i = 0. We observe that, 
for this frequency, the building-up is effectually accomplished 
in about one cycle, and that the transient distortion is only 
appreciable during the first half-cycle. 

The case is very much different when a higher frequency is 
applied. Fig. 8 shows the building-up of the alternating current 
in the cable when an e.m.f. sin is applied at time f = 0. The 
frequency is so chosen that w'== aw/4 = lOr. The outstanding 
features of this curve are that the initial current surge is very 
large compared with the final steady state, and that the transient 
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distortion is relatively very large. It is evident that the fre¬ 
quency here shown could not be employed for signaling purposes. 
This curve has been computed from the steady state formulas, 
and equations (160) and (161) for the transient distortion. 

If the applied frequency co/2t is very high, the steady state 
becomes negligibly small, and the complete current is obtained 
to a good approximation by taking the leading terms of (160) 



Fig. 8.‘ 


-Non-inductive cable (G^—0), Building-up of alternating current. 

4 


Applied G.m.f. sin w ^ IOtt 


xmC 


and (161). Thus if the applied e.m.f. is sin and co is suflS.ciently 

large, the cable current is 

2 1 d 

by (160) and (170) while, if the impressed e.m.f. is cos o)t, it is 

2 /1 

^/txR Vr 

by (161) and (170). Here and r = 4i/a. 




CHAPTER VII 


THE PROPAGATION OF CURRENT AND VOLTAGE IN 
THE TRANSMISSION LINE 

We now take up the more important and difficult problem of 
investigating the propagation phenomena in the transmission line. 
The transmission line has distributed series resistance R and 
inductance L, and distributed shunt capacity C and leakage con¬ 
ductance G. It is the addition of the series inductance L which 
makes our problem more difficult and at the same time introduces 
the phenomena of true propagation with finite velocity, as 
distinguished from the diffusion phenomena of the cable problem. 
The cable theory serves very well for the problems of trans¬ 
oceanic telegraphy^ but is quite inadequate in the problems of 
telephonic transmission. 

If I denotes the current and V the voltage at point x on the line, 
the well known differential equations of the problem are: 

Replacing d/dt by p, these become 

(.cr+s)y—f/- 

From the second of these equations 

dx Cpl^Gdx^^ 

1 With the installation of the new submarine cable, continuously loaded 
with permalloy, this statement must be modified. In this cable, the induc¬ 
tance plays a very important part, and is responsible for the greatly increased 
speed of signaling obtainable. (See paper by Buckley, A. I. E. E., Aug., 
1925, “The Loaded Submarine Telegraph Cable.’O 
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(195) 


(196) 
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and substitution in the first gives 

{Lp+R){Cv+G)I = ^J. (197) 

Similarly if we eliminate J, we get 

{Lp+R){Cv^G)V=j^J. (198) 

If we assume a solution of the form 


where A and B are arbitrary constants, substitution shows that 
the solution satisfies the differential equation for V provided 


y^={Lp+R){Cp+Q). 

(199) 

From equation (196) it then follows that 


Lp-\-R 

(200) 




Now restricting attention to the infinitely long line extending 
along the positive x axis, with voltage Fo impressed directly on 
the line at = the reflected wave vanishes and we get 



F=Foe-r* 

/ = ^^±^Foe-T'*, 
y 

(201) 

Now let US write . 

y^=(Lp+R)iCp+G). 



y^=^2[(p+py~<T^] 

(202) 

where 

v = l/\/'LC, 

R , G 
" 2L'^'2C’ 

R G 
"" 2L 2C 


Then setting Vq = 
become 

1, the operational equations 

of the problem 



(203) 


{cA‘ 

\ vr w {p-rpr—<r-^ 

(204) 
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Now consider the operational equation, defining a new variable F: 

e~ j V(p+p)'“— 0-2 


F-p~ 


(205) 


V ip+py~<r^ 

It follows at once from our operational rules, and (203) and (204), 


that 


I = v{ C-\~G 


F = 





(206) 


(207) 


Our problem is thus reduced to evaluating the function F, from 
the operational equation (206). This equation can be solved by 
aid of the operational rules and formulas already given. The 
process is rather complicated, and there is less chance of error if 
we deal instead with the integral equation of the problem 




2-£r2 

Jo 




(208) 


V {v+py 

Now let us search through our table of definite integrals. We 
do not find this integral equation as it stands, but we do observe 
that formula (m) resembles it, and this resemblance suggests 
that formula (m) can be suitably transformed to give the solution 
of (208). We therefore start with the formula 

g—X'v/ps-fi - 

vp+i"i « 

This, regarded as an integral equation, defines a function which it 
zero for t<\ and has the value Jo('\/i^^X^) for ^>X; Jo being tht 
Bessel function of order zero. We now transform (m) as follows 

(1) Let \p — q and i/X — h. Substituting in (m) we get 

Now, in order to keep our original notation in p and replaci 
g to p and h by t; we get 

^00 .. 

(m.l 

(2) In (m.l) make the substitution p = q+fx and then in th 
final expression replace g by p; we get 

” . p-V(F+/i)*Tx® 

- (m.2 

/I Vip+py+\^ 


1 “ 
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X X 

(3) In (in.2) make the substitution p = -q and and 


ultimately replace g by p and ia by if; we get 




V(p+/Xi)2 + Xi2 


(m.3) 


V V' 

whereXi = -X and (They are, of course, as yet, arbitrary 

X X 

parameters, except that they are restricted to positive values.) 

(4) Now if we compare (m. 3) with the integral equation (208) 
for Pj we see that they are identical provided we set 


Xi = 2<r =(r'\/ —1. 


Introducing these relations, we have 



.6“^^ Jo (cr-\/if2 — X^fv^)dt == 


-'“V'(p+p) 2 -«r 2 
e ^ 

V (p+py—o-^ 


(m.4) 


Here Jo denotes the Bessel function of imaginary argument: thus 

It follows from (m.4) and the integral equation (208) that 

F(^) =0 for Kxjv^ (209) 


Having now solved for 7 ’=jP( 0, the current and voltage are gotten 
from equations (206) and (207). Thus 

7 = 0 for i<x/Vj (210) 


= JrF(t)+vG f FiOdt fort^x/v. 
yi j-j Jx/v 


The corresponding voltage formula is 
7 = 0 for t<x/Vi 


== Q-PX/VJ^.. 


crx r e 
V J%/v 






’X^/v^ 


dTiovt>x/v, (211) 


Here hia's/r^—x^/v^) is the Bessel function of order 1: thus 
—iJi {iz) ^Ii{z). The function is entirely real. The derivation 
of formula (211) is a little troublesome, owing to the discontinuous 
character of the function F) the detailed steps are given in an 
appendix to this chapter. 

The preceding solution depends for its outstanding directness 
and simplicity on the recognition of the infinite integral identity 
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(m), into which the integral equation of the problem can be trans 
formed. When such identities are known their value in connectioi 
with the solution of operational equations requires no emphasis 
On the other hand, we cannot always expect to find such an iden 
tity in the case of every operational equation; and, particularb 
in such an important case as the transmission equation it wouh 
be unfortunate to have no alternative mode of solution. For 
tunately a quite direct series expansion solution is obtainable fron 
the operational equation, and this will now be derived. As i 
matter of convenience we shall restrict the derivation to th 
voltage formula 


As a further matter of mere convenience we shall assume tha 
(7= 0, so that a=p and (203) becomes 

V = e-rV'p 2 + 2 pp (203~a 

where r — x/v. 

The method holds equally well for the current equation (204) an 
for the general case o-T^p. 

Write (203»a) as 

y — (1+2p/p) 1/2 


and expand the exponential factor (l+2p/p)^^^ by the binomii 
theorem; thus 


(l+-2p/p)i/^ = l + ™+ai2 

r 


d-aai " I +• 

\V 


so that 


F=. exp 


a2 rp ^ rp a.\ rp ^ 


Now the operational equation 

«;=exD(^- 
\ V 

can be expanded in inverse powers of p; thus 

P P^ 

the power series solution of which is 

^ ^ 1! ^ 2! ^ 3! ^ 
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It follows at once from the preceding and Theorem VII that 
F(<)=0forf<T 

for # >r. 

If the coefficients . . . are evaluated, a simple matter of 

elementary algebra, the foregoing expansion in the retarded time 
t—T will be found to agree with the solution (211) when cr is put 
equal to p. 

We shall now discuss the outstanding features of the propaga¬ 
tion phenomena in the light of equations (210) and (211) for the 
current and voltage. We observe, first, that we have a true 
finite velocity of propagation ?; = 1/'\/LC. No matter what the 
form of impressed e.m.f. at the beginning of the line (a; = 0), 
its effect does not reach the point x of the line until a time i — xjv 
has elapsed. Consequently v — xli\^ the velocity with which the 
wave is propagated. This is a strict consequence of the dis¬ 
tributed inductance and capacity of the line and depends only on 
them, since y = l/^/LC. It will be recalled that in the case of the 
cable, where the inductance is ignored, no finite velocity of propa¬ 
gation exists. 

The question of velocity of propagation of the wave has been 
the subject of considerable confusion and misinterpretation when 
dealing with the steady state phenomena. It seems worth while 
to briefly touch on this in passing. 

As has been pointed out in preceding chapters, the symbolic or 
complex steady state formula is gotten from the operational 
equation by replacing the symbol y by ?co where i — and 

a)/27r is the frequency. If this is done in the operational equation 
(203) for the voltage, the symbolic formula is 


If the expression '\/(ta)+p)^~or2 is separated into its real and 
imaginary parts we get an expression of the form 




where 




2 _ p2(^2 _ p2) 2 q. 4^‘ip?. 


2co" 


and 


a=^p/^v. 
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Now if we keep the expression t—p- constant, that is, if we move 

along the line with velocity dxfdt — v(p, the phase of the wave will 
remain constant. This is interpreted often as meaning that the 
velocity of propagation of the wave is v/l3. Now since p is 
greater than unity and only approaches unity as the frequency 
becomes indefinitely great, the inference is frequently made that 
the velocity of propagation depends upon and increases to a 
limiting value y, with the frequency. This velocity, however, is 
not the true velocity of propagation, which is always v, but is the 
velocity of phase propagation in the steady, state. This distinction 
is quite important and failure to bear it in mind has led to serious 
mistakes. 

Returning to equation (210) and (211) we see that after a time 
interval has elapsed since the unit e.m.f. was impressed on 

the cable, the voltage at point x suddenly jumps from zero to the 
value-while the current correspondingly jumps to the value 



The exponential factor px/v is 



\2\L^2'VC/ 


= ax 


which will be recognized as the steady state attenuation factor for 
high frequencies. Similarly V'C/L is the steady state admittance 
of the line for high frequencies. The sudden jumps in the current 
and voltage at iim.Qt-x/v are called the heads of the current and 
voltage waves. If, instead of a unit e.m.f., a voltage f{t) is 
impressed on the line at time t = 0, the corresponding heads of the 
waves are /(O)^"®® and 's/C/L f{0)e'~^^ for voltage and current 
respectively. These expressions follow at once from the integral 
formula 


=fmit)+j‘f(t-TMr)dr. 


The tails of the waves, that is, the parts of the waves subse¬ 
quent to the timet —x/vj are more complicated and will depend on 
the distance x along the line and on the line parameters p and a. 
The two simplest cases are the non-dissipative line, and the 
distortionless line. 
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The ideal non-dissipative line, quite unrealizable in practice^ is 
one in which both R and G are zero. In this case p=o- = 0, and 
formulas (210) and (211) become 

7 = 0 for t<x/vy 



for t>x/Vy 


F==0 for t<x/Vj 
— 1 for t>x/v. 


Both current and voltage jump, at time t = x/Vy to their steady 
values. If an e.m.f. /(i) is impressed on the line at time t — 0, the 
corresponding current and voltage waves are 
7 = 0 for t<x/Vy 

— — for t>x/v, 

V = 0 for t<x/Vy 
=f(t — x/v) for t>x/v. 

Consequently the ideal non-dissipative line transmits the waves 
with finite velocity v, without attenuation or distortion. Such a 
line is, of course, the ideal transmission system. 

The non-dissipative line is, of course, purely theoretical and 
unrealizable in practice; the distortionless line is, however, approx¬ 
imately realizable, and as the name implies, transmits without 
distortion of wave form. The distortionless line is one in which 
the line constants are so related that 

2L -20 


If this condition is satisfied, formulas (210) and (211) become 
7 = 0 loxi<x/Vy 



for t>xfvy 


F = 0 for i <x/Vy 

fQj. i^x/v. 


Furthermore, if the impressed e.m.f. is /(i), the corresponding 
current and voltage waves are: 

7 = 0 for 



fii—x/v) for t>x/vj 


F = 0 for 

— e'^^^f(t'--x/v) iQTi>x/v, 
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The distortionless line, therefore, transmits-the waves without dis¬ 
tortion of wave form, but attenuates the waves by the factor 
Such a line is an ideal transmission system as regards preservation 
of wave form, but introduces serious attenuation losses. For 
example, if a line has normally negligible leakage, and leakage is 
introduced to secure the condition i2/L=G^/C, the line is thereby 
rendered distortionless but the attenuation is doubled. 

One of Heaviside^s most important contributions to wire 
transmission theory was to point out the properties of the dis¬ 
tortionless line, its approximately realizable character, and to 
base on it a correct theory of telephonic transmission. 

The character of the wave propagation when the parameters 
p and a are not restricted to special values, can only be roughly 
inferred from inspection of the formulas, and then only when 
the properties of the Bessel function Iq and Ii have been studied. 
Fortunately these functions have been computed and tabulated 
for small values of the argument, and have simple asymptotic 
expansions for large values. It is therefore a simple matter to 
compute and graph a representative set of curves which show 
the current and voltage waves for various values of p, <t and x. 
For this purpose it is convenient to introduce a change 
of variables and write: 

r^^vt 

a=pfv 

b—crfv 


whence the formulas for current and voltage become: 




V = 






(210-a) 


(211-a) 


Figs* (9) to (18) give a representative set of curves illus¬ 
trating the form of the propagated current and voltage waves 
for different lengths of line, and different values of the line 
parameters a and b, or p and cr. 

The curves of Figs. (9) and (10) show the current entering the 
line in response to a unit e.m.f. applied at time i=0. The line 
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is assumed to be non-leaky (& = 0) and is computed for two 
different values of the parameter a. We see that the current 
instantly jumps to the value \/CIL and then begins to die 



away, the rate at which it dies away depending on and increasing 


with the parameter a — 




If we now consider a point x out on the line, the current is 
zero until t = Xj at which time it jumps to the value \^CIL 
It then begins to die away provided x and a are such that ax<2. 
If, however, we are considering a point at which ax>2j the 
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current begins to rise instead of fall after the initial jump^ and 
may attain a maximum value very large compared with the head 




Fig. 12. —Propagated curront in lino; a?200; »0.02645; G =0. 

Multiply ordinates by 


before it starts to die away. This is shown in the curves of 
Figs. (11), (12) and (13), also computed for the non-leaky line 
(6 = 0). From these curves we see that, as the length of the line 
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and the parameter a increase, the relative magnitude of the tail, 
as compared with the head of the wave, increases. Finally 
when the line becomes very long, the head of the wave becomes 




Fia, 14.—Propagated current in cable. 
Multiply ordinates by 2/Rx, 


negligibly small, and the wave, except in the neighborhood of 
its head, becomes very close to that of the corresponding non- 
inductive cable. This is shown in curves (13) and (14), for the 
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line and the corresponding cable, which are plotted to the same 
time scale and ordinate scale to facilitate comparison. Curve 
(15) shows the effect of leakage in eliminating the tail. This 
line is not quite distortionless but nearly so. 
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Fia. 15.—Propagated current in line; x =200. 


Multiply ordinates by v/ 



An interesting feature of both current and voltage waves is 
that when a sufficient time has elapsed after the arrival of 
the head of the wave, the waves become closer and closer to the 
waves of the corresponding non-inductive cable; that is, to the 
cable having the same B, C and (?. Consequently the inductance 
plays no part in the subsidence of the waves to their final values. 
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Curves (16), (17) and (18) illustrate the voltage wave for 
several conditions. After the arrival of the head, the wave 
slowly builds up to its final value. Curve (18) represents the 



case where the line is very nearly distortionless, showing how 
completely the distorting tail of the wave is eliminated. 

So far we have confined attention to the current and voltage 
waves in response to a unit e.m.f. applied at time ^==0 to the 



line terminals. Of much greater technical importance is the 
question of the waves in response to a sinusoidal e.m.f. suddenly 
applied to the line terminals. In order to investigate this 
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important problem it is convenient to divide the expressions for 
the current and voltage waves as given by equations (210-a) 
and (211-a) into two components. We write for T>Xy 

/ = (210-b) 

7=e-«+Tr(i), (211-b) 

where, by definition, J{t) and W{t) are the differences between 
the total waves and their heads. The advantage of analyzing 
the waves into these components is that the distortion of the 
waves is due to J{t) and W(t) respectively, while the first compo¬ 
nents of (210'-b) and (211-b) introduce merely a delay. Thus, 
if the e.m.f. impressed at time ^ = 0 is /(O, the corresponding 
waves for f>x/v or r>Xj are 



(212) 



+ 

I 

(213) 


where /'(«) = f/CO and W(t). 

The integrals of (212) and (213) can be computed and analyzed 
in precisely the same way as discussed in connection with the 
non-inductive cable problem, and are of very much the same 
character as the corresponding waves of the cable. In the 
total waves, however, as given by (212) and (213), a very essen¬ 
tial difference is introduced by the presence of the first terms, 
which represent undistorted waves propagated with velocity v. 
Thus, if the impressed e.m.f. is sin cot, (212) and (213) become 

I = sin c>o{t-'x/v) 

+ r sin co(^ —for (214) 

Jx/v 

V sin co{t — x/v) 

+ r ^mcc(t--ti)W'(tx)dti,fort>x/v, (215) 

Jx/v 

Now the first terms of (214) and (215) are simply the usual 
steady state expressions for the current and voltage waves when 
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the frequency is sufficiently high to make the steady state atten¬ 
uation constant equal to a and the phase velocity equal to v. 
Furthermore the integral terms become smaller and smaller 
as the applied frequency a)/2x is increased. It follows, there¬ 
fore, that for high frequencies the waves assume substantially 
their final steady value at time f==xfv, and that the tails of the 
waves, or the transient distortion, becomes negligible. This is 
a consequence entirely of the presence of inductance in the line, 
and shows its extreme importance in the propagation of alter¬ 
nating waves and the reduction of transient distortion. 

It should be pointed out, however, that if the line is very long 
and the attenuation is very high, the integral terms of (214) and 
(215) are not negligible unless the applied frequency is corre¬ 
spondingly very high. For example, on a long submarine cable, 
the a.c. attenuation is so large that the first terms of (214) and 
(215) are very small, and J(t) is very large compared with 
's/Cjh Consequently here there is very serious transient 

distortion and alternating currents are therefore not adapted 
for submarine telegraph signalling. 

This discussion may possibly be made a little clearer, without 
detailed analysis, if we recall the discussion of alternating current 
propagation in the non-inductive cable of the preceding chapter. 
From that analysis it follows that, when the applied frequency 
co/27r is sufficiently high, the integral term of (214) becomes 
approximately 

and the complete current wave is 

sin (216) 

and similarly the voltage wave is 

e"”®® sin o){t—x/v)+^W\t). (217) 

Now if the total attenuation ax is large the last terms of (216) 
and' (217), before they ultimately die away, may become very 
large compared with the first terms, which represent the ultimate 
steady state. 

APPENDIX TO CHAPTER VII 

Derivation of Formula (211).—The only troublesome question 
involved in deriving (211) from (207) and (209) is that we have to 
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differentiate with respect to in accordance with (207), the 
discontinuous function F{t). To accomplish this we write (209) 
in the form 

F(t) = <j>(t ~xfv)e-'‘‘ (209-a) 

where <l>(t) is defined as a function which is zero for tKxJv and 
unity for i>x/v. Clearly this is equivalent to (209) and permits 
us to deal with F{t) as a continuous function. Now, in accord¬ 
ance with (207), perform the operation of differentiation upon 
(209-'a): we get 

x/ I o(o-\/ F' — x'^/v^) . 

The first expression follows from the fact that 

We observe also that "^4>{t'-x/v)—0 except at i — when it is 
infinite. We also obseiwe that, for t>xjvj 
J ^^(t>it—x/v)dt^l 

and that the whole contribution to the integral occurs at i^x/v. 
With these points clearly in mind, the expression 



reduces to (211) without difficulty. 



CHAPTER Vlir 


PROPAGATION OF CURRENT AND VOLTAGE IN 
ARTIFICIAL LINES AND WAVE FILTERS 

The artificial line here considered is a periodic structure, 
composed of a series of sections connected in tandem, each section 
consisting of a lumped impedance in series with the line, and 
a lumped impedance 02 in shunt across the line. In the artifi.cial 
line which we shall consider it will be assumed that the voltage 
is applied at the middle of the initial or zeroth section, as shown 
in Fig. 19. This termination is chosen because of its practical 
importance, and because also of the fact that the mathematical 
analysis is simplified thereby. Furthermore any other termi- 



Fig. 19. 


nation can be regarded and dealt with as an additional terminal 
impedance, so there is no essential loss of generality involved. 

A study of the properties of the artificial line is of practical 
importance for several reasons; 

1. The artificial line is often used as a model of an actual 
transmission line and it is therefore of importance to determine 
theoretically the degree of correspondence between the two. 

2. The solution for the corresponding transmission line with 
continuously distributed constants is derivable from the solution 
for the artificial line by keeping the total inductance, resistance, 
capacity and leakage constant or finite, and letting the number 
of sections approach infinity. 

3. The artificial line is very closely related, in its properties 
and performance, to the periodically loaded line, and its solution 
is, to a first approximation, a working solution for the loaded 
line. 


117 





118 


ELECTRIC CIRCUIT THEORY 


4. The structure is of great importance in its own right, and 
when the impedance elements are properly chosen, constitutes 
a ^^wave filter/' 

We shall now derive the operational and symbolic equations 
which formulate the propagation phenomena in the artificial 
line. Let In denote the mesh current in the section of the 
line; In-^i the mesh current in the section etc. Now 

write down the expression for the voltage drop in the section; 
in accordance with Kirchhoff^s law we get: 

=0 (218) 

where, of course, the impedances have the usual significance. 

Now this is a difference equation, as distinguished from a 
differential equation, but the method of solution is essentially 
the same. We assume a solution of the form 

(219) 

where A, B and T are independent of n, and substitute in (218). 
After some simple rearrangements we get 

{ (Z1+2Z2) —2z2 cosh T]'[Ae'^^+Be^^} =0 (220) 

Equation (218) is clearly satisfied by the assumed form of solu¬ 
tion, and furthermore leaves the constants A and B arbitrary 
and at our disposal to satisfy any boundary conditions, provided 
r is so chosen that 


cosh r = 


2Z2 

-l+2p 


( 221 ) 


where p^ZxjAz^, 

Now by reference to equation (219) it is easily seen that V is the 
propagation constant of the artificial line, precisely analogous to 
the propagation constant 7 of the smooth line. In terms of the 
impedances Zi and the propagation constant of the artificial 
line is determined by (221). This equation may either be 
regarded as an operational equation or a symbolic equation, 
depending on whether the impedances are expressed in terms of 
the operator p or in terms of fco, where co is 27r times the frequency. 

Now suppose in (221) we write e^^x) the equation becomes 


a;+l/aj = 2(l+2p) 
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and solving for x we get 

x = e^ = {l +2p)+V(l+2p)2-l 

_ ^ _ _ ( 222 ) 

= (Vl+io+\/p)^ = ("S/l+P — Vp)"^ 


which is an explicit formula for r. 

Now return to equation (219) and let us assume that the line 
is either infinitely long, or, what amounts to the same thing, 
that it is closed by an impedance which suppresses the reflected 
wave. We assume also that a voltage Fo is impressed at mid¬ 
series position of the zero"* section (n = 0). Equation (219) 
becomes 

and the currents in the zero"* and 1“‘ sections are 
Jo=4, Ix=Ae-^. 


Now, by direct application of Kirchhoff’s law to the zero"* 
section, we have 

Fu = (j 2 i +22) Jo — zJi, 

whence 

A{^Zi-{-Z2(l—e^)} = Vo. (223) 


But 


Jo A 

In J^e , 


Vo, 


where K is the characteristic impedance of the artificial line (at 
mid-series position). Hence by (223) and (222) 


1 1 1 . 

= .1 1 ^ __ 1 _ 1 
222\/p-)-p^ "s/zi 02-v/l-f-P 


(224) 


By aid of the preceding the direct current wave can be written as 


■ = .Vj_ [Vi+p-Vp]^" 
VzA Vi+p 


(225) 


This formula is not so physically suggestive as its equivalent 
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but is useful when we come to the solution of the operational 
equation. 

Before proceeding with the operational equation, and the 
investigation of transient phenomena in artificial lines, it will be 
of interest to deduce from the foregoing the unique and remark¬ 
able properties of wave filters in the steady state. For this 
purpose we return to equation ( 221 ) 

cosh r = l+ 2 p. 

Now suppose that the series impedance Zi is an inductance L and 
the shunt impedance a capacity C, so that, symbolically, 

1 cc^LG 

Zi=tC^L, ^2^^^ P- - - 4 “ ’ 

and 

cosh r = 1 —^(jo^LC. (226) 

Now let us write r==i(9, where i — ^ ~ 1 ; the preceding equation 
becomes 

cos (9 = l-ia)2LC (227) 

and the ratio of currents in adjacent sections is Consequently 
if Q is a real quantity the ratio of the absolute values of the currents 
in adjacent sections is unity, and the current is propagated without 
attenuation. 

Inspection of equation (227) shows that 9 is real provided the 
right hand side lies between +1 and — 1 ; or that co lies between 
0 and 2I\^LC. Consequently this type of artificial line trans¬ 
mits, in the steady state, sinusoidal currents of all frequencies 
from zero to If'jr's/LC without attenuation. It is known as the 
low pass filter. 

If we invert the structure, that is, make the series impedance Zi 
a capacity C and the shunt impedance Z 2 an inductance L, so 
that 

we get, corresponding to (226) and (227), 

cosh r = 1 - (228) 

cos 0 = 1(228-a) 

This type of artificial line transmits without attenuation currents 
of all frequencies for which the right hand side of (228-a) lies 
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between +1 and —1; that is, all frequencies from infinity to a 
lower limiting frequency I/ItVI/C, while it attenuates all 
frequencies below this range. It is known, on this account, as the 
high pass filter. 

It is possible by using more complicated impedances to design 
filters which transmit a series of bands of frequencies. We 
cannot, however, go into the complicated theory of wave-filters 
here, which has been covered in a series of important papers.^ 
One point should be noted, however: transmission without 
attenuation implies that the impedance elements are non- 
dissipative. Actually, of course, all the elements introduce 
some loss, so that in practice the filter attenuates all frequencies. 
Careful design, however, keeps the attenuation very low in the 
transmission bands. 

We shall now derive the indicial admittance formulas for 
some representative types of artificial lines and wave-filters 
from the operational formula 

= .[Vi+p+Vp]"'"- (229) 

V(l + p)^1^2 

This equation follows directly from (225) on putting 70 = 1. 

We start with the so-called low pass filter on account of its 
simplicity and also its great importance in technical applications. 
This type of filter consists of series inductance L and shunt 
capacity C. The general case which includes series resistance 
R and shunt leakage G has been worked out (see Transient 
Oscillations, Trans. A.I.E.E., 1919). The solution is, however, 
extremely complicated and will not be dealt with here. We shall, 
instead, consider the important and illuminating case where the 
series and shunt losses are so related as to make the circuit 
quasi-distortionless. We therefore take, operationally, 

^,^pL+iJ==L(p+X) (230) 

1/^2 “PO -!-(?=C(p 4“ k) 

where \ — R/L—G/C. 

We then have 

z^z,^LlC, 

^iA2-LC(p+X)2, (231) 

p = _(p + X)2. 


^ See Bibliography. 
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Now by reference to formula (229) we see that A„ is a function of 
(p+X); thus 

Z„(p+X) V 

Now write 

AO _L_. 

" Znip) 


It follows at once from reference to theorem V that 



(232) 


so that the problem is reduced to the solution of the operational 
equation for A^. Writing Uc = ‘2f^/LC, we have 


k-4 


i Vi+(pAo)^ 


Vl + (p/tJc) ^+p/Wc 


-2?i 


(233) 


Ic CO. 



COc 


Now refer to formula (n) of the table of 
■\/L/C=1s, we see by Theorem II that 



integrals; writing 
(234) 


where Jinir) is the Bessel function of order 2n and argument r. 
We note also that this is the indicial admittance of the non- 
dissipative low pass wave-41ter; that is, the current in the n”' 
section in response to a unit e.m.f. applied to the initial section 
(n=0). From (232) and (234) it follows at once that 



Integrating the second member by parts and noting that 
(0) =0, this reduces to 

= e o>IJ 2 „(r)dT (235) 

which is the indicial admittance formula for the quasi-distortion- 
less low pass filter, or artificial line. 
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Before discussing these formulas, it is of interest to derive the 
formula for Al by power series expansion. Formula (233) can be 
written 

10_l/^cY" + ^ 1 1 

k\p / •\/i + (&5c/p)*[lH-'\/l + (Ws/p)^]^’* 

This can be expanded in a series in inverse powers of p; thus 
4“ — ^ f { 2«+i_2n+2/w<,Y"+3 

, (2n+3) (2w+4)/co,\="+' 1 

" 2^2! (p) ■■■■]• 

Replacing l/p” by i^/n! in accordance with the Heaviside Rule we 
get 

2f 1 /o3,ty’^+^ 2n+2 

^”~k\ (2n+i)'\2 ) l!(2n+3)!\2 / 

(2n+3) (2n+4)/w„fY'*+'‘ 

2!(2n+5)! 2 j 

This can be recognized as the power series expansion of (234). 

The artificial cable is also of interest and practical importance. 
In this structure the series impedance is a resistance B and the 
shunt impedance is a capacity C, so that 

Zi=B, l/Zi=pC, 

z,z^=BlpC,Zilz, = pRC, (236) 

p=pRC/4:. 


(235-a) 


Now let us return to formula (229), and expand in inverse powers 
of p: we get 


An- 


pziZi 


2n+2 1 , 

2*1! p”+i”^ 


(271+3) (2w+^ 

2 ^ 2 ! p"+^ 


( 237 ) 


Now since V* pziZi=\> we have 


JL\( Y 277+2/ 2 Y^^ 

^^~2’'B\\RCp) ”2.1! \RCp) 

, (277+3) (271+4)/ 2 V+2_ I 

2*21 \RCp) I’ 
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Replacing 1/p” by V'/n\ we get finally 

2 \l/2t\ (2n+2) /2 iY+i 

2”;? 2.l\{n+l)\\Rc) 

(2n+3)(2n+4)/2i! Y+“ 1 

22.2!(n+2)! \Rc) ""j- 


(238 


For large values of n and t this series is difficult to conapute c 
interpret. It can, however, be recognized as the series expansio: 
of the function 

2_21 / 2t\ 

4„=^e (239 

where In{2t/R'C) is the Bessel function I of order n and argumen 
(2t/BC). This solution, it may be remarked, can be derives 
directly by a modification of the integral formula (n). 

It is beyond the scope of this paper to consider other types o 
artificial lines and wave“filters;for a fairly extensive discussion th 
reader is referred to '‘Transient Osillations in Electric Wave 
Filters/^ B. S. T. J., July, 1923. The low pass wave-filter, how 
ever, both in its own right and on account of its close relation t< 
the,periodically loaded line, deserves further discussion. 

For the non-dissipative low pass wave-filter, we have 

= (234; 

while for the quasi-distortionless low pass wave-filter 

(235; 

where m = \/o)e ^R/Lo)c = Rl2vL. 

Computation and analysis of these formulas involve an ele 
mentary knowledge of Bessel functions. The properties neces 
sary for our purposes are briefly discussed in an appendix to thif 
chapter. 

The indicial admittances for the non-dissipative low pass filter 
that is, the current in response to a steady unit e.mi. applied ai 
time ^=0, are shown in the curves of Figs. 20, 21 and 22, for th( 
initial or zeroth, the 3rd and the 5th sections, respectively. These 
curves together with the exact and approximate formulas giver 
above are sufficient to give a reasonably comprehensive idea oj 
the general character of these oscillations and their dependence or 
the number of sections and the constants of the filter. 
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It will be observed that the current is small until a time approx¬ 
imately equal to 2nfcoc — n'\/LiC 2 has elapsed after the voltage is 
applied. Consequently the low pass filter behaves as though enr¬ 
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Fig. 20.—^Low pass wave filter. Indicial admittance of initial section (n=0). 
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Fig. 21. —Low pass wave filter. Indicial admittance of third section (w =2). 
Multiply ordinates by 's/'CJL. 
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0 

Fig. 22.—^Low pass wave filter. Indicial admittance of fifth section (n —4). 
Multiply ordinates by C/L. 



rents were transmitted with a finite velocity of propagation 
o)cI2 = II\/LC sections per second. This velocity is, however, 
only apparent or virtual since in every section the currents are 
actually finite for all values of time >0. 
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After time t = n\/LC has elapsed the current oscillates about 
the value l/k with increasing frequency and diminishing ampli¬ 
tude. The amplitude of these oscillations is approximately 

1/fc Ij2 

Uet)^\ ITU ct 

and their instantaneous frequency (measured by intervals between 
zeros) 

~-\/l-(2n/coc0^ 

The oscillations are therefore ultimately of cut-off or critical 
frequency w<,/2;r in all sections, but this frequency is approached 
more and more slowly as the number of filter sections is increased. 



Fia. 23.—^Low pass wave filter. Indicial admittance of 100th section (n =99), 
Multiply ordinates by 's/ C/L. 

Figs. 23, 24, 25, give the indicial admittance in the 100th, 500th 
and 1000th section of the low-pass wave-filter. The filter itself 
seldom embodies more than 5 sections. The case of a large 
number of sections is of interest, however, because it represents 
a first approximation to the periodically loaded line. While the 
non-dissipative line is ideal and unrealizable, its study is of practi¬ 
cal importance because in this type of line the effect of the dis¬ 
continuous character of the loading of the periodically loaded 
line is isolated and exhibited in the clearest possible manner. 

The dotted curves represent the current in the corresponding 
smooth line, For the smooth line, the current, as we have seen, is 
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discontinuous, being identically zero for a time vt=^n and having 
an instantaneous jump to its final value ^s/CjL at vt^n. The 
current in the artificial or periodically loaded line differs from that 



Fig. 24.—Low pass wave filter. Indicial admittance of 600th section (n =499) 
Multiply ordinates by ”>/ C/L. 



Fig. 25.—Low pass wave filter. Indicial admittance of 1000th section (n =999). 
Multiply ordinates by *\/ C/L, 

in the corresponding smooth line in three important respects: (1) 
the absence of the abrupt discontinuous wave front, (2) the 
presence of superposed oscillations; and (3) the absence of a true 
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finite veloctiy of propagation. It will be observed, however, 
that the current in any section is negligibly small or even sensibly 
zero until vt = nj so that the current is propagated with a virtual 
velocity lI's/LC per section. The presence of a well marked wave 
front is also evident although this is not abrupt, as in the smooth 
line. The effective slope of the wave front becomes smaller 
as the current wave travels out on the line, decreasing noticeably 
as the number of sections is increased. When the number of 
sections becomes large, however, the decrease in the slope is not 
rapid, being in the 500th section about 60 per cent, of that in the 
100 th section. 

The superposed oscillations are of interest. These are initially 
of a frequency depending upon and decreasing with the number of 
sections, n, but in all sections ultimately attaining the frequency 

_ 1 _ 

tta/LC ^ 

which is the critical or cut-off frequency of the line, above which 
steady state currents are attenuated during transmission and 
below which they are unattenuated. When vt is large compared 
with n the amplitude of these oscillations becomes ^/T/irv't so 
that they ultimately die away and the current approaches the 
value '\/CIL for all sections. The current in the loaded line is 
thus asymptotic to the current in the corresponding smooth line 
and oscillates about it with diminishing amplitude and increasing 
frequency. 

Since the abscissas of these curves represent values of 2vt = 
2tl\/LC 3 and the ordinates are to be multiplied by ■\/C/L to trans¬ 
late into actual values, the curves are of universal application for 
all values of the constants L and C. 

The investigation of the building up of alternating currents in 
wave filters and loaded lines is very important. It depends for 
the non-dissipative case on the properties of the definite integrals 

£ 

sin ivr Jn(r)(lT, 

J cos WT Jn{r)dr, 

where w — (^fo)o and co = 27r times the applied frequency. The 
mathematical discussion is, however, quite complicated and will 
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not be entered into here. The reader, who wishes to follow this 
further, is referred to Transient Oscillations, Trans. A. I. E. E., 
1919 and Transient Oscillations in Electric Wave Filters, 
B. S. T. J., July, 1923. 

APPENDIX TO CHAPTER VIII 

Note on Bessel Functions.—The Bessel Functions of the first 
kind, Jn{x) and are defined, when n is zero or a positive 

integer, by the absolutely convergent series 

I ^ ~2(2n+2)■^^(2n+2)^+4) 

a;' 4- 1, 

2.4.6(2ra+2)(2n+4)(2ra+6)‘^ j 

[ /y»2 /v»4 

J _ .,_5_1 I _I___ 

2’‘.n! I ^2(2«+2)^2.4(2??,+2)(2n+4) 

._?!_J_ \ 

^2.4.6.(2«+2)(2n+4)(2n+6)^ |‘ 

In the following discussion of the properties of these functions it 
will be assumed that the argument a; is a pure real quantity. 

For large values of the argument (a: large compared with n), the 
behavior of the functions is shown by the asymptotic expansions: 
, . . e® f- 4n^—1 . (4n^—l)(4n^—9) 

1 I'KSa:) ^ 2 !(8a;)^ 

(4w2-l)(4?^=-9)(4?^2-25) 

-“ 'sTCSa;)®' ^ 

r / ^ I a / 2n+l r\ ■ r 2n+l . I 

«f^«(a;) = '^— Pn cos (x— ^-“7r)-Q„sm (a:- -^—t) |» 

where 

p _ (4n2-l)(4n2-9) 

“2!(8a:)2 

(4?i2 -1) - 9) (4n2 - 25) (4n^ - 49) 

^ ” 4l '(8^ , '-•••> 

_ 4„2_1 (4n2-l)(4n2-9)(4ra2_26) 

‘^"“'8a:'~ 3!(8a:)3' “ "'■■■■■ 

We thus see that I„ increases indefinitely and behaves ulti¬ 
mately as 

_e®_ 

'\/2xTr 
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The function Jn{x), however, is oscillatory and ultimately 
behaves as 


For all orders of n 



2n+l 

4 




The properties of Jn{x) may be described qualitatively as 
follows:— 

When the argument is less than the order (0^a;<?^) the func¬ 
tion is very small and positive, and is initially zero (except when 
7 ^ = 0). In the neighborhood of x — n, the function begins to 
build up and reaches a maximum a little beyond the point x^n. 
Thereafter the function oscillates with increasing frequency and 
diminishirig amplitude, and ultimately behaves as 

lY ( 2n+l \ 

IT-"/ 

When n = 0, the initial value is unity, but the subsequent behavior 
of the function is as described above. 

A more precise description of the function is gotten from the 
following approximate formulas.’- 


where 




Jn(oc) — Bn(x) cos for x>n 


t -.o I o 


y 2 (1 —m' 


1 _ 


r, / \ r /i , m . ,[m\ 1 

n{x) X yj ^2 4a;‘l (1 — 


2n+l 


(^) — ^ 


•i- 


^ 1 

x^ 2 x"^ {i—m^fx^y^ 


1 Th«se formulas are due to Mr. Zobel (See ‘^Transient Oscillations in 
Electric Wave-filters” Bell System Technical Journal, July, 1923.) Some¬ 
what similar,formulas had previously been developed by Graf and Gubler 
(Einleitung in die Theorie der Besselschen Funktionen), and by Nicholson 
(Phil. Mag., 1910, p. 249), 
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and 

This approximate formula is valid only where x>7ij its accuracy 
increasing with x and with n. For all orders of n it is quite 
accurate beyond the first zero of the function. 

The ^instantaneous frequencyof oscillation is approximately 

1 


2t 




-iW'- 


, 3^ 1 

" 0:2 ^2'^ (i - m^/x^) 2 ‘ 


By this it is meant that at any point x{x>n) the interval between 
successive zeros is approximately 7r/0'(o:). Otherwise stated, in 
the neighborhood of any point Xj the function behaves like a 
sinusoid of amplitude Bn{x) and frequency (jo/2t where w == 12'n(x). 

The following approximate formulas, while not sufficiently 
precise for the purpose of accurate computation except for quite 
large values of Xj clearly ^exhibit the character of the functions 
for values of the argument^a;>n, and of the order n>2. 

cos {qn^-Bn), 




/: 


JU^) = -g»hnyl- 


2 

TX 

2 

TTX 


sin 0„), 

sin (qni^-Bn), 


where 




y/i 


l + — 
4x^ 


Sn- 

On- 


= •\/l—nVa;^, 
2^1 
4 


TT—n sin-^ (n/x). 



CHAPTER IX 


THE FINITE LINE WITH TERMINAL IMPEDANCES 


So far in our discussions of wave propagation in lines and wave- 
filters, we have confined attention to the case where the impressed 
voltage is applied directly to the infinitely long line. We have 
found that, by virtue of this restriction, the indicial admittance 
functions of the important types of transmission systems are 
rather easily derived and expressible in terms of well known 
functions, and the essential phenomena of wave propagation 
clearly exhibited. In practice, however, we are concerned with 
lines of finite length with the voltage impressed on the line 
through a terminal impedance Zi and the distant end closed 
by a second terminal impedance Z%, We now take up the problem 
presented by such a system. 

Let K — K{p) denote the characteristic operational impedance 
of the line, and y=y{v) the operational propagation constant of 
the line. We have then 

V=Ae-'y^+Bey^, 


where A and B are so far arbitrary constants. To determine 
these constants we assume an e.in.f. E impressed on the line at 
a; = 0 through a terminal impedance Zi and the lino closed at x = s 
by a second terminal impedance ^ 2 - At x = s we have therefore 


Z2l = 7 


whence from (240) 


^ e-^” ex^B = A e-'^+Be'"’ 


and 


B= ~\j^e-’^->'A 
1 + P2 


where p^-ZifK. 

At X = 0 we have 


( 241 ) 


V=E-ZxI 

132 
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whence 

a+b^e-^a+^b, 

{l-\-pi)A-\-(l—pi)B=E, 

where pi=ZilK. 

From (241) and (242) we get 

A _1+^2_ 

(1+Pi) (1+P2) - (1-Pi) (1 - 

_g ^_ -(1- P2)e-2T» 


(242) 


' (1+Pl)(l+P2) - (1 -pi)(l -P2)e-2r« 


E 


and finally 




fi—y X _|_ 1 - Ela—y(. 2 «—id) 

E ^1+P2^ 


^+^11 _ luElp-iy^ 
i+pi 1+P2 


(243) 


If we replace F by a unit e.m.f.'we get the operational formula 
for the indicial admittance Ax; thus 




A =---- 
* K 1 




-278 


1__ 

Zxip) 


(244) 


where 


\=K/(K+Z0, 

1-px^K-Z, 
’i+Pi K+Zi 
_l-p^K-Z^ 
I'Fpb 


Pi 


P2 = 


K,y,Zi,Zi,pi and are, of course, functions of the operator p. 

The integral equation corresponding to the operational formula 
(244) is 


dx(p) Jo 


e-^^Ax{i)dt. 


(245) 


Now by (244) we can expand IJZxiv); it is 
1 , e-T'* 


Zxiv) 


^—7(28—a) 
g--7{2«-|-a:) g~7(48—a;) 

+Xpi Pi --l-XpiP22-^- 

^—7(4^4. a:) 

+ Xpi®P2*- r 


( 246 ) 


K 
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Now we observe that is simply the operational formula 

for the indicial admittance at point x of an infinitely long line 
with unit e.m.f. impressed directly on the line at a; = 0. This will 
be denoted by Similarly /K is the operational for¬ 

mula for the indicial admittance at point (2s—a;) with unit e.mi, 
impressed directly on the line at x == 0. This will be denoted by 
etc. 

Recognition of this fact allows us to derive a formal solution in 
terms of a series of reflected waves. For let a set of functions 
V 0 jVi,V 2 fVsj . . . satisfy and be defined by the operational 
equations 


yo = X(39)=X 
Vi—\fX2 
V% == X/XijU2 

— etc. 


(247) 


It then follows from the preceding and theorem II that 

[ +2^2(i —T)a2s-|-a;(T)+• ■ • 

If, therefore, we know the indicial admittance of the infinitely 
long line with unit e.mi. directly applied and if we can solve the 
operational equations (247), then Ax{t) is given by (248) by inte¬ 
gration. This solution may well present formidable difficulty in 
the way of computation. It is, however, formally straightforward 
and the numerical computation is entirely possible, the only ques¬ 
tion being as to whether the importance of the problem justifies 
the necessary expenditure of time and effort. Without any 
computations, however, the solution (248) admits of considerable 
instructive interpretation by inspection. The first term repre¬ 
sents the current at point x of an infinitely long line in response 
to a unit e.m.f. impressed at rc==0 through an impedance Zi; 
PQ=VQ(t) is the corresponding voltage across the line terminals 
proper. The second term is a reflected wave from the other 
terminal due to the terminal irregularity which exists there. The 
third term is a reflected wave from the sending end terminal, etc. 
The solution is therefore a wave solution and is expanded in a 
form which corresponds exactly with the sequence of the 
phenomena, which it represents. 
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The solution takes a particularly instructive form when 
Zi=kiK and Z 2 ~h 2 K where hi and h^. are numerics. Then 


^;o = 

Vi = 

V2^ 


1+h 
1 l-k2 
1+ki 1+fe 

_1_ 1—^2 l—fcl 

i+&i l+fc2 


j etc. 


and 


A,(t) = 


1+h 


a*(0+jip|^a24-A=(^) 


(249) 


(250) 


If ki = 0, k 2 = 1 we have the case of the e.m.f. impressed directly on 
the sending end of the line and the distant end closed through its 
characteristic impedance; the solution reduces to 


Axit)=ax{t) 


aS; of course, it should by definition. 

If ki = 0 and fc 2 = we have the case of the line open-circuited 
at the distant end, and the solution reduces to 

Ayi(i) ^ {^a;(0 <^2S“-a;(0 0^28+x{^) + + * * * }• (251) 

Finally, if both ki and /c 2 are zero, the line is shorted and 

As(t) = {ax(t)+a2,-x(t)+ai,^:,(t)+au-xit)-\ -}. (252) 

The operational equations (247) admit of further interesting 
and instructive physical interpretation without computation. 
Consider a circuit consisting of an impedance Zi in series with an 
impedance K. Let a unit e.m.f. be applied to this circuit and let 
Vt, be the resultant voltage across the impedance K. Then, 
operationally, 

so that Vo, thus defined in physical terms, is the Vo of equations 
(247). 

Now let this voltage be impressed on a circuit consisting of 
an’impedance 2 Z 2 in series with an impedance K Z 2 so that the 
total impedance is K+Z^. Let the resultant voltage drop across 
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the impedance element K—Zi be denoted by vi; then operation 
ally 

K K-Zi^ 

"" K+ZxK-{-Zi 

which agrees with Vi as given by equation (24.7), 

Similarly if voltage Vi is applied to a circuit consisting of a; 
impedance 2Z i in series with an impedance K--Zi and if Vi denot 
the voltage drop across impedance K—Zi, then 

V2 = 


We can thus see physically what the voltages VQyVi,V 2 , . , 
mean in terms of simple circuits consisting of K and Zi i 
series and K and Z^ in series respectively. 

I shall now work out a specific problem exemplifying tli 
preceding theory. The example is made as simple as possible fc 
two reasons. First because its simplicity makes it more instru( 
tive than when the phenomena depicted and the essentials ( 
the mathematical methods are obscured by complicated formuk 
and extensive computations. Secondly while the general methc 
of solution illustrated is thoroughly practical we cannot hope 1 
arrive at the numerical solutions of the complicated problen 
without a large amount of laborious computations, Problen 
involving transmission lines with complicated terminal imp 
dances are among the most difficult, as regards actual numeric 
solution, of any which present themselves in mathematic 
physics. On the other hand, the formal solution (248) gives at 
glance the essential character of the phenomena involved. 

The specific problem we shall deal with may be stated as follow 
A unit e.m.f. is directly impressed on the terminals of a transmi 
sion line of length s, the distant end of which is closed by a co 
denser Co. The line is supposed to be non-dissipative, i 
constants being inductance L and capacity C per unit lengt 
Required the current at any point a;(a;<s) of the line, 

W e write \/L/C ^k,ll\/LC—v: then by virtue of the precedi: 
analysis of transmission line propagation the indicial admittan 
a 33 of the infinitely long line is given by 


for t<x/v, 
for t>xfv. 
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The operational characteristic impedance is, of course, h = -\/LfC, 
and the terminal impedances Zi and Z 2 are given by 

Zi = 0, 

Z2 = llpC,. 

Referring now to equation (244) we have: 

X = l, Mi = l, 

_k-l/pCo_kCap-l 
^““/c+l/pCo kC,p+l 

Consequently, referring to equations (247), we have, operationally. 


Vo = l 

Vl=V2 


kCop — l 
fcCopH”! 


V3 = Vi = 


TbCop-r 

icCoP+1 


VcCop-r 

kCoP+i 


2 


s 


In order to determine these functions we have therefore to solve 
the general operational equation 


kCop-l 

kCop+l 


where F„ denotes either v^n-i or v^n- 
In order to eliminate the coefficient kCo we make use of theorem 
VI, and write 



In accordance with that theorem 

Vn(t) = <i>n{t/kCo). 

We therefore start with the operational equation 



Now the solution of this operational equation is very easy and can 
be expressed in a number of ways. We require it expressed in the 
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form most easily computed. The following appears best adapted 
for our purposes. Consider the auxiliary operational equation: 

_ _/p-2V 


The explicit solution is gotten by placing l/p™ by and p" by 
whence 

n2t n(n-l){2ty , 1 ^( 20 ". 

~ flfl'''' 2! 2! ^ nl 


But writing 


it follows that 


P~2Y 1 
P ) H{v) 


Y>+1/ iy(p+l) 

P+1 _ P _ 

p p+lif(p+l) 

. -2. -?— 


p) p+lHip+1) 


Referring now to theorem V, we see that 

= ^1+J' dt^anit) • e“‘. 

Since we have already solved for <r„(i), this determines (j>nit) 
and hence F„(<)- The functions vo,Vi,V 2 . . . are therefore 
determined. 

Now refer back to equation (248) giving the required current in 
terms of vt.iv^ . . . and the admittances axit),a 2 »-xitn), ■ • • 
It follows at once by substitution of the preceding that 

• ■ • j 

the functions va,V},V 2 being zero for negative values of the argu¬ 
ment. This result may possibly require a little explanation. 
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Consider the expression 

where 1(0 denote a function which is zero for t<h and unity for 
t>tQ, It is evidently identical with the admittance pro¬ 
vided the proper value is assigned to ^o. 

Now since 1 (0 = 0 for ^ and unity for t>h the preceding may 
be written as zero for t<U and 

ioxt>U 

which is equal to 

If we set x^Oj we get the. current entering the line; thus 

-I I 

+27,j. 

This has been computed for the case where '\/sLCo = 10 and is 
shown in Fig. 26. Referring to this figure we see that the current 



Fig, 26 . —Curreiit ontering non-dissipativo line terminated by capacity Co unit 
e.m.f. applied to line. 

jumps at i==0 to the value ^/C/L^l/k, and keeps this constant 
value for a time interval 2s/v. At this instant the first reflected 
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wave arrives and the current takes another jump of 2/k, There¬ 
after it begins to decrease very slowly until time t-is/v at which 
time it takes another jump of 2/k. Thereafter we have a series 
of jumps of 2/k at time intervals 2s/y, the current decreasing 
between successive jumps. The smooth curve is the indicial 
admittance of an oscillation circuit consisting of an inductance sL 
in series with a capacity Co. We see, therefore, that the current 
in the line oscillates with discontinuous jumps about the current 
in the corresponding oscillation circuit. Since the whole circuit 
contains no resistance, the oscillations never die away, but con¬ 
tinue to oscillate, as shown, about the curve 

which is the indicial admittance of the corresponding oscillation 
circuit. 

I shall now discuss a method of solving circuit theory problems, 
quite generally applicable to complicated networks, and particu- 



Fig. 27. 


larly useful in dealing with transmission lines terminated in 
impedances. I have found it particularly useful in arriving at 
numerical solutions where other methods prove far more laborious. 
It is also of mathematical interest, as it applies another type of 
integral equation to the problems of electric circuit theory. 

Suppose that we have a network with two sets of terminals as 
shown in Fig. 21} Now suppose that terminals 22 are short 
circuited and a unit e.m.f. inserted between terminals 11. Let 
the resultant current flowing between terminals 11 be denoted by 
Sn{f)=8n and that between terminals 22 by 
Sn is the driving point indicial admittance with respect to 
terminals 11 and Su the transfer indicial admittance of terminals 
22 with respect to 11 under short circuit conditions. 

Similarly if terminals 11 are shortcircuited and a unit e.m.f. 
inserted between terminals 22 the current flowing between termi- 
1 Regarding conventions as to signs, see the Appendix to this chapter. 
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aals 22 is denoted by S 22 {t)—S 22 and that flowing between ter¬ 
minals 11 by Si 2 {t) =Si 2 ^ If the network is passive, i.e., contains 
no internal source of energy, it follows from the reciprocal theorem 
that 821 = Si 2 * As far as the two sets of terminals are concerned, 
the network is completely specified by the indicial admittances 
^U;S22,^S21 = ^12' 

Now let a voltage Vi(t) = Fi be inserted between terminals 11, 
and a voltage 7*2(0 ="^2 between terminals 22. The current 
flowing between terminals 11, denoted by Ji is 

hit) =j^jji{T)Sii(t-r)dr+jJj2(r)Sx,(t-r)dT (253) 

while the corresponding current between terminals 22 is 

hit) =1 JV2(r)S22(4-r)dr+|J^Vi(r);S2i(t-r)dr. (254) 

Now consider two networks of indicial admittances SnjS 22 ,Si 2 
= 821 and Tu,T 22 , T 12 - T 21 respectively and let them be connected 



Fig. 28 . 


in tandem as shown in Fig. 28 to form a compound network. 
We require the indicial admittances of the compound network in 
terms of the indicial admittances of the component networks. 

Short circuit terminals 22 of the compound network and insert 
a unit e.m.f. between terminals 11. Let Vzi(t) denote the result¬ 
ant voltage between terminals 33 measured in the direction of the 
arrow, and h the current flowing between the networks. We 
have then the two following expressions for the current 13 , 

h^Snit)-^^£v3iir)S,2it-r)dr (255) 

and 

h = j^£Vnir)Txx(t-T)dr. (256) 

Equating, we get 

jJjziir)[ 822 it-r) + Tnit-r)]dT^Snit). (257) 
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By precisely similar reasoning, if terminals 11 are short circuited 
and a unit e.m.f. inserted between terminals 22, and the corre¬ 
sponding voltage across terminals 33 denoted by 732 , we have'- 

I Jj732(T)[522(i-r)-l-riia-r)]dr= ri2(i). (258) 


Equations (257) and (258) are integral equations of the Poisson 
type which completely determine Fsi and 732 in terms of the 
indicial admittances S and T. We shall discuss the solution of 
these equations presently. 

If Un,U 32 ,U 2 i=lJn denote the indicial admittances of the 
compound network we have at once 


and 


C^n = /Su(i) “ 

--I 

dt] 

^V3l(T)SJ2(^- 

'0 

-T)dr 

(259) 


_d 

dt^ 

rV32(T)T2l(i- 

/o 

-T)dr 

(260) 


U2i=Ul2- 


j‘Vn(r)T2i(t-r)dr 

■j'v,2ir)8Ut-ryir. 


d 
di 

A 

dt 


(261) 


If, therefore, equations (257) and (258) are solved for 731 and 
732 , the required indicial admittances of the compound network 
are given by (259), (260) and (261) in terms of the indicial admit¬ 
tances of the component networks. 

A simple example will now be worked out illustrating the 
method of solution just discussed. Suppose that a unit e.m.f. is 
impressed on a transmission line (infinitely long) of distributed 
constants through a terminal resistance i?o. Required 

the terminal line voltage 7. 

The operational equation of this problem is gotten in the usual 
manner. The current entering the line is 


7 


'Lp+B 


'' 7n being opposite to Vn in direction. 
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It is also obviously equal to ^ (1 — 7): equating the two expres¬ 
sions, and rearranging we get : 

iLp^ 

Lp+R 

Writing R/2L=p and setting Rt, = \/L/C, this becomes 

7 = (262) 

l + '\/l+2p/p 

This operational equation can, of course, be solved in a number 
of ways, though, as a matter of fact, its numerical solution is 
quite troublesome. This point will be returned to later; we shall 
first formulate the problem in accordance with the method just 
discussed. 

The indicial admittance of the line is known; it is 

Consequently the current entering the line is explicitly 

1 

But the current is also equal to (1 —7(i)); equating, we get 


V(t) = l-Roj^j‘v(T)A it-T)dr. 
Performing the indicated differentiations 

= l-i?o A(0)7(0 -Roj‘v(r)A’(t- 


V(t) = l-RoAiO)Vii)-Ro V(T)A'(t-T)dr 


NowA(0)-- 


and Ro^^/LfC; therefore the equation becomes 

V(t) = ^+^£v(t-r)[h{pT)-hipT)]e-^^dr. 
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As a matter of convenience we change the time scale to pt^ and get 


^i_i n 

2 Vo 


r)[7o(r)-Ji(r)]e-'dr 


drV{t-T)~e-Io{T) 

dr 


(263) 


where it is understood that t is actually pt. This is the integral 
equation of the problem and is in the canonical form of Poisson's 
integral equation. 

Before solving this equation numerically I shall show how a 
simple approximate solution is obtainable immediately; an advan¬ 
tage often attaching to this type of integral equation. 

The function ^ is equal to —1 for ^ = 0 and converges 

rapidly to zero. F(i^) has, as we know'from the operational equa¬ 
tion, the initial value 1/2 and the final value 1. Neither function 
changes sign. It follows from the mean value theorem that the 
equation can be written as 


where aSl, Integrating 

and 


(264) 


The correct initial and final values of V{t) result for all final 
values of a g 1; so that approximately 


F(0 = 


1 + e 


This equation, while not exact, except for ^ = 0 and t very large, 
shows faithfully the general character of V(t) and the way it 
approaches its final value unity. For large values of t 


whence 




e'"^Io(0 = l/'\/27r^ 
1 


1 + 1/V27r{' 


«S:8. 


(264-a) 
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Approximations of the foregoing type are not always possible 
and may not be of sufficient accuracy. I shall therefore give 
next a method of numerical solution which is generally applicable 
to integral equations of this type and works quite well in practice. 
We shall write the integral equation in the more general form 

w(^) =/(•'«) + u(x—y)k(y)dy (265) 

where/(o;) and k{y) are known and u(x) unknown. The method 
depends on the numerical integration of the definite integral. 
Let us divide the x scale into small intervals d and for convenience 
write 

u(nd) = Un 
find) =/„ 
k(nd) ===kn» 


Now from the integral equation we have at once 

^^(0)=^4o =/o 

u(d) = ui=fi+J u{d-y)k(y)dy. 

Now if d is taken sufficiently small 

^ d 

uid—y)Hy)dy = -[uiU+uo /ci] 


X 


whence 


and 


Ui =/i +2 [^1 ko^i-Uo ki] 




which determines Ui since wo is known. Similarly 
W2==/2+cZ[|W0 7<^2 + Wi/Ci+|U2/Co] 


which determines U 2 . Proceeding in the same manner 

^^3-/3+ d[iU q fcs +ui h + U2 ki + 1^3 kolj etc. 

In this way we determine the value of u(x)j point by point from 
the recurrence formula 


Wn = 


fn'i^d[2'Uo /Ctji —1 

..■‘■i- 


+ ^2^71-2+* 

‘ ■g'ifco d 


» -j-Un-^lki] 


(266) 
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The result of the application of numerical integration, in 
accordance with formula (266), to the integral equation (263) is 
shown in Fig. (29). The dotted curve is a plot of the approxi- 



Fig. 29.—^Lino terminal voltage unit e.m.f. impressed on line througH resistance. 

Ro«Vl7c, 


mate solution as given by equation (264), for ot = l. We see 
that the voltage starts with the value 1/2 and slowly reaches its 
ultimate value, unity, its approach to unity, for large values of 
being in accordance with the formula 

7 ( 4 ) =- 

l + llV2Tt 

The application of the foregoing method to the transmission 
line problem proceeds as follows. Let Suit), S 22 {t) and Suit) 
be the short indieial admittances of the line. Sii{t) is the current 
entering the line (at rr = 0) with unit e.mi. directly impressed 
and the distant end short circuited. Suit) is the current at 
x-s under the same circumstances. Consequently from (252) 


/Su(i) =ao(i)+2a2a(i)4-2a4®(i)+• ■ • 

Suit) — 2{aXt) + a3«(^)+aBa(i) + • • •} • 

(267) 


S 22 is clearly equal to Snhy symmetry. 

Now suppose that an e.m.f. E = Eii) is impressed on the line 
at oj-O, ^=0, through a terminal impedance Zi^ and the distant 
end (rr=s) closed through an impedance Z 2 . We suppose these 
terminal impedances and the actual impressed e.m.f. replaced 
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by the actual line voltages Vi and Fs, impressed directly on the 
line at a; = 0 and Sit x==s; the terminal currents are 


j Sl2(t — T)Vl(r)dT 

S22(t — r)V2(T)dT. 


(268) 


(269) 


But the current at tc = s is also equal to the current in the terminal 
impedance Z 2 in response to the terminal voltage F 2 : denoting 
by 0 ( 2(0 indicial admittance of Z 2 it is 

I.(t) = j^£a2(t-T)V2ir)dr. (270) 

Similarly the current entering the line at a;=0 is the current 
flowing in the terminal impedance Zi in response to the e.m.f. 
E — V^. Denoting by ai{t) the indicial admittance of Zi, it is 

{Eir) - Fi(r) }dr. (271) 


Equating equations (268) and (271) and (269) and (270) we 
eliminate Jo(0 and L{t) and get 


rkn(^~T)+t^i(«~r)]Fi(r)^?r-- rSi2(t-r)V2(r)dT 

Jo Jo 

= J* ai{t--r)E(j)drj 

- r!si 2 (^-r)Fx(T)dT+ r[S22(i-r)--a2(^-T)]F2(T)(^r-0. 

Jo Jo 


(272) 

(273) 


These two equations are simultaneous integral equations of the 
Poisson type in Vi and F 2 , which completely determine these voltages 
provided the admittances and the impressed voltages are known. 
They therefore represent the application of a new type of integral 
equation to the problem of electric circuit theory. 

The numerical solution of the general case, either by (248) or 
(272-273) is necessarily laborious when the terminal impedances 
are complicated and is only justified when the technical impor¬ 
tance of the problem is considerable. I wish, however, to empha* 
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size two points in this connection: the numerical solution is 
always entirely possible and, compared with other and older 
forms of solution, enormously simpler. One has only to inspect 
the classical forms of solution of problems of this type to realize 
the truth of this last statement. 

I shall now give two applications of equations (272-273) to 
specific problems, in one of which an approximate solution of the 
integral equation can be gotten, and in the other of which numer¬ 
ical integration is applied. 

Problem L Given a non-inductive cable of distributed con¬ 
stants C and R and length s, with unit e.m.f. applied at a: = 0, 
while at 05 = 8 the cable is closed by a condenser Co. Required 
the terminal voltage F(0 across the condenser Co. 

We first write down the short-circuit indicial admittances of 
the cable; from equation (168) of a preceding chapter and equation 
(267) they are:— 

Snit)^S22{t) 

4(3 IG^ 36/3 ^ 

l~f-2c i ~j-2e ir ^ ~(" ■ * |? (274) 

where /3 =s“jBC'/4. 

N ow the current at a: = s is equal to 

Si,(t)-^^£v(T)S2,(t-r)dr. 

It is also the condenser current due to the voltage V (t) ; that is 

c,^yw. 



Equating the two expressions and integrating we get 

CoV(t)= rSi 2 (r)dT- rV(r)S 2 Ut~r)dr (27()) 

Jo Jo 

which is the integral equation of the problem. In order to get 
an approximate solution without detailed computation we assume 
that the cable is long. In this case the leading terms of (274) 
and (275) are large compared with the terms following: Further¬ 
more Sis(t) builds up very slowly while Siz(i) is a rapidly varying 
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function. A good approximation therefore results if we take 
V (r) outside the integral sign in (276) and write 


whence 





V(t) 322(r)dr 



(277) 


This approximation is quite good for long cables and shows the 
way y (t) builds up quite truthfully. We see that V is initially 
zero, and builds up ultimately to unity. For large values of t, it 
becomes 


y(f) = 



(278) 


This is the approximate formula also for the open circuit voltage, 
as may be seen by setting Co = 0 in (277). 

In electric circuit problems, it is often sufficient, as implied 
above, to Icnow qualitatively the behavior of an electric system 
without going through the labor of detailed computation. For 
this purpose the formulation of the problem as a Poisson Integral 
Equation is particularly well adapted. A simple example 
will be given, which can be checked from the known solution. 
Suppose that we required the voltage V at point x of an infinitely 
long transmission line (L,R,C) in response to a unit e.m.f. 
impressed at x=0. This is, of course, known from formula 
(211-a): we shall here be concerned, however, with the approxi¬ 
mate solution from the Poisson integral equation of the problem. 

If ax(t) denote the indicial admittance of the line at point x, 
then the current at point x is simply ax(t), which is given by form¬ 
ula (210-a). But if V(t) is the voltage at point x, the current is 
also given by 



Equating these two expansions, we get the integral equation of 
the problem 
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xB IC 

Now if we write r = pf—Awherep = 22/2LandA==—-J|-)then 

a«=^|e-tT+A)7„v’r(2^1), T>o, 


and in terms of the relative time T, the integral equation is 
reducible to 

4 rV(T-T)e--Io(r)dT = e-(^+^)Io(VT(^Ij 

dl Jo 

while the exact formula for V is by (211-a) 

7{T) = e~^+Ae-^ 

Jo T(r+2A) 


From the integral equation it is easy to establish superior and 
inferior limits for V{T); it is 


7(r) 


rT 

.Jo 


^hiT)Va(T)dT 


X T 

e-^Io{T)dT 


= n(r). 


Both formulas give the correct initial and final values of V ; 
namely and unity. Since V lies between Va and Vt, the 
mean value (7a+y6)/2 also has correct initial and final values 
and should be a better approximation than either. The table 
given below shows the orders of approximation obtainable from 
the case where A =3. It is evident from this table that the 
foregoing approximate formulas exhibit the form of V{T) qualita¬ 
tively in a quite satisfactory manner. 


T 

Fa 

n 

UVa+V,) 

7 

0 

0.05 

0.05 

0.05 

0.05 

2 

, 0.25 

0.12 

0.18 

• 0.17 

4 

1 0.39 

0.19 

0.29 

0.26 

6 

0.60 

0.23 

0.36 

0.32 

8 

0.57 

0.27 

0.42 

0.37 

10 

0.64 

0.31 

0.47 

0.41 

12 

0.69 

0.34 

0.51 

0.44 

15 

0.74 

0.38 

0.56 

0.48 

18 

0.78 

0.41 

0.60 

0.52 
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Prohlem IL Our second illustrative problem may be stated 
as follows:—A unit e.mi. is impressed on a transmission line of 
length s and distributed constants At = s the line is 

closed by a resistance i?o in parallel with an inductance Lq. 
Required the current in the terminal resistance. If Vit) denotes 
the terminal voltage, the current at a: = s is given by 

It is also equal to the current flowing into the terminal impe¬ 
dance; that is 

1 1 Vv(r)dr. 


Equating and rearranging 




+^22(0) 


F(0=5i2(0- 




+^^22 


dr. (279) 


N’ow the short circuit admittances 3^^ and 3 12 are given by formula 
'210-a) of a preceding chapter, and Sn{0) — '\/C/L. In order 
;o apply numerical integration to (279), numerical values must 
36 assigned to the constants. We take 

i?o = V" L/C = 1935 ohms, 

I/O = 0.4 henry, 

I-,-29^ 

^; = l/Vi:-O = 1.105Xl0^ 
s = 100. 


The results of the numerical evaluation of equation (279), with 
hese values inserted, is shown in Fig. 30. The voltage is identi- 
ally zero until t)i=100; i = 100/ti is the time of propagation of 
he line. At that instant it jumps to the value = 

,nd then begins to die away rapidly due to the draining action 
f the inductance. The effect of secondary reflection is insignifi- 
ant and therefore not shown. The current in the terminal 
esistance is 7/i?o so that it is given by the same curve. 

I have reserved until the last the exposition of the expansion 
'leorem solution as applied to transmission lines with terminal 
npedanees, for the reason that it is the least powerful and the 
lost restricted, although most closely resembling the classical 
irm of solution. Furthermore, it does not represent the 
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sequence of physical phenomena, in fact it is not a wave solution, 
but a solution in terms of normal or characteristic vibrations. In 
practical application its usefulness is restricted to the non- 
inductive cable. 



Fig. 30.—Voltage across terminal impedance on smooth line. 


It will be recalled that the expansion theorem solution is 
formulated as follows:— 

If 


A = l/Z{v) 


is the operational equation of the problem, then the explicit 
solution is 


1 


where Pi, 2>2 * . • are the roots of the equation Z(p) =0. 

Let us apply this formula to the case of a line of length s, with 
unit e.m.f. directly applied at and line short circuited at 
a? = 5. Referring to equation (244) and putting X=mi=P 2 = 1 
we get 


1 cosh 7 (^ — 0 ?) _ ^ 

K sinh ~ Zx{v) 


(280) 


as the operational formula of the problem, 
as 


A,-((7p+(?) 


c osh y{s—x) 
7 sinh 7 S 


This can be written 




( 281 ) 
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where in the general case, 

T = (282) 

The values of y for which Zx{'p) vanishes are the roots of the 
transcendental equation 

sinh 7 S = 0 

excluding zero. These roots are infinite in number: Let ym be 
the root; then 

977 TT 

= i - 7 m = lj2, .... 00. (283) 

5 


The corresponding values of are then gotten by solving (282) 
for p and writing y —ym- 

The explicit solution of the operational equation (281) is then 




cosh ym^ 


“’'•‘dp. 


(284) 


Let us apply this tojth^non-inductive, non-leaky cable in which 
L=G=0 and y = \/RCp, so that 

and 

dy m jKO 


Also Zx(0) ==5i2. We thus get 


Aaj(j/) 



WTT 

cos— x.e 
s 


smc 


w = 1 


(285) 


This is a thoroughly practical formula for computation^ owing 
to the rapid convergence of the series. In fact, for this particular 
line termination chosen, it is probably the simplest and most 
easily computed form of solution. These advantages depend, 
however, strictly on two facts^ First, the fact that the line is 
taken as non-inductive and secondly, that the terminations 
chosen are those of a short circuit. In fact, as we shall see, it is 
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only in the case of the non-inductive cable that this type of 
solution is of any practical value. 

There is one other point which should be carefully observed 
in connection with this solution (285). This is that it is not 
expressed in terms of a series of direct and reflected waves, 
corresponding to the sequence of physical phenomena, but in 
terms of normal or characteristic vibrations. This point will be 
returned to later. 

Let us now attempt to apply this type of solution to the 
transmission line, L^R,CjG, Writing 

R , G 

R G 
2L 2C’ 

v = llVW 

we have 

whence 



Setting G = 0 for simplicity and substituting in (284) we get, 
after easy simplifications, 






mT 


s 


If we write 
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(286) can be written as 


, 1 ,vC' 


Q-pl 


rmr 


sin —{fx^vt- 


s I 


Mm" 


(287) 


This type of solution is often referred to as a wave solution and the 
component terms of the series regarded as travelling waves. As 
a matter of fact it is a solution in terms of normal or characteristic 
vibrations, each of which is to be regarded as instantaneously 
produced at time t = 0. The solution in terms of true waves has 
been fully discussed in the preceding. 

Formula (287) is practically useless for computation on account 
of the slow convergence of the series (the series are only condi¬ 
tionally convergent), and cannot be interpreted to bring out the 
existence of the actual direct and reflected waves and the physical 
character of the phenomena it formulates. In fact, as stated 
above, this form of solution is useful only in connection with the 
non-inductive cable. 

In the cases considered above we have taken the simplest 
possible terminations—these of short circuits in which case the 
roots of Z{'p) are easily evaluated. If, however, the line is 
closed by arbitrary impedances, the case is quite different, and 
the location of the roots becomes, except for simple impedances, 
and then only in the case of non-inductive cable, practically 
impossible. While, therefore, the expansion theorem solution 
can be formally written down, its actual numerical evaluation is 
a practical impossibility, except in a few cases. For this reason 
it will not be considered further here. 

The physically artificial character of the expansion solution, 
as applied to transmission lines, may be seen from the following 
considerations. When a wave is sent into the line, for a finite 
time equal to the time of the propagation of the line, it is inde¬ 
pendent of the character of the distant termination. Yet in the 
expansion solution every term involves and is dependent upon 
the impedance constituting the distant termination. Evidently, 
from physical considerations, the series of component vibrations 
making up the complete solution must therefore so combine as 
to annihilate the effect of the distant termination for a finite 
time. The solution is, therefore, mathematically correct but 
physically artificial. 
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Note on Integral Equations. 

An integral equation is defined as an equation in which the 
unknown function occurs under a sign of integration; the process 
of determining the unknown function is called solving the 
equation. 

•Integral equations are of great importance in mathematical 
physics and in recent years very considerable work has been done 
on them from the standpoint of pure analysis. 

The types of integral equations with which we are concerned 
in the present work are Laplace^s equation 


F(p)^- 




-pt 


mdt 


and Poisson’s equation 

cj}(x)=fix)+4>{y)K{x-y)dy. 

But little work has been done on Laplace^s equation from the 
standpoint of pure analysis; its most extensive and useful appli¬ 
cations appear to be in connection with the operational calculus. 
Practical methods of solution are extensively discussed in the 
text. 

We shall now briefly discuss the solution of Poisson’s equation. 
The formal series solution, which is absolutely convergent, 
is obtained by successive substitution. Thus suppose we write 

— • • • 

and define the terms of the series in accordance with the scheme 

0fl(a;) =/(x), 

<l}i{x) = jo 4>a{y)K{x-y)dy, 

X X 

4>i{y)K{x-y)dy, etc., 

the resulting series satisfies the integral equation and is absolutely 
convergent. It is, however, practically useless for computation 
or interpretation. 



TEE FINITE LINE WITH TERMINAL IMPEDANCES 157 


A power series solution, when it exists, can be gotten by 
repeated differentiation; thus 

^( 0 )=/( 0 ), 

^'(x) (x)+- y)K(y)dy, 

^'(0) =/(0)+.^(0)^(0) 

In this way all the derivatives at a; = 0 are calculable; let them 
be denoted by 0o, <#>2 . . . Then 

This form of solution, also, is of limited practical usefulness, 
except for small values of x, 

A number of mathematicians, including Whittaker and 
Bateman, have studied the question of numerical solution and 
suggested other processes. After quite extensive study of the 
question, however, the writer is of the opinion that point-by¬ 
point numerical integration like that discussed in the text is, in 
general, the most practical, rapid and accurate method of numer¬ 
ical solution. This judgment is confirmed by G. Prasad who, in 
a paper on the Numerical Solution of Integral Equations 
before the International Mathematical Congress (Toronto, 
1924), discusses the whole question and arrives at the same 
conclusion. 

In the text, numerical integration is carried out in accordance 
with Simpson^s Rule. It is possible, of course, to employ more 
complicated and refined formulas for approximate quadrature. 
It is the writer^s opinion that this is hardly justified in practical 
problems and that the required accuracy is more simply obtained 
by employing smaller intervals. 

APPENDIX TO CHAPTER IX 
Note on Conventions as to Signs in Networks. 

In the network shown on page 140 the arrows indicate the 
directions chosen as positive in the network itself, quite regardless 
of the presence of any e.m.fs. and currents. 

The sign attributed to a current, an e.mi., or a voltage is 
positive if the current, e.m.f., or voltage is in the positive direc¬ 
tion; otherwise the sign is negative. 
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Stated more fully; 

A current at a specific point (at a specific instant of time) is 
positive if it is flowing in the positive direction; negative if 
flowing in the negative direction. 

An e.m.f. or a voltage between two points is positive if the 
potential increases in the positive direction between the two 
points; negative if the potential increases in the negative direc¬ 
tion. (It may be noted that this convention makes the sign of a 
voltage the same as the sign of that e.m.f. which could be inserted 
between the two points without producing any effects in the 
network.) 



CHAPTER X 


INTRODUCTION TO THE THEORY OF VARIABLE 

CIRCUITS! 

In the preceding chapters it has everywhere been assumed 
that the networks are invariable: that is to say, that the constants 
and connections of the network do not vary or change with 
time. In many important technical problems, however, we 
wish to know, not merely what happens when an electromotive 
force is applied to an invariable network, but the effect of sud¬ 
denly changing a circuit constant or of introducing a variable 
circuit element. In the present chapter we shall show that this 
type of problem can be dealt with by a simple extension of the 
methods discussed in the preceding chapters. 

The simplest and at the same time one of the most technically 
important problems of this type is the effect of sudden short 
circuits and sudden open circuits on an energized network or 
system, This type of problem will serve as an introduction to 
the more general theory. 

The Sudden Short Circuit.—Consider the network shown in 
Fig. 31. This network, which for generality is supposed to consist 


1 


a 

- r - 


E(t) 

A 

|vtt) 

B 

1_ 


— J — 1 



Fig. 31. 


of two parts A and B, indicated schematically, is energized by an 
electromotive force E{t) which produces a voltage y(0 between 
the points ab. The voltage V{t) is calculable by usual methods 
from Eit) and the constants and connections of the network, 
supposed to be specified. 

1 The material in this chapter is largely taken from a paper by the wi'iter on 
Theory and Calculation of Variable Electrical Systems,^’ Phys. Rev. 
Feb. 1921. 
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We now suppose that, at reference time i = a short circuit is 
suddenly placed across ab; and require the effect of this short 
circuit on the distribution of currents in the network. The 
solution of this problem is based on the following proposition: 

The effect of the short circuit is precisely the same as the insertion 
at time t = 0 of a voltage—V{t), equal and opposite to V{t)j between 
points a and b. 

The resultant currents in the system for i ^ 0 are then composed 
of two components:— 

(1) The currents which would exist in the invariable network, 
in the absence of the short circuit, due to the impressed source 
E{t), These are calculable by usual methods. 

(2) The currents due to the electromotive force V{t) inserted 
at time i — O, between the points a and b. These are also calcu¬ 
lable by usual methods, since V{t) is itself known from the 
primary distribution of currents and charges. 

By the preceding analysis , we have succeeded, therefore, in 
reducing the problem of a sudden short circuit, to the determina¬ 
tion of the currents in an invariable network in response to a 
suddenly impressed electromotive force: that is, the problem 
to which the preceding chapters have been devoted. 

The Sudden Open Circuit.—The problem of a sudden open 
circuit in any part of a network can be dealt with in a precisely 



Fig. 32. 

analogous manner, although the actual calculation of the result¬ 
ant current and voltage distribution is mathematically more 
complicated. Consider the network shown in Pig. 32. Here 
the network is supposed to be energized by an electromotive 
force E(t) which produces a current 1(0 in the invariable network 
in branch ab. We require the effect of suddenly opening this 
branch. The solution of this problem depends on the following 
proposition. 

The effect of opening branch ab at reference time t — Q is the same as 
suddenly inserting at time i = a voltage V(t) which produces in 
branch ab a current —^7(0 equal and opposite to the current which 
would exist in the branch in the absence of the open circuit. 
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While this proposition is precisely analogous to the correspond- 
ng proposition in the case of a sudden short circuit, it does not 
'xplicitly determine the voltage V{t)j which must be calculated 
bs follows: 

Let the driving point indicial admittance of the network, as 
leen from branch ab be denoted by Aah{t)^ Then, from the 
)receding proposition, it follows at once that Vit) is given by 



rhis is a Poisson integral equation in V{t), from which V(t) is 
3alculable. With V(t) determined, the currents in any part of 
}he network are calculable by usual methods, and consist of two 
components:— 

(1) The current distribution in the network due to the 
impressed source E{t) in the absence of the open circuit. 

(2) The current distribution due to the electromotive force 
V(t) inserted in branch ab at time i = 0. 

As in the case of the sudden short circuit, we have thus reduced 
the problem of a sudden open circuit to the determination of the 
current distribution in an invariable network in response to a 
suddenly impressed electromotive force. 

Variable Circuit Elements*—In the preceding cases of sudden 
open and short circuits it will be observed that the network 
changes discontinuously from one invariable state to another. 
A more general case, and one which includes the preceding as 
limiting cases, is presented by a network which includes a 
variable circuit element: that is, a circuit element which varies, 
continuously or discontinuously, with time. A network which 
includes such a variable circuit element will be called a variable 
network. Variable circuit elements of practical importance are 
the microphone transmitter, which consists of a variable resist¬ 
ance, varied by some source of energy outside the system; the 
condenser transmitter, which consists of a condenser of variable 
capacity; and the induction generator, in which the mutual 
inductance between primary and secondary, or stator and rotor, 
is varied by the motion of the latter. The case of a variable 
resistance will serve as an introduction to the general theory 
of such variable networks. 

Consider a network, energized by a source E{t) in branch 1, 
and containing a variable resistance element r(f) in branch n. 
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The functional notation r{t) indicates that the resistance r varies 
with time. Let In{t) denote the current in branch and assume 
that the network is in equilibrium prior to the reference time 
i = The mathematical theory of this network depends on the 
following proposition:— 

The network described above can be treated as an invariable network 
by eliminating the variable resistance element r{t) and inserting an 
electromotive force lhat iSj an electromotive force equal 

and opposite to the potential drop across the variable resistance 
element. Consequently the current in the variable resistance 
branch is determined analytically by the integral equation 

J.(i) {r)h{r)AUt-r)dT. ( 288 ) 

The first component is simply the current I^{i) which would exist 
in the variable branch if the variable element were absent; 
hence, dropping the subscript n for convenience, the current in 
the variable branch is given by the integral equation 

I(i)=/o(t)-|J^‘r(T)7(r)A(«-T)<iT • (289) 

and the voltage across the variable element by 

v{t)^r{t)I{i), (290) 

Having determined I{t) and v{t) from this integral equation, the 
distribution of currents in the network is calculable as that due to 
a source E{t) in branch 1 and a source v{t) in branch n of the 
invariable network: that is, the network with the variable resist¬ 
ance element eliminated. 

A very simple example will serve to illustrate the foregoing:— 
Into a circuit of unit resistance, and inductance L = l/a, in 
which a steady current /o is flowing, a resistance r is suddenly 
inserted at time t=0: required the resultant current I(;t), In 
this case we have: 

A(0=indicial admittance of unvaried circuit 
r{t)=r, 

and the integral equation of the problem is: 

m =-’’I (< - y)‘^y 

=Iq— raj' I{t—y)e~‘‘’‘dy. 
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If the solution is carried out as indicated by (291) below, and if 
the notation is introduced, we get without difficulty 

I(t) = Jo { l—r{l — ei{x)e-^)+r^(l — e 2 (a;)e“®)-r^(l —e 3 (aj)e~*) + 

■ ■ ■ }> 

where the function e„(a:) is defined as: 

en{x) = H-a:/l!+a:V2!+a;V3H- ■ • • +a:”-V(n-1)! 

= first n terms of the exponential series. 

For all finite values of the resistance increment r the series can be 
summed by aid of the identity 

1—e„(a:)e“''= f 


Substitution of this identity gives 


1(0 = Jo(^l-r 


= J0' 


l+re’ 


-(l+r) X 


1+r 


Equation (289) is an integral equation of the Volterra type, 
which includes the Poisson integral equation as a special case. 
Its formal series solution is obtained as follows:—Assume a series 
solution of the form 

• (291) 

and define the terms of the series by the scheme 
Ji(0=-^j' r{T)Io{T)A{t-T)dT 

. (292) 

I:fc+i(i) =1 JjKr)//.(T)A(f-T)dr. 

Direct substitution shows that this series satisfies the integral 
equation. Furthermore, it is easily shown that it is absolutely 
convergent. 

While this series solution is not, in general, well adapted for 
numerical calculations, it throws a good deal of valuable light on 
the ultimate character of the oscillations in the important case 
where E{t) and rif) both vary sinusoidally with time. In this 
case, if the frequency of the applied e.m.f. be denoted by F and 
that of the resistance variation by /, it is easy to show that the 
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current Jo(0 in the unvaried circuit is ultimately^ a steady state 
current of frequency F. This follows from the fact that the 
definite integral which defines the current Io(i5) is resolvable 
into the ultimate steady state current corresponding to an 
applied force of frequency F, and the accompanying transient 
oscillations which ultimately die away. The fictitious e.m.f. 
which may be regarded as producing the component current 
Ii{t) is this is ultimately the product of the two fre¬ 

quencies F and and therefore resolvable into two terms of 
frequency and F-'f respectively. Carrying through this 
analysis, it is easy to show that each component current is 
ultimately a steady state but poly-periodic oscillation, as indi¬ 
cated in the following table: 

Component Current Frequency 

/o. F, 

h . 

I, . F+2/,F,F-2/, 

h . F+3/,F+/,F--/,F-3/, 

U . F+4/, F+2/, F, F - 2/, F ~ 4/. 

It is of importance to observe that the component currents 
involve, from a mathematical standpoint, multiple integrals of 
successively higher orders, the nth component In{t) involving a 
multiple integral of the nth order with respect to Conse¬ 

quently the successive currents require longer and longer inter¬ 
vals of time to build up to their proximate steady state values, 
so that the time required for the resultant steady state to be 
reached cannot be inferred from the time constant of the unvaried 
circuit. 

From the preceding table it will be seen that the ultimate 
steady state current is obtained by rearranging the series lo+Ii-f 
I 2 and is of the form 
- 1-00 

An cos Bn sin (O-l-nct))^ 

n= — 00 

where 12= 2TrF and co = 27r/. 

lit hardly seems necessary to remark that the reference time ^=0 is 
purely arbitrary and that the resistance variation may start at such a time 
thereafter that Io(t) may be regarded as steady state during the entire 
time interval in which we are interested. Going farther, if we confine our 
attention to sufficiently large values of the whole process may be^treated 
as steady state. 
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It is interesting to note that this series comes within the defini¬ 
tion of a Pourier series only when F = 0 or an exact multiple of /. 
The steady state solution is of very considerable importance and 
is considered in more detail in a succeeding section. 

From the foregoing we deduce an outstanding distinction 
between the variable and invariable networks. In the latter the 
currents are ultimately of the same frequency as the impressed 
e.m.f., whereas in the former they are ultimately of an infinite 
series of frequencies. 

In the preceding example, the variable impedance element is a 
resistance T(t), If the variable element is taken as an inductance 
\{t) the voltage, corresponding to equation (290) is 

The case of a variable capacity element is handled as follows: 
Let 1/(7 = 5 and assume that S is variable; thus, AS = ASo+cr(i). 
The drop across the variable condenser element is then 

Similarly a variable mutual inductance m( 0 between branches m 
and n produces the voltages 

in branch m, and 

in branch n. This case may be illustrated by: 

The Induction Generator ProhUyn 
In a sufficiently general form, this problem, which includes the 
fundamental theory of the dynamo, may be stated as follows: 

Given an invariable primary and secondary circuit with a 
variable mutual inductance Mf{t) which is an arbitrary but 
specified time function, and let the primary be energized by an 
e.m.f, E{t) impressed on the circuit at the reference time / -O: 
required the primary and secondary currents. 

In operational notation the problem may be formulated by the 
equations: 

“-pM/(^)Ji+2’22 l2 = 0, 
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in which Zu and Z 22 are the self impedances of the primary and 
secondary respectively; Mf(t) is the variable mutual inductance; 
E{t) is the applied e.na.f. in the primary; and p denotes the 
differential operator d/dt By aid of the fundamental formula 
these equations may be written down as the following simultane¬ 
ous integral equations; 

hit) ^j^jUyA,,it-y) (siy) + M£[fiy)hiy)]^ 
hit) =Mf^j*dyA22it-y)^[fiy)hiy)]. 

In these equations, Aii{t) and A 22(0 denote the indicial admit¬ 
tances of the primary and secondary circuits respectively (when 
M = 0): that is, the currents in these circuits in response to a 
unit e,md. (zero before, unity after time t — Q), We assume, of 
course, that they are known or can be determined by usual 
methods. 

It follows at once that the formal solution of these equations is 
the infinite series: 


Ii(t) =Xo(i5) +X2(^) +X4(0 + '' ' +X2w(0d” * • • 
l2(0 = I"l(0+J"3(0 + y'5(0+-** 

in which the successive terms of the series are defined as follows: 


Xo{t) 


=-r 

^Uo 


dyA nit-y)Eiy)=hit), 


Fi(<) =m| J Vn(i-2/) |[/(2/)Xo(2/)], 
X,ii) =M^J‘dyAi,it-y) j- [/(y) Fi(2/)], 


Y^it)=M- 


X 


dyA 22 it-v)-^ylfiy)X 2 iy)], etc. 


In the light of formula 

I it) = “ y'^^y 

the physical interpretation of the series solutions follows at once: 
thus, Xnit) is equal to the current loit) flowing in the isolated 
primary in response to the applied e.m.f. EQ :); the first component 
current Yiit) in the secondary is equal to the current which 
would flow in the isolated secondary in response to the applied 
e.m.f. Mid/dt)fit)Xait)-, XaCO, the second component current 
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in the primary, is equal to the current in the isolated primary in 
response to the applied e.m.f. M(d/dt)f(t)Y^{t); etc. The 
resultant currents are thus represented as built up by a to-and- 
fro interchange of energy between primary and secondary, or 
by a series of successive reactions. In the important case where 
the applied e.m.f. and the variation of mutual inductance are both 
sinusoidal time functions, of frequency F and / respectively, it is 
easy to show that each component current becomes ultimately 
equal to a set of periodic steady state currents. Thus the com¬ 
ponent Xq is ultimately singly periodic, of frequency F; Yi is 
ultimately doubly periodic, of frequencies F+f and F—f; 
triply periodic, of frequencies i^+2/, F and F—2f; Ys quadruply 
periodic, of frequencies F+3/, F+/, F—f, F—Sf; etc. 

The Solution for the Steady State Oscillations.—^For the very 
important case of periodic applied forces and periodic variations 
of circuit elements we are often concerned exclusively with the 
ultimate steady state of the system, and not at all with the 
mode in which the steady state is approached: that is, attention 
is restricted to the periodic oscillations which the system executes 
after transient disturbances have died away. In this case, if 
the periodic variations of circuit elements are sufficiently small, 
the required steady state is obtained in the form of a series by 
replacing each term of the complete series solution by its ultimate 
steady state value; a process which is very simple in view of the 
physical significance of each term of the latter series. The 
appropriate procedure will be briefly illustrated in connection 
with the variable .resistance element. In view of the fact that 
we are concerned only with the ultimate steady state oscillations, 
we can base the solutions on the symbolic equation 

(293) 

Hero r(i) is the variable resistance element; To is the current which 
would flow in the absence of the resistance variation; and Z is a 
generalized impedance of the network, as seen from the variable 
branch. Its precise significance and functional form is given 
below. 

We now suppose that lo is given by 

(real part) 

= K-^oe’“+/Qe“'”‘) 


(294) 

(295) 
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where the bar indicates the conjugate imaginary of the unbarred 
symbol, so that (295) is entirely real. Correspondingly the 
variable resistance will be taken as 



— ioyt^ 


— (real part) 

= r cos o)t. 


(296) 


Here r is taken as a pure real quantity, which fixes the size of the 
resistance variation. No loss of generality is involved in this, 
since it merely involves referring the time scale to the zero of the 
resistance variation. 

The symbolic impedance Z, as employed in the theory of alter¬ 
nating currents, will depend on the frequency and is, in general, a 
complex quantity. Its value at frequency Q/27r will be denoted 
by 

z{im) ==-Zq 

while its value at frequency (0“t“^^)/27r will be written as 
2 ('l ( 0-f- '^^w) ) =-Zn. 

We now assume a series solution of (293) of the form 


J=J0+/Z + J2+-- 


where the terms of the series are defined by the symbolic 
equations 



T{t) 

z 


0 


/n-hl — 


r{i) 

Z 


In. 


(297) 


Substitution shows that this series formally satisfies the equation. 

Starting with the first of (287) and substituting (295) and (290) 
we get 


"" ~lz i ' (298) 



—<*))£ 1 


( 299 ) 
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In (299) it is to be understood that the real part is alone to be 
retained. 

Proceeding in a similar way with the equation 

r _ ’"(0 r 

i 2- h- i 1 


we get 


Jo 


gi(Q4-2£«>){ QiQ,t / 2 


(300) 


In this way the steady state series solution is built up term by 
term, the component currents being poly-perioclie as indicated 
in a previous table on p. 164. 

For sufficiently small impedance variations this method of solu¬ 
tion works very well, and leads to a rapidly convergent solution. 
In, other cases, however, the solution so obtained may be diverg¬ 
ent, even when the complete series solution from which it is 
derived is absolutely convergent. The explanation of this lies 
in the fact that the steady state series so obtained is the sum of 
the limits (as t approaches infinity) of the terms of the complete 
series solution, whereas the actual steady state is the limit of the 
sum. These are not in general equal; in particular the former 
may be and often is divergent when the latter is convergent. 

In view of the foregoing considerations it is of importance to 
develop another method of investigating the steady state oscilla¬ 
tions which avoids the difficulties in the formal series solution. 
The following method has suggested itself to the writer and works 
very well in cases where the previous form of solution fails. It 
should be stated at the outset, however, that the absolute con¬ 
vergence of the solution to be discussed, while reasonably 
certain in all physically possible systems, has not been established 
by a rigorous mathematical investigation, which appears to 
present very considerable difficulties. 

We start with the problem just discussed and, in view of the 
results of the formal series solution there obtained, assume a 
solution of the form: 


-V 

N 

= m '*’’”"** (real part). 


~N 


(301) 

(302) 
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Here the series is supposed to extend from m = +N tom ——N. 
Ultimately, however, N will be put equal to infinity. As before, 
the bar indicates the conjugate imaginary of the unbarred symbol 
and (301) is therefore entirely real. 

If we now substitute (301) in the symbohc equation (293) we 
get, by (295) and (296), 




■iM 


_r 

"2Z 






(303) 


Simplifying this equation and dropping the conjugate imagi- 
naries gives:— 






(304) 



gi(£2 + Cm — l)a))i^ 


Finally, if we write 
and 




(305) 


and equate terms of the same frequency on the two sides of the 
equation, we get 

An — — KnAn—i 

~*^m(A7n-l + -^»n+l) 0<|m|<iV* (306) 

A 0 = «/o ~ Ao (-4-i l) • 

It will be observed that, by (305), starting with An each coeffi¬ 
cient is determined in terms of the coefficient of the next lower 
index. Thus: 


An^ '^HnAn^I 

An--1= --/l7/_l(Ajv-2 + ^i7) 
hN^t ^iV—2 


Similarly 


-^AT—2 — 


1 ~ Jin^Jin 

Hn-zA n^z 
1 — AiV-2 


(307) 


^Ar-2 n.N-1 7 j. TT 
I — /ljv-1 Hn 

Continuing this process it is easy to show that, for positive indices 
(m positive), 

An.^-hmCmAn.^1 (308) 
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where Cm designates the continued fraction 


1 


1 




_^1_ 

1 /Iat—iAjv 


The procedure for the coefficient is precisely similar. For 
convenience we write and == In 

this notation we get by precisely similar procedure 


Am — K mC mA m—l 

where C^m designates the continued fraction 

1 


m •*' m-fl 


m4-2» 


(309) 


1 _ 


We now put the index N equal to infinity and the continued 
fractions Cm and C'm become infinite instead of terminating 
fractions. 

Collecting formulas we now have 


and 

whence 


j4.t) 


■ CmAm—l 
-V O' A' , 

f*' 771 mAX TM—l 


A 


Aq — Jq — (ho AIA ho'A'i) 
Jo 


1-hohiCi-ho'hi' Cl' 


(310) 


The coefficients are thus all determined in terms of Jo. 

The practical value of this method of solution will depend, of 
course, on the rate of convergence of the continued fractions. 
While no rigorous proof has been obtained, it is believed that they 
are absolutely convergent for all physically possible systems, but 
this question certainly requires fuller investigation. Neverthe¬ 
less any doubt regarding the convergence of the solution need 
not prevent the use of the method in a great many problems 
where physical considerations furnish a safe guide. For example 
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this method of solution, when applied to the problem of the 
induction generator, discussed above, leads to the usual simplified 
engineering theory of the induction generator and motor, besides 
exhibiting effects which the usual treatment either ignores or 
fails to recognize. 

Non-linear Circuits.—In the previous examples discussed, the 
variations of the variable circuit elements are assumed to be 
specified time functions, which is the same thing as postulating 
that these variations are controlled by ignored forces which do 
not explicitly appear in the statement and equations of the problem. 
We distinguish another type of variable circuit element, where 
the variation is not an explicit time function but rather a func¬ 
tion of the current (and its derivatives) which is flowing through 
the circuit. For example, the inductance of an iron-core coil 
varies with the current strength as a consequence of magnetic 
saturation. The equation of a circuit which contains such a 
variable element (provided it is a single valued function) may be 
written down in operational notation 

ZJ+<^(J)=i7(0, 

or 

(311) 

In this equation Z is, of course, to be taken as the impedance of 
the invariable part of the circuit, the indicial admittance of which 
is denoted by the usual symbol A(t). 

Equation (311) may be interpreted as the equation of the 
current I{t) in a circuit of invariable impedance Z when subjected 
to an applied e.m.f. consequently, by aid of our 

fundamental formula, I{t) is given by 

The first integral is simply the current in the invariable circuit of 
impedance Z in response to the applied e.m.f. E{t); denoting this 
by I(t), we have 

i(t) =loit) ~jJ^A(t-y)my)]dy. 

This is a functional integral equation, the solution of which is 
gotten by some process of successive approximations. For 
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example, provided the sequence converges, I{t) is the limit as n 
approaches infinity of the sequence 

-^2(0? • ‘ •? 

where the successive terms of the sequence are defined by the 
relations: 

li(0 =/o(0 -j^^\{t-y)4>[h{y)]dy. 


=Io(0 - A it-y)4>{IM]dy. 

We shall not pursue the discussion of non-linear circuits 
further, in view of their mathematical complexity and their 
relatively specialized technical interest. The reader who is 
interested may, however, consult the writer’s paper on Variable 
Electrical Systems,^ for a fuller treatment of the subject. 

^ Phys, Rev. Feb., 1921. 




CHAPTER XI 


THE APPLICATION OF THE FOURIER INTEGRAL TO 
ELECTRIC CIRCUIT THEORY 

The application of Pourier^s series in electrotechnics is a com¬ 
monplacethe use of the Fourier integral, however, has largely 
remained in the hands of professional mathematicians. An 
outstanding distinction between the series and the integral, 
from which the greater power of the latter may be inferred, is 
that the series represents only a periodic regularly recurrent 
function, whereas the integral is capable of representing a non¬ 
periodic function: in fact all types of functions, subject to certain 
mathematical restrictions which are usually satisfied in physical 
problems. 

Before taking up the application of the Fourier Integral to 
Electric Circuit Theory, we shall very briefly review the elemen¬ 
tary mathematics of the series and integral; for a fuller treatment 
the reader is referred to Byerly, Fourier^s^^Series and Spherical 
Harmonics.^ 

Consider a function <^(Q, which in the region 0<t<T is finite, 
single-valued and has only a finite number of discontinuities or 
of maxima or minima. In this region it can then be expressed 
as the Fourier series 

00 

4>(t) = jAo+'^, 

1 

where 

= f <t>it)-QOS 

IJo \J- / (313) 

=I0 (0 • sin t^dt. 

1 In this chapter the Fourier integral is approached from the view-point of 
its physical application and no completeness or rigour is claimed for the 
treatment. The mathematical theory of the Fourier integral is, of course, 
completely developed in treatises on the subject. The object of this chapter 
is merely to outline some of its applications. 
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An COS 




(312) 
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An equivalent series is 

00 

<!>{{) =5^’o+2^” 

where 

F. = VAl+Bi, (316) 

en = tail“K5n/An). 

This expansion is valid in the region 0<i<r, irrespective of 
the form of the function elsewhere. Let us, however, assume 
that the function repeats itself in the period T: that is 

Then the expansion represents the function in the region 
-NT<t<{N+ 1)T. 

Pinally if N is made infinite, the function is truly periodic and 
the Fourier series represents it for all positive and negative 
values of time. 

It follows from the foregoing that, if the Fourier series repre¬ 
sents the function for all positive and negative values of time, 
the function must be periodic for all positive and negative values 
of time; otherwise the expansion is valid only over a restricted 
range of time. 

Now let us suppose that is non-periodic. For convenience, 
in connection with subsequent applications we shall suppose that 
it is zero for all finite negative values of time, that it converges 
to zero as it 00 , and that 

/^co 

\<l>{t)\dt 

exists. Such a function obviously cannot be represented by the 
usual Fourier series for all finite positive and negative values of 
time; it can be represented, however, by the limiting form 
assumed by the series as the fundamental period T is made 
infinite. That is, we can assume that the function is periodic 
in an infinite fundamental period and this will not affect the 
expansion for finite positive and negative values of time. Pro¬ 
ceeding in this way and putting the fundamental period T 
equal to infinity in the limit, the Fourier series (314) becomes an 
infinite integral and we get 

(hU) = 1 r F{cjo) ' cos (wi— d(cjo))dco 


(316) 
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where 


and 


'co) = I <^(i) cos coi di j 

oo ^ CO 

</)(i) sin o)t dt-^ I 4>ii) cos coi dt (318) 
0 Jo 


This is the Fourier integral identity of the function (j){t) and is 
valid for all finite positive and negative values of time. 

In physical applications, particularly those to electric circuit 
theory, it is often convenient to employ exponential instead of 
trigonometric functions. The required transformation follows 
easily from the relation 

e*^ = cos 6+i sin (9, ^ = '\/ —1. 


Thus if we write — m the Fourier series (312) is easily 
reduced to the form 

-j-oo 


<#>(«) 

(319) 

— Oo 

where 


1 n 

F{incao)=-^ 1 

(320) 

In precisely similar manner the Fourier integral (316) 
written as 

can be 

0(«) = j 

(321) 

- ^00 poo 

<^(r)dT 

Stt Jo J - CO 

(322) 


Applications to Electric Circuit Theory,—^Let us assume that 
at time t^—NT, an electromotive force E{t)j periodic in funda¬ 
mental period T, is impressed on a circuit of complex impedance 
Z{io)), where w denotes 27r times the frequency. Required the 
resultant current /. 

For values of t>—NT the electromotive force (see formula 
(319)) can be expressed as the Fourier series 

oO 

E{t) == 
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where 

1 

F(mcoo) == yj E(r)e-^^^ordT, 
The resultant current for t> —NT is therefore 


CO 



Ziinoio) 


transient oscillations 
initiated at time 
-NT. 


If we are concerned with the current for values of i>0, and if 
NT is made sufficiently large, the initial transients will have 
died away and the complete current for t>0, will be given by 


CO 



Z(m(x)Q) 


(323) 


This formula implies the periodic character of E(t) for suffi¬ 
ciently large negative values of time. If, however, E{t) is zero 
for negative values of time, we can employ the Fourier integrals 
(321) and (322) in precisely the same way and get, as the com- 
plete expression or the current for positive or negative values 
of time:— 



(324) 

si 

(325) 


The infinite integrals (324) and (325) formulate the current in the 
network, specified by the impedance function Z(io)), in response to 
an electromotive force E(t) impressed at time t-Q; they therefore 
mathematically formulate, by aid of the Fourier integral identity, 
the fundamental problem dealt with in the preceding chapters and 
solved by aid of the operational calculus. 

No attempt will be made here to discuss the solution of the 
infinite integral (325), which is usually a problem presenting formi¬ 
dable difficulties, even to the professional mathematician. The 
general method of solution is by contour integration in the com¬ 
plex plane and the calculus of residues. By this process it has 
been successfully applied to the solution of special problems, and 
also to deriving some general forms of solution such as the expan- 
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sion theorem solution.^ Compared, however, with the opera¬ 
tional calculus, it has no advantages from the standpoint of rig¬ 
our, and lacks entirely the remarkable simplicity and directness 
of the Heaviside method. 

In the direct solution of circuit problems, therefore, it is 
believed that the application of the Fourier integral is attended 
by few if any advantages, and presents formidable mathematical 
difficulties. On the other hand, there are certain types of prob¬ 
lems encountered in circuit theory, where the Fourier integral 
is a powerful tool. These will be briefly discussed. 

The Indicial Admittance Expressed as a Fourier Integral.— 
Without recourse to the general formulas (324) and (325), it is 
a very simple matter to formulate the indicial admittance A(i) 
(and the Heaviside Function h(t)) as a Fourier Integral, To 
this we now proceed, starting with the known identity 

,>0 (3.6) 

=0 

Evidently this expression represents the unit function, zero 
before, unity after, time 

Now let us suppose that the unit e.m.f,, as formulated by 
(326), is impressed on a network of impedance and let us 

write 

SO that o:(co) and l3(o)) are the two components of the steady state 
admittance at frequency oi/2T, 

Now the first term, of (326) is of zero frequency and pro¬ 
duces a current ^ aj(0). The integral term contains components 
of all frequencies, that lying between co and co + dw, being 

" “ sin cotdo) 

TC 0} 

and the corresponding current is 

2 1 

—sin (co^+^(co)^ico. 

Trcjo\Z{to))\ 

^ Bush, “Summary of Wagner^s Proof of Heaviside^s Formula,’^ Proc, 
Inst, of Radio Engineers. Oct., 1917. Fry. “The Solution of Circuit 
Problems.^^ (Phys. Rev. Aug., 1919.) 
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Consequently the total current^ which is the indicial admittance, 
is given by 


A.{t)=\oL 



1 f“sin (co<+0(co)). 

(328) 

, 1 f “«(co) . , , 

(0) H— 1 —^ sm ut dco 

^Jo “ 

(329) 

,1 f"/3(co) , , 

_|— '' V. cos coi dco. 

^Jo 



The same process applies, obviously, to the Heaviside function 
/i(0, representing current, charge, potential drop, or any of the 
circuit variables with which we may be concerned. 

Thus if the generalized steady state impedance is represented 
by H(ico), and if we write 


we have 

/i(i) = 0^(0)+^ r sin oyt doo (330) 

2 ttJo ^ 

ml . 7 

+ 1 - - COS cot dco, 

TTjo ^ 

Formula (329) for the indicial admittance is then simply a 
special case of formula (330) for the Heaviside function h{t); 
that is, the response of the network to a unit e.m.f, impressed 
at time ^ = 0. 

By (330) we have, changing the sign of t 

h{ -t) = L(0) A f AA sin (331) 

2 TTjo CO 

4-”“ I —" cos cot dco 
ttJo 

By subtraction and addition of (330) and (331) we get at once 
h(t)-h(-t) A sin cof (332) 

TCjo 0) 

2 C^Hco) 

hit)+hi-t)==aiO) + A cos cot (333) 

ttJo 
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Now let us suppose that i>0; then since h{t) is identically 
zero for negative values of it follows that /i(*—/) is zero in equa¬ 
tions (332) and (333), and we have 

/i(^)=" f sin Oil ■ (334) 

ttJo ^ 

==a(0)+" r cos coJ5 (335) 

ttJo ^ 


provided i5>0. It is therefore only necessary to evaluate one of 
the infinite integrals of equation (330). As an interesting corol¬ 
lary, the behavior of the network under all circumstances is completely 
determined if either the real or imaginary component of the complex 
steady state admittance is specified over the entire frequeny range. 
The analytical evaluation of h{t) from (334) or (335) is usually 
quite complicated, unless the infinite integral can be recognized, 
and involves advanced mathematical methods. Its numerical 
evaluation can sometimes be effected without prohibitive labor 
by the following process, applied to either (334) or (335). Select¬ 
ing (334), it can be written as 


and when t becomes very large, this approaches 
aiO),- I ;r sin \dX = a(0) 

A 

which is, of course, the ultimate steady value of h{t). Now 
write (336) in the form 


(336) 


(337) 


h{t) 


IX 


sin X.dX 
sin X.dX 


(338) 


where x may have any finite real value we please. Now having 
assigned some particular value to a;, the first integral of (338) 
can be evaluated mechanically or numerically. The second 
integral can then be expanded asymptotically by repeated 
partial integrations as follows: 
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Repeating this process we get for the second integral term of 
(338) the series 


2 

TT 


cosx 1 — 


2 . 

— sin X 

TT 


dx^ dx^ 
dx dx^ 



(339) 


which may be used for computation for sufficiently large values 
of t and X, 

This process for evaluating (338) is usually laborious but is 
often quite practicable. It will be observed that the computa- 
tation becomes easier for large values of t On the other hand, 
the evaluation of ]i(t) by the operational calculus is usually 
simpler for small values of t. This suggests that in some prob¬ 
lems the numerical solution may often be advantageously 
effected by a combination of the operational calculus and the 
Fourier integral, using the former to calculate h{t) for small 
values of and formula (338) to calculate it for large values of 
t. In practical problems, therefore, the Fourier integral may 
be of great value in supplementing the operational calculus. 

The Building-up of Alternating Currents,—The preceding 
analysis can be generalized and extended to the case of sinusoidal 
impressed forces by employing the expressions 

^1,1 r"^sin \t 1 /o^n\ 

cos o>t I — d\ (340) 

ttJo a J 

and 


sin cot] 


11 f'sinj^i 




(341) 


Formula (340) represents a function which is zero for negative 
values of time i, and equal to cos xt for positive values of time 
t, while (341) represents a function which is zero for ^<0 and 
sin 03 t for t>0. They can therefore be applied to formulating and 
investigating the behavior of networks in response to suddenly 
applied sinusoidal electromotive forces. 

The application suggested by the preceding is of considerable 
practical importance because in long transmission systems the 
calculaton of the building-up of alternating currents is extremely 
complicated and laborious. On the other hand an investigation 
of the mode of building-up is necessary in long communication 
systems since on the character and duration of the building-up 
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process, the excellence of a signal transmission system largely 
depends. 

To formulate the problem of the building-up of alternating 
currents, let us suppose that an e.m.f. E cos ot is suddenly applied 
at reference time ^ = 0, to a net work or transmission system of 
complex steady state impedance 

It is then always possible to write the resultant current in the 
form, 

■^(0=^ {(1+p) cos (o)t-B)+a sin — (342) 

= |V(l+p)'+^^ B(co) - d) (343) 

where 6 = tan*”^::— 

1+p 

In these formulas p, o* and d are time functions which it is our 
problem to evaluate. Evidently p and cr must be *-1 and 0 
respectively for negative values of time and approach the limit 
1 and 0 respectively as 

The function (l+py+a- of equation (343) is proportional 

to the envelope of the oscillations. In an engineering study of 
the building-up process we are often not concerned with the 
instantaneous values of the current but merely with their 
envelope. In such cases an investigation of this function is 
sufficient. 

The mathematical formulation of p and <r as Fourier integrals 
is based on formula (340), Starting with this formula as repre¬ 
senting a suddenly applied sinusoidal e.m.f., cos of frequency 
a)/27r, and proceeding by reasoning strictly analogous to that 
employed in the preceding section, it can be shown that p and 
cr are given by the formulas 

P=^J {P.(X)+P.(-X)}sin t\ (344) 

r“{Qo.(X)-Q.(-X)} cosfX 
cr=Jj^”{Q„(X)+Q„(-X)} sin a 


(345) 
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where 

P.iS) = cos [(-BSTx) -B(co)] (346) 

= |Z(Sh-X)| (347) 

By virtue of the fact that p and <t are identically — 1 and 0 
respectively for negative values of time t, it follows at once (See 
derivation of (332) and (333) from (330)), that we have also for 
'positive values of time (t>0) 

1+p {PM+P.(-^)} sin a (348-a) 

1-p = 1Q.(X)-Q.(~X)} cos fK (348-b) 

and 

9 w 7X 

<r= {Q.(X)+Q.(-X)} siniX (349-a) 

^Jo ^ 

= !P.(X)-P.(~X)} cos a (349~b) 

\Jo ^ 

For important types of transmission systems, particularly 
the periodically loaded line, these formulas have been success¬ 
fully dealt with and solutions of a satisfactory approximate 
character obtained. For further details the reader is referred to 
a paper on ^^The Building-up of Sinusoidal Currents in Long 
Periodically Loaded Lines^’ (Bell System Technical Jouimal, Oct., 
1924). 

The Distortion of Signals in Communication Circuits,—In 
electrical communication, the transmission system must be so 
designed that the received signal must be a more or less faithful 
copy of the transmitted signal; the departure from this require¬ 
ment is termed distortion. The Fourier integral provides a 
simple method of formulating the conditions for distortionless 
transmission, and the factors on which distortion depend. 

Let/(it) be the voltage impres>sed on the sending end of a trans¬ 
mission system of transfer impedance 

Z{io>) - Z\(io>)\e^^^‘^\ 

and let us suppose that f{t) is expressed as the Fourier integral 

/(/;) = ^ r^i?(u)-cos (co<-0(co)) (350) 

^Jo 
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where u) = 2jr times the frequency. The received current is then 
given by 

1(0 =- cos (cot-e(a,)-5(«)). (351) 

Now f(t) of formula (350) contains components of all fre¬ 
quencies from zero to infinity; in other words F{(jo) exists for all 
values of w. Practically, however, P(co) is negligibly small 
except in a range depending on the character of the signal. 
For example in slow speed telegraphy, the essential frequencies 
lie in a band between zero and 10 to 20 cycles per second. In 
telephony the essential frequencies occupy a band extending 
from a low frequency, say 50 c.p.s., to an upper limiting fre¬ 
quency of approximately 4000 c.p.s. 

Let us suppose, therefore, that the essential frequency range 
of f(t) lies between coi/27r and aj2/27r. Then equations (360) 
and (351) can be replaced by 

fit) r”V(a.) cos (a><-5(w)) (352) 

(cot-0(co)-5(to)). (353) 

Now let us suppose that the transmission system is so designed 
that over the range o;i<co<co 2 

\Zii(^)\ = K 

B (co) = COT, 

where K and r are constants. Formula (353) then becomes 

Iii)=^ r>(a,) cos io>it-r)-eim)). (354) 

Comparison of (352) and (354) shows at once that 

I{t)^yit-r). (355) 

That is, the received current has precisely the same wave form 
as the impressed e.m.f.; it is simply changed in amplitude by 
the factor 1/N, and delayed in time by the ^^transmission time^^ r. 
The received signal is therefore a faithful copy of the transmitted 
signal and the transmission is distortionless. We thus arrive 
at the following important proposition. 
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The necessary and sufficient conditions for the distortionless 
transmission of signals in communication systems is thatj over the 
essential range of frequencies contained in the transmitted sig^ial, 
the transfer impedance of the transmission circuit he equalized both 
as regards amplitude and phase; that iSj the amplitude must be 
constant and the phase angle linear in the frequency. 

The case where the amplitude alone is equalized over the 
frequency range corresponding to the spectrum of the signal is 
covered by the following proposition. 

The necessary and sufficient condition that the transmitted and 
received signals have the same relative energy content with respect 
to frequency is that the transmission system be equalized as regards 
amplitude over the range of frequencies contained in the signal. 

It is to be observed that this condition alone, however, does 
not insure preservation of wave form. In types of signalling 
systems, therefore, where preservation of wave form is essential, 
it is necessary to equalize both phase and amplitude. 

The Energy Absorbed from Transient Applied Forces.—In 
many technical problems, the complete solution for the instan¬ 
taneous current due to suddenly applied electromotive forces, 
although formally straight-forward, involves a prohibitive 
amount of labor. In yet others, the applied forces may be 
random and specified only by their mean square values. In 
such problems a great deal of useful information is furnished by 
the mean power and mean square current absorbed by the net¬ 
work, and to the calculation of these quantities, the Fourier 
integral is ideally adapted. Its application depends on the 
following proposition, due to Rayleigh (Phil. Mag., Vol. 27, 
1889, p. 466), and its corollary. 

Let a function c/)(i), supposed to exist only in the epoch 0</< jT, 
be formulated as the Fourier integral 


f |/(c*j)| . cos [o)t^d{cv)]do) 
^Jo 

where 

fM 


■-i[r 


(/)(/) COS o:t dt\ +1^ j (h{t) sin ot dt 
<j){t) cos o)t df. 


2l I 


rr rT 

tan 0(co)= I </)(0sm oi dt-^ \ ' 

./o Jo 

m'rdw, 


Then 
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whereby the time integral is transformed into an integral with 
respect to frequency. 

A corollary of this theorem is as follows: 

If two functions supposed to exist only in the epoch 

0<t<Tj are formulated by the Fourier integrals 

' |/i(co)|*COS 

^Jo 

02(0='" I |/2(w)l-COSM~“^2(w)]dc0, 

^Jo 

then 

f f |/i(“)i-|/2(w)|‘ cos (0i-02)dw. 

Jo '^Jo 

The application of these theorems to circuit theory proceeds 
as follows:— 

If an electromotive force supposed to exist only in the 
epoch 0<t<T, is applied to a network of complex impedance 
Z(fco) = we laiow from the preceding discussion of the 

Fourier integral^ that the electromotive force JS(t) and current 
I(t) are expressible as the Fourier integrals 


f |/(aj)|- cos (o)t—6(a)))do), 
^Jo 

1(0=^ f 

^Jo 


(356) 


COS (cot^d(o:)—l3(o)))doi}. 


It follows at once from Rayleigh’s theorem that 


^ OQ 1 Z' ° 

Pdt=- 
Jo ^Jo 


\Z{m)\ 


idco. 


(357) 


Now let In be the current absorbed in branch n; let z(ico) — 
| 2 (fco)|e^“^“^ be the impedance of that branch and let En{t) be the 
potential drop across that branch. It follows at once from the 
corollary to Rayleigh’s theorem that 

W= (358) 

Jo ^Jo 


Formulas (357) and (358) formulate the mean square current 
and mean power absorbed by the branch of the network under 
consideration, and enable us to calculate these quantities, even 
in the case of complicated networks, with a minimum of labor. 
Formula (357) is particularly well adapted to computation 
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because the integrand is everywhere positive, pennitting, in 
most problems, of easy numerical integration, whereas the 
analytical solution may be complicated. 

Selective Circuits and Interference.—An important applica¬ 
tion of the preceding formulas is to the theory of selective circuits 
with particular reference to the protection they afford against 
irregular and transient disturbances such as ‘^noise” in telephone 
systems and '^static interference’^ in radio systems. In these 
cases the disturbance or interference is essentially random in 
character; its instantaneous wave form is unknown, and it can 
be specified only by its ‘^frequency distribution.” Under these 
circumstances the problem is a statistical one and the only 
information deducible is the mean power and mean square current 
absorbed by the circuit from the disturbance. This is calculable 
from the preceding formulas when the statistical or average 
frequency distribution of the disturbance is determined or 
specified. A brief sketch of the analytical theory follows. 

Let a disturbance /(O, for convenience supposed to exist only 
in the epoch 0<i<r, be impressed on a network of impedance 
Z(zco). If the disturbance is specified for all values 'of t in this 
epoch, it is theoretically possible to calculate jP(co) and 0{oT) and 
to express fit) as the Fourier integral 

/(O = ” rT(«). cos («<-0(a,))dw. (359) 

'^Jo 

In any case even if insufficient information is supplied to calcu¬ 
late F((xi) and 6{(x) we know that/(<) is expressible as an integral 
of this form. Prom formula (358) we have 

I FT I 

Now let the epoch T become indefinitely large; if the disturbance 
is random in character, the left hand side of this equation will 
approach a limit; equal to the mean square value of /(i), and we 
can write 

= ~ rT(a))d£0 (360) 

^Jo 

where i2(w) is the limit of (w) as T becomes very large. 

li(co) will be termed the staiistical frequency spectrum of the 
disturbance. 
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Corresponding to the above the mean square current absorbed 
by the network is given by (See equ* (357)) 

72 _1 r n 


Formula (361) suggests a method for the determination of the 
frequency spectrum of the disturbance as follows: Let us suppose 
that the network is a selective circuit which absorbs current only 
in a narrow frequency band corresponding to a)i<co<co 2 . Then 
if 0 ^ 2 —coi is made sufficiently small^ formula (361) becomes 



where is defined by the equation 

_ co.dw do) 

"" Jo ‘ Jo 


(362) 


(363) 


It follows from the preceding that measurement of the mean 
square current and calculation or measurement of the impedance 
characteristics of the selective circuit supplies the information 
necessary to compute the statistical frequency spectrum of the 
disturbance. 

Returning to formula (361) let us suppose that the selective 
circuit is designed to receive a band of frequencies a)i<co<w 2 
and excludes all other frequencies. The mean square current 
is then given by 



where wi<co'<wo. 

If therefore the selective circuit is required to receive a band 
of frequencies corresponding to coi<oj<co 2 , formulas (364) and 
(365) give the irreducible minimum of interference. That is, the 
selective circuit must necessarily receive the energy of the dis¬ 
turbance which lies in the frequency range it is designed to receive. 
A definite limit is therefore set to the protection against inter¬ 
ference which can be obtained by means of selective circuits. 

This limitation is inherent in the nature of the interference 
and the frequency requirements of the signalling system. 
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For further and more detailed application of formulas (357) 
and (358) the reader is referred to the following papers: 

Transient Oscillations in Electric Wave Filters, Bell System 
Technical Journal, July, 1923. 

Selective Circuits and Static Interference, Trans. A. I. E. E.^ 
1924. 

An Application of the Periodogram to Wireless (Burch and 
Bloemsma), Phil. Mag., Feb., 1925. 

The Theory of the Schrotteffekt (Fry), Journal Franklin 
Institute, Feb., 1925. 

The foregoing must conclude our very brief account of the 
Fourier Integral and its applications in Electric Circuit theory: 
an adequate treatment of this subject would require a treatise 
in itself, and is beyond the scope of the present work. All that 
has been attempted is to give a very brief introduction to its 
significance in physical problems and a few of its outstanding 
applications in circuit theoiy. The reader who is interested in 
pursuing this subject further is referred to a paper by T. C. Fry 
on “The Solution of Circuit Problems’^ (Phys. Rev., Aug., 
1919), which gives a rigorous discussion of the solution of the 
Fourier Integral by contour integration, together with some 
general fornxs of solution of the circuit problems. 
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INDEX 


A 

Activity, equation of, 5, 6 
Alternating currents, building up of, 
77 

application of Fourier integral to 
theory of, 181-183 
Artificial cable, indicial admittance 
of, 124 
line, 117-129 

characteristic impedance of, 119 
definition of, 117 
indicial admittance of, 124 
operational equations of, 118 
propagation constant of, 118 
Asymptotic expansion, 29 
scries, definition of, 54 
solution of operational equations, 
62-84 

generalization of Heaviside rule, 
78-84 
rule 1, 64 
example, 65, 67 
theory of, 63, 64 
rule, II, 67, 68 
examples, 70-74 
theory of, 68-70 

B 

Bessel functions, 129-131 
Borers theorem, 41 
Boundary conditions, 11 

C 

Cable, non-inductive, 85-99 
building up of alternating cur¬ 
rents in ,96-99 

current and voltage waves in 
infinite cable, 87-89 


Cable, non-inductive, differential 
equations of, 50-53 
leakage in, theory and formulas, 
93-96 

operational equations for current 
and voltage, 85 

solution of operational equations 
in infinite cable, 86, 87 
Circuit elements, definition of, 2, 3 
Compiementarj^ solution, 10, 11 

D 

Determinant of network, 11 
Distortion of signals, theory of, 183- 
185 

Dynamic theory of electric circuits, 2 
E 

Electric circuit theory, definition of, 

1 

differential equations of, 3, 4, 5 
fundamental formulas of, 16,21, 22 
fundamental problems of, 13, 14 
Electric network, definition of, 1 
Energy absorbed from transient 
forces, 185-187 
Equation of activity, 5, 6 
Expansion theorem solution, 32-37 
Exponential solution, 6, 7, 8, 9 

F 

Finite transmission line, 132-155 
expansion solution of operational 
equations, 152-155 
open-circuited line, 135 
operational equations for current 
and voltage, 132, 133 
for indicial admittance, 133 
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ELECTRIC CIRCUIT THEORY 


Finite transmission line, short- 
circuited line, 135 
wave solution of operational equa¬ 
tions, 133-136 

Fourier integral, application to elec¬ 
tric circuit theory, 176-178 
mathematical theory of, 17^176 

H 

Heaviside function, 22 
rule, for asymptotic solutions, 62- 
84 

generalization of, 78-84 
I 

Impedance function, 10, 18, 19, 21 
generalized, 22 

Indicial admittance, definition of, 14 
expressed as a Fourier integral, 
178-181 

integral equation of, IS, 19 
reciprocal property of, 16 
Induction generator, theory of, 165- 
167 

Infinite integrals, table of, 39-40 
Integral equations, 156-157 
of indicial admittance, 18, 19 
of Laplace, 19 
Interference, 187-188 

K 

Kirchhofi^s laws, 2 

law of cable telegraphy, 
89-92 

generalizations of, 92 
L 

Laplace integral equation, 19 
Loaded lines, 117, 128 

N 

Non-inductive cable, differential 
equations of, 50, 53 
Non-linear circuits, 172-173 
Normal vibrations, 32, 34 


0 

Operational calculus, fundamental 
theorem of, 23 
theorems of, 42-50 
Operational equation, 22, 23 
asymptotic solution of, 62-84 
expansion theorem, solution of, 
32-37 

power series solution of, 21-32 
P 

Partial fraction expansion, 33, 34 
Poisson integral equation, 142-151 
application to electric circuit 
theory, 142 
to finite line, 146, 147 
approximate solution of, 144, 149, 
150 

numerical solution of, 145, 151 
Power senes solution, 21-32 

R 

Rayleigh^s theorem, 185 
Reciprocal property of mutual ele¬ 
ments, 8 
theorem, 8 

S 

Selective circuits, theory of, 187-188 
Steady-state oscillations in variable 
electric networks, 167-172 
solution, 9, 10 

Sudden open circuit, theory of, 160, 
161 

short circuit, theory of, 159, 160 
T 

Table of infinite integrals, 39-40 
Terminal impedances, 132 
Transient distortion, asymptotic for¬ 
mulas for, 77, 78 
definition of, 18 
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Transmission line, 100-116 
attenuation factor, 106 
building up of alternating cur¬ 
rents in, 114, 115 
current and voltage waves in in¬ 
finite line, lOS-113 
differential equations of, 59, 100 
distortionless, 108 
non-dissipative, 106, 107 
operational equations for current 
and voltage, 101 

solution of operational equations 
for infinite line, 102-105 
velocity of phase propagation in, 
lOG 

of propagation in, 105 


V 

Variable circuit element, 161 
electric circuits, 159-173 
networks, mathematical theory 
of, 161-165 

steady-state oscillations in, 
167 t-172 

W 

Wave filters, 120-128 
building up of alternating currents 
in, 128 

current waves in, 125-128 
indicial admittance of low pass, 
121, 122 



